





CAL 


UNIVERSITY 
OF MICHIGAN 


APR 20 1954 


ENGINEERING 
LIBRARY 








| Votume 21 APRIL, 1954 NuMBER 4 





CONTENTS 


The Ultimate Strength of Aluminum-Alloy Formed Structural Shapes in Compression. 
Rosert A. NEEDHAM 


Linearized Theory of the Oscillating Airfoil in Compressible Subsonic Flow 
R. TIMMAN 


Calculation of Compressible Cascade Flow by the Method of Flux Analysis 
SHIGEO UCHIDA 


A Vector Method Approach to the Analysis of the Dynamic Lateral Stability of Aircraft. 
L. STERNFIELD 


Results of Some Base Pressure Experiments at Intermediate Reynolds Numbers with 
L. L. KAVANAU 


The Flow over a Flat Plate with a Small Angle of Attack at Mach Number 1 
GOTTFRIED GUDERLEY 


Readers’ Forum: ‘An Exact Solution for Two-Dimensional Linear Panel Flutter at Supersonic 
Speeds,” by Martin Goland and Yudell L. Luke; “Nose Shapes for Minimum Pressure Drag at 
Supersonic Mach Numbers,” by Leland H. Jorgensen; ‘Streamline Curvature Effects for Finite 
Wings in Wind Tunnels,” by Alan Pope; ‘“Wind-Tunnel Model Position for Transonic Testing,” 
by Alan Pope; “Similar Solutions of Compressible Laminar Boundary-Layer Equations,” by 
Clarence B. Cohen; “The Method of Superposition of Planar Wave Systems,” by Wallace D. 
Hayes; ‘Comments On ‘Ripple-Type Buckling of Sandwich Columns,’” by Paul Seide; “Spin 
Effect on Base Pressure of Cone-Cylinders,”” by R. Lehnert and S. M. Hastings; ‘Buckling Loads 
for Columns of Variable Cross Section,” by B. E. Gatewood 


Published by 
INSTITUTE OF THE AERONAUTICAL SCIENCES, INC. 











JOURNAL OF THE 
AERONAUTICAL SCIENCES 





Associate Editor: Ropert R. DEXTER 


Aerodynamics 


Francis Clauser 
Howard Emmons 
Antonio Ferri 
Alexander Flax 

R. Paul Harrington 
R. T. Jones 

" Carl Kaplan 

H, W. Liepmann 
C. C, Lin 

H, F. Ludloff 
Harold Luskin 
John R. Markham 
Clark B, Millikan 
Shatswell Ober 

W. Bailey Oswald 
Allen E. Puckett 
Elliott G. Reid 

H. J. E. Reid 
Russell Robinson 
G. B. Schubauer 
W. R, Sears 

John Stack 
Theodore Theodorsen 
Theodore von Karman 


Air Transport 


Charles Froesch 
J. A. Herlihy 

R. D. Kelly 
Otto E. Kirchner 
Jerome Lederer 
John C. Leslie 
R. Dixon Speas 


Flight Testing 


W. K. Ebel 

R. R, Gilruth 
Melvin N. Gough 
E, G, Stout 


Editor: Hucu L. Drypen, Director, NACA 


Managing Editor: Berneice H. Jarcx 





Editorial Committee 


Structures and Materials 


W. B. Bleakney 
E, W. Conlon 
L. H. Donnell 
E, C. Hartmann 
N. J. Hoff 
Samuel Levy 

E, E. Lundquist 
Alfred S. Niles 
W. Ramberg 

F. R, Shanley 
H. W. Sibert 

C. R. Strang 

S. Timoshenko 
H. S. Tsien 
Chi-Teh Wang 


Rotating Wing Aircraft 


F. B. Gustafson 

K. H. Hohenemser 

Bartram Kelley 
B. Maloy 


R, 
Ralph McClarren 


R. H, Miller 
R. H. Prewitt 
Igor I. Sikorsky 


Vibration and Flutter 


Lee Arnold 
M. A. Biot 
R. L. Bisplinghoff 
William Bollay 
Chieh-Chien Chang 
I. P, Den Hartog 

. E. Garrick 
Martin Goland 
Theodore Theodorsen 


Aircraft Design 


W. K. Ebel 

Hall L. Hibbard 
Richard Hutton 
C. L. Johnson 

A, A, Kartveli 
C. J. McCarthy 
Fred E. Weick 
Robert J. Woods 


Electronics in Aviation 


E, F, Herzog 
David S. Little 
D..K. Martin 
George Rappaport 
Paul Rosenberg 
L. M. Sherer 


Physiologic Problems 


H. G. Armstrong 

Louis H. Bauer 

Robert J. Benford 
Eugene Du Bois 

W. Randolph Lovelace, II 
Ross A. MacFarland 


Fuels and Oils 


D. P. Barnard 
J. H. Doolittle 
Graham Edgar 
R. T. Goodwin 
S. D. Heron 

A. M. Rothrock 





SUBSCRIPTION RATES 


Journac or THE Asronauticat Sciences, Published Monthly.—United States and Possessions: 1 Year, $12.00; Single Copies, 
$1.50. Foreign Countries Including Canada (American Currency Rates): 1 Year, $13.00; Single Copies, $1.50. 

AERonavTicaL Encinsarinc Review, Published Monthly.—United States and Possessions: 1 Year, $3.00; Single Copies, $0.50. 
Foreign Countries Including Canada (American Currency Rates): 1 Year, $3.50; Single Copies, $0.50. 

Notices of change of address should be sent to the Institute at least 30 days prior to actual change of address, 

Manuscripts for publication, proofs, and all correspondence should be addressed to the Editorial Office of the Journat, 
Correspondence regarding subscriptions may be addressed to Publication Office, 20th and Northampton Sts., Easton, Pa., or to 
the Editorial Office of the JournaL or THE AERONAUTICAL Sciences, 2 East 64th Street, New York 21, New York. 


Copyright, 1954, by the Institute of the Aeronautical Sciences, Inc 


Flight Propulsion 


George W. Brady 
Frank W. Caldwell 
L. S. Hobbs 
E, E. Johnson 
ee Keenan 

. P. Kroom 
Hans Reissner 
Abe Silverstein 
C, Fayette Taylor 
E. S, Taylor 
P. Taylor 
E, S, Thompson 
H. S. Tsien 


Instruments 


Allan G. Binnie 
J. R. Bradburn 
Charles H. Colvin 
C. S. Draper 

S. G, Eskin 

O. E. Esval 

H, W. Holzapfel 
W. L. Howland 


. Hugo Schuck 
. L. Seward, Jr. 
D. Statham 
C. Sylvander 
K. Warner 


Meteorology 


C. F. Brooks 
D. M. Little 
E. J. Minser 
F, W. Reichelderfer 
H. C, Willett 


Entered as Second Class Matter at the Post Office, Easton, Pa., May 1, 1937. Acceptance for mailing at a special rate of 
Postage provided for in the Act of August 24,1912. Authorized April 29, 1937. 








An 
stress 
preset 
that a 


variot 


ing 
is obt: 
gle sec 


struct 
Em 

basic 

mately 





JOURNAL OF THE 
ERONAUTICAL SCIENCES 





—_—_— 


VOLUME 21 


APRIL, 1954 


NUMBER 4 





The Ultimate Strength of Aluminum-Alloy 
Formed Structural Shapes in Compression’ 


ROBERT A. NEEDHAM? 


Unwersity of California at Los Angeles 


ABSTRACT 


Anew method of predicting the maximum average (crippling) 
stress of formed sheet-metal structural shapes in compression is 
presented. The proposed method is based on the assumption 
that a formed structural shape consisting of a series of flat plate 
elements can be treated as a series of angle sections possessing 
various degrees of edge support parallel to the direction of load 
ing. A weighted crippling stress of the formed structural shape 
is obtained by summing the crippling loads of the individual an 
gle sections and dividing this summation by the gross area of the 
structural shape. 

Empirical equations for predicting the crippling stress of the 
basic angle sections were obtained from the analysis of approxi 
mately 200 tests of formed angles and channels of 24S-T3 and 
75S-T6 aluminum-alloy sheet. The proposed method of anal 
ysis was applied to 24 typical aircraft The 
agreement of the predicted and test crippling stresses was found 


stringer shapes 
to be very good. 


NOTATION 


width of channel or angle flange, in. 

area, sq.in. 

width of flat plate or half width of channel web, in 
equivalent width, (a + 6)/2, in. 

effective width, in. 

crippling coefficient for angle section 

Young’s modulus, Ibs. per sq.in 

stress at junction of flat plate elements, Ibs. per sq.in 
crippling or maximum average stress, Ibs. per sq.in 
critical plate buckling stress, lbs. per sq.in 
compression yield stress, Ibs. per sq.in. 

plate buckling coefficient 

bend radius, in. 


sheet thickness, in 
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* A thesis submitted in partial satisfaction of the requirements 
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effective modulus ratio, dimensionless 


Poisson's ratio, dimensionless 


INTRODUCTION 


Soho PRIMARY STRUCTURE of the modern airplane 
usually contains many wide, thin plate elements 
loaded in compression. Thin-walled round, square, 
and rectangular tubes and flat or curved sheet to which 
are riveted formed or extruded structural shapes are 
but a few of the examples of this type of construction. 
These members are analyzed as columns. 

As the length of a column of compact cross section 
decreases, the stress at primary failure increases, and 
the failure stress is accurately predicted by the Euler 
equation employing the tangent modulus.” If, how- 
ever, the cross section of the column is made up of thin 
plate elements, the stress at failure will be modified by 
localized plate buckling of these elements. When the 
column slenderness ratio is large, the stress for general 
column failure is less than that at which plate buckling 
takes place, and the Euler equation predicts correctly 
the failure stress. As the slenderness ratio decreases, 
a point is reached at which the stress for general col- 
umn failure and the stress at which a plate element of 
the unstable will be 
equal. 
umn curve begins to depart from the theoretical col 


column cross section becomes 


It is at this point that the experimental col 


Failure above this point is an interaction 
In the short column range 


umn curve. 
of general and local failure. 
(slenderness ratios in the order of 20), failure is due 
almost entirely to the plate buckling phenomenon 
i.e., the column strength is not materially affected by 
length. The stress in this region is the upper allow 
able of the column stress and is variously referred to as 
the crippling, crushing, or maximum average stress. 


217 
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PLATE BUCKLING 


When flat plates are subjected to edge compression, 
the critical stress at which the plate becomes unstable 
and buckles is expressed by the well-known relation- 
ship 


F../n = kw*Et?/12(1 — p?)b? (1) 


The nondimensional coefficient, 7, takes into account 
the reduction of the modulus of elasticity for stresses 
above the elastic range. Methods of evaluating 7 are 
discussed in reference 1. Eq. (1) may be used to pre- 
dict the stress at which a flat plate will buckle, provided 
the degree of edge fixity of the plate is known. Thin- 
walled equal angles and square tubes are examples of 
structural shapes with known boundary conditions. 
Thin-walled zee and channel sections are also examples 
of structural shapes for which the buckling coefficients 
are known. In these cases, however, the buckling 
coefficients vary with the ratio of the flange width to 
the back width, and Eq. (1) is written as 

Ct pera = ky’ E(t/b)? (2) 

n 1201 — u2)b? 
Values of k, are obtained from charts such as shown in 
Fig. 1.? 

The National Advisory Committee for Aeronautics 
has conducted extensive investigations in order to de- 
termine plate compressive strengths of a number of 
aircraft structural materials. The results of these 
tests are summarized in reference 3. An obvious and 





i i as 
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FIG.| ~- Kw FOR BUCKLING OF Z-— AND C -SECTIONS 


highly significant fact brought out by these tests was 
that the maximum stress developed by the sections 
tested was always higher than the stress at which plate 
buckling occurred. In the discussion that follows, the 
maximum average stress developed by a structural] 
shape will be referred to as the crippling stress. Exist. 
ing methods of predicting this stress will be discussed, 
and, finally, a new method of crippling stress prediction, 
based on the results of over 200 tests, will be presented, 


EXISTING METHODS OF PREDICTING CRIPPLING 
STRESSES 


farly in the development of thin-sheet structural 
analysis techniques, a method of crippling stress pre- 
diction proposed by several authors was that of sub. 
dividing a given cross section into a series of flat plate 
elements. The elements were assumed to possess vari- 
ous degrees of edge support, usually hinged. The buck- 
ling load for each element was calculated from Eq. (| 
(n = 1.0) and the area of the element. The crippling 
stress of the section was taken as the summation of the 
buckling loads of the individual elements divided by 
the total area of the section. (See, for example, refer 
ence 4.) There are two obvious discrepancies in this 
technique. After buckling of contiguous flat elements 
of a structural shape, narrow strips along the junc- 
tions of the plate elements continue to support load 
and, actually, may permit the section to carry load 
several times greater than the theoretical buckling load 
of the strongest flat plate element of the cross section. 
Further, the assumption that adjacent elements pro- 
vide mutual support equivalent to simply supported 
(hinged) edges is often overconservative. 

A method of crippling stress prediction that has found 
wide usage in the aircraft industry was proposed by 
Crockett® in 1942. The method is empirical and is 
based on design curves that are trial-and-error vari- 
ations of the theoretical to give the best possible agree- 
ment with test data (see reference 5, page 501). Crock- 
ett’s method of crippling stress prediction is somewhat 
the same as the technique outlined above—that 1s, 
the method consists of summing the crippling loads of 
the components and dividing by the section area to 
obtain the crippling stress. In the Crockett method, 
however, the crippling stress of curved elements, 
formed at the junctions of the flat plate elements, 1s 
accounted for. The principal disadvantage of this 
method lies in the fact that a considerable amount oi 
computation time is involved in analyzing complex 
shapes. 

Schuette® has made a comprehensive study of the 
results of the NACA channel- and zee-section crippling 
strength tests. According to his findings, the crip- 
pling stress of formed sheet-metal sections of 
and 24S-T aluminum alloys may be expressed by 


Fee] Fer = K(Fer/Fey)~°™ 


17S-T 
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Eq. (3) was derived from logarithmic plots of F,, Po 
against F.,/F.,. Here, as before, the critical stress, 
F.,, is the stress at the start of plate buckling and may 
be calculated from Eq. (1). The coefficient, A, was 
introduced to account for the fact that the ratio F,./ F-, 
isalso a function of the b/t ratio of the section (Fig. 1). 


The variation of AK with }/t is given in Table I. 


TABLE I 
(From Reference 6, Page 6) 
100 
0.63 


b/t 20 30 40 50 
K 0.90 0.82 0.75 0.71 


The values of A given in this table apply to plates 
formed at the edges to an inside bend radius of three 
times the sheet thickness (see reference 6, page 6). 

Schuette’s equations were used by The Dow Chemical 
Company in preparing a set of crippling stress curves 
for magnesium-alloy formed structural shapes.’ 
CRIPPLING STRESS PREDICTION 


A New MetTHOD OF 


The various methods of crippling stress prediction 
just discussed are all based on the assumption that a 
structural shape can be treated as a series of isotropic 
flat plate elements, the strengths of which can be pre- 
dicted by the classical buckling equations. It is well 
established, however, that substantial increases in the 
compression yield stress occur in the vicinity of bend 


lines of formed shapes. In some aluminum alloys this 


increase may amount to as much as 25 per cent over the 
yield stress of the unformed sheet.* A flat element of a 
formed structural shape then becomes, in effect, a 
plate possessing variable mechanical properties normal 
to the direction of loading, and there are no provisions 
in the plate buckling equations for this effect. 

It occurred to the writer that a more rational ap- 
proach to the problem of predicting crippling stresses 
could be effected by considering the basic structural 
element as an angle section rather than a plate. Fur- 
ther, it was believed that the strength of these angle 
sections could be established experimentally in such a 
way that the local increase in yield stress at the bend, 
brought about by cold-forming the shape, could be 
accounted for. Observations of tests of formed sheet- 
metal shapes in compression indicate that initial in- 
stability starts with the flat plate elements and that 
complete failure of the section occurs when the buckle 
pattern passes through the bends. Hence, if the 
edges of the flanges of the basic angle section are taken 
to the centers of contiguous flat plate elements (see, 
for example, Fig. 5b), the ‘‘cut’’ is made at the least 
sensitive region of the plate. 

A series of tests was conducted on formed 24S-T3 
plain (unclad), and 24S-T3 and 75S-T6 clad, alumi 
num-alloy equal and unequal angles and channels as 
part of a test program to establish crippling stress 
equations. Because of the large number of tests con- 
ducted (over 200), tabulated results have not been 


included in this paper. All specimens were formed 
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on a leaf brake and had inside bend radii of approxi- 
mately three times the sheet thickness. The specimens 
were cut to a length corresponding to a_ pin-ended 
slenderness ratio of approximately 20. The ends of the 
specimens were ground flat and parallel, and areas were 
calculated from the weights and lengths of the speci- 
mens using a material density of 0.100 Ib. per cu.in. 

The tests of this series were made in a Tinius-Olsen 
hydraulic testing machine having a maximum capacity 
of 120,000 Ibs. and an intermediate load range of 20,000 
Ibs. The machine had been calibrated just prior to 
the start of the test program, and it is believed that the 
errors in the load readings were less than +1 per cent 
of the load range used. The specimens were tested 
flat-ended with the ends of the sections bearing di- 
rectly on the guided platens of the machine. Maxi- 
mum care was exercised in obtaining uniform contact 
of the ends of the sections and the testing machine 
platens, and in some instances it was necessary to re 
grind the ends of the specimens to ensure good contact. 
The maximum load reached in the test was recorded, 
and this load divided by the gross area of the section 
was taken as the experimental crippling stress. 

The first group of sections tested consisted of 78 equal 
Fee, 
The 


angles. Logarithmic plots of crippling stress, 


against the parameter }/t are shown in Fig. 2. 


solid lines on the plot represent least square fits of 
equations of the form 
a =k (b $)®.% (4 
where k, the stress intercept, is as follows: 
24S-T3 clad: k 
24S-T3 plain: k 
795-T6 clad: k = 


200,100 Ibs. per sq.in. 


225,000 Ibs. per sq.in. 
259,000 Ibs. per sq.in. 


In Eq. (4) the flange width } is measured to the center- 
line of the adjacent web. The experimental scatter 
appears to be random and is characteristic of com 
pression tests of thin-sheet sections. 

The majority of structural shapes that are of interest 
are made up of angle sections with unequal leg widths. 
Here, a suitable buckling parameter may be defined by 


b’/t = (a + b)/2t 


Thus, 6’/t will be referred to as the ‘‘equivalent” 4! 
and is simply the average flange width, (a + ))/2, 
divided by the sheet thickness. Fig. 3 is a plot of the 
tests of 24 unequal angles of 24S-T3 clad aluminum. 
Here, the abscissa is the equivalent /t of the sections 
as defined above. Included also is the curve for 245- 
T3 clad equal angles of Fig. 2. The agreement of the 


two sets of data is good and definitely establishes 
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ULTIMATE STRENGTH OF 
(a + 6)/2t as a suitable parameter for the buckling 
of unequal angles. 

The excellent agreement of the two sets of tests just 
discussed immediately pointed to the possibility of 
treating a symmetrically formed channel as two un- 
equal angles supported at the midpoint of the web. 
The effective 6/t of the section would then be, as in 
the case of the unequal angle, (a + 6)/2t. Here, the 
dimension 0 is the half-width of the web of the channel. 
Fig. 4 includes logarithmic plots of the channel tests 
for the three materials investigated. In all three cases 
the test data yield equations analogous to Eq. (4), the 
only difference being in the constants k, which, in the 
case of the channels, are higher than the corresponding 
k-values for the angles. This is obviously due to 
the fact that the ‘‘channel angles’’ were restrained along 
one unloaded edge, whereas the angle sections were 
as follows: 


not. The k-values for the channels are 


221,300 Ibs. 
244,000 Ibs. 
272,000 Ibs. 


24S-T3 clad: k = 
248-T3 plain: k = 
79S-T6 clad: k 


per sq.in. 
per sq.in. 
per sq.in. 


The above results are based on the results of 127 tests. 
The 245-T3 and 75S-T6 clad specimens (98 tests) were 
conducted by the writer, while the 24S-T3 plain channel! 
test data were taken from reference 9. The range of 
sheet gages covered in both the angle and channel tests 
is listed in Table IT. 


FORMED 


bo 
to 
pa 


STRUCTURAL SHAPES 


TABLE II 
Material Shape Sheet Gages Tested 
248-T3 Equal angles 0.016, 0.020, 0.025, 0.032, 0.040, 
Clad 0.051 
Unequalangles 0.032, 0.040, 0.051 
Channels 0.016, 0.020, 0.082, 0.040, 0.051 
248-T3 Equal angles 0.025, 0.032, 0.040, 0.051 
Plain Channels 0.102 
758-T6 Equal angles 0.040, 0.051, 0.072 
Clad Channels 0.025, 0.040, 0.051 


Formed sheet-metal structural shapes appear in a 


variety of forms such as zee’s, lipped channels, and 
hats, and presenting empirical crippling stress curves 
for such a variety of structural sections is a practical 
impossibility. It is interesting to note, however, that 
many of these shapes can be considered as being made 
up of a series of unequal angle sections with various 
for example, the 


degrees of edge support. Consider, 


unbalanced zee-section shown in Fig. 5b. This shape 


is made up of three unequal angles. Elements | and 3 
are unequal angles that are supported on one unloaded 
edge and free on the other, while element 2 is an un- 
equal angle that is supported on both unloaded edges. 
(Included also in Fig. 5 is the method of subdividing 
the shape as proposed by Crockett in reference 5.) If 
it is assumed that the crippling load of the entire shape 
is equal to the summation of the crippling loads of the 
three angles, the weighted crippling stress will be given 
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by 
F > (crippling loads of angles) 
sl > (area of angles) 

The above discussion assumes that equations are 
available to predict the crippling strength of unequal 
angles possessing support along both unloaded edges 
(no edge free). Empirical equations for this condi- 
tion of support could be obtained from tests. Such a 
procedure is not necessary, however, if the rational as- 
sumption is made that the percentage increase in crip- 
pling strength of the no-edge-free angle over the one- 
edge-free angle is the same as the increase of the one- 
edge-free angle over the angle with both unloaded 
edges free. For example, consider the angle and 
channel tests of the 24S-T3 clad alloy of Figs. 2 and 4. 
Thus, 


Rno edge free = 221,300 + (221,300 — 200, 100) 
= 242,500 Ibs. per sq.in. 


The stress intercept values for the three materials in- 
vestigated in this paper are listed in Table IIT. 


TABLE III 
Stress Intercept (&) Values for F.. = k/(b’/t)®-75 


k 


Material (Two Edges Free) (One Edge Free) (No Edge Free) 


248-T3 200 , 100 221 ,300 242,500 
Clad 

248-T3 225,000 244,000 263 ,000 
Plain 

75S-T6 259 , 000 272 ,000 285 ,000 
Clad 


DISCUSSION 


The fact that the experimental plots of F,. against 
b'/t of Figs. 2, 3, and 4 all follow a 0.75 power law is 
intriguing. One possible explanation for this may be 
found by considering the sections in the light of effective 
width theories. 

According to von Karman,'® the major portion of the 
load supported by a thin plate in edge compression is 
carried by narrow strips along the edges parallel to the 
direction of loading of the plate. These ‘‘effective 
widths” of sheet, acting at a uniform edge stress, f,,, take 
the place of the total sheet width in considering the 
ultimate strength of the plate. The original equation 
suggested by von Karman for calculating these effective 
widths is 


b, = RLV E/ fre (5) 


The constant, k, in Eq. (5) depends upon the degree 
of edge support along the unloaded edges of the plate. 
If this equation is applied to the equal-angle section of 
Fig. 6 and if it is assumed that the edge at the junction 
of the two legs of the angle is hinged, k = 0.60, and the 
total load supported by the angle is 
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FIG.8 -CRIPPLING OF FS-Ih MAG. ALLOY ZEE'S 


Poe = fuc2bet = 1.20fecl?WE/ fe 


and since the area of the angle is 2b/, the average crip- 


pling stress is 
Fee = Pee/2bt = 0.60(t/b)V fF 


or 
F.. = k,/(6/t) (6) 


where k, = 0.60 V faelé. In reference (11) Marguerre 


proposes the equation 
3 


2b, = b WV Fer! fu (7) 


as the expression for the total effective width. Here, 
F., is the critical edge stress of the plate. If it is again 
assumed that the edges of the plate elements of the 
angle of Fig. 6 are simply supported, F,, = 0.4#(t/b)?, 
and the crippling stress is given by 
F.. = ke/(b/t)*” (8) 
3 - 
where ky = 0.37 V fk. Thus Eqs. (6) and (8), based 
on the effective width equations of von Karman and 
Marguerre, yield exponents for the )/¢ ratio that 
bracket the experimentally derived value of 0.75. 
Eq. (3) of this paper is an expression for the crippling 
strength of zee and channel sections derived by Schuette 
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Material Fey E 
24S-T3 Clad sheet ¢ < 0.051 $2,000 10,000,000 
24S-T3 Plain (unclad) 45,100 10,500,000 
75S-T6 Clad sheet ¢ > 0.032 66,200 10,000 ,000 


from NACA test data. The coefficient A is a function 
of the 6/t ratio of the section (here, 6 is the width of the 
web of the zee or channel section). Values of A for 
various b/t ratios are given in Table I. The variation 
of K with b// given in this table can be expressed 
closely by the equation 

K = 


1.90/(b/t)®-* (9) 


Substituting Eq. (9) in Eq. (3) and solving for F.., 


Foe = 1.90 Fe9:?9/(b/t)9> Fey 7 (10) 


But, from Eq. (2), Fe, = k’nk(t/b)*, and Eq. (10) be- 
comes 


C’/(b/t)®* (11) 


Fee = 


The coefficient C’ is a function of the effective modulus 
of elasticity, the flange-web ratio, and the mechanical 
properties of the material. It is interesting to note, 
however, that Schuette’s equation for crippling is of the 
same general form as Eq. (4) obtained by the writer. 
In the preceding discussion, certain similarities be- 
tween the empirical crippling equation proposed by the 
writer and those of other investigators have been 
pointed out. In the light of this discussion, it may be 
reasonably concluded that the variation of crippling 
stress with the 6’/t ratio for angle sections possessing 
various degrees of edge support is correctly expressed 
by equations of the form 
Fre = k/(b'/t)” 


(12) 


The exponent » has been assigned a value of 0.75, 
which, while not necessarily the value for best fit of the 
test results, is at least sufficient for the present pur- 
There remain, however, two questions to be 


poses. 
The first question 


answered with regard to Eq. (12). 
is that of the physical significance of the constant k. 
The effective width investigations based on the equa- 
tions of von Karman and Marguerre [Eqs. (5) and 
(7) |] suggest the form of k to be 


m 


k= CV fp" E (13) 
Assuming a value of m of 2 and letting f,. = Fe, 
k=CVF.,E 
and Eq. (12) becomes 
F.. = CV Fe, E/(b'/t)?® (14) 


The writer has studied the results of compression 
stress-strain tests of several investigators, and as a 


Values of C. 


ae Two edges One edge No edge 

V FiyE free free free 
653 ,000 0.307 0.339 0.379 
694,000 0.324 0.352 0.379 
814,000 0.318 0.334 0.350 
Average values 0.316 0.342 0.366 


result of these studies, average values of F,, and E for 
the three materials covered in this paper were found to 
Table IV. Values of C in Eq. (14) are 
also listed in Table IV for the various degrees of edge 
support. It is significant that the variations of C are 
in no case greater than 8 per cent for given edge condi- 
tions. If it is assumed that this variation is attribut- 
able to inaccuracies in arriving at the values of & from 


be as listed in 


test, a single set of dimensionless crippling stress equa- 
tions is obtained by taking the average values of ( 
listed in Table IV. The equations, applicable to all 
three materials, are 


Foe] V Felt = C,/(6'/t)* (15 


where 

C, = 0.316 
0.342 
C. = 0.366 


(two edges free) 
(one edge free) 


Cy 
II 


(no edge free) 


The second question brought out by the foregoing 
discussion arises from the assumption upon which the 
no-edge-free curves are based—-namely, that the per- 
centage increase in crippling strength of the no-edge- 
free angle over the one-edge-free angle is the same as 
the increase of the one-edge-free angle over the angle 
with both unloaded edges free. Evidently, for this 
assumption to be correct, correlation between the no- 
edge-free curve and crippling tests of square and rec- 
tangular tubes should exist. Table V contains the re- 
sults of crippling tests of square tubes obtained from 
references 12 to 15, together with additional tests on 
both square and rectangular tubes of 24S-T aluminum 
alloy conducted by the writer. Fig. 7 is a plot of Eq. 
(15) for the case where C, = 0.366 (no edge free). In- 
cluded on this plot are the test results listed in Table 
V. It will be observed that the agreement of Eq. (15) 
with the test results is very good for b’ ¢ values above 
20. Below this value, Eq. (15) exhibits a considerable 
degree of conservatism with respect to the tests. It 1s 
interesting to note, also, that good agreement exists for 
materials and alloys other than those upon which Eq. 
(15) was developed. 

An additional check (Fig. 8) on Eq. (15) was ob- 
tained by a comparison with the crippling tests on 
FS-lh magnesium-alloy formed zee sections conducted 
by Gallaher.'® Here, C, = 0.342 (angles with one edge 
free), and for the FS-1h material, “ = 6,500 kips pet 
sq.in. and F,, 26.3 kips per sq.in. Again, good 
agreement between test and theory exists. 
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Design Curves 

A design chart based on Eq. (15) is shown in Fig. 9. 
It will be noted upon examination of the logarithmic 
plots of the angle and channel tests (Figs. 2, 3, and 4) 
that few of the sections tested possessed b’/t values 
below 10. Further, the power equation, Eq. (15), 
is of hyperbolic form and is asymptotic with the zero 
h'/t ordinate at an infinite stress. Since the available 
test data furnished no clue as to the nature of the 
true crippling stress curve below 6’/t values of about 
10, the upper limit of the crippling stress was arbi- 
trarily fixed at the yield stress of the material (0.2 per 
cent offset method). The yield stress cutoffs for the 
dimensionless crippling stress curves of Fig. 9 would 
be obtained by substituting F,, for F,, in Eq. (15) and 
solving for the cutoff (b’/t)o. Thus, 


(b’/t)o = C.? (E/Fiy)*” (16) 


and it is seen that the yield stress cutoff varies some- 
what with the edge support condition. 
this variation is small, an average cutoff is used in lieu 
of the three cutoffs dictated by Eq. (16). 

The design chart of Fig. 9 can be placed in a more 
useful form by plotting dimensionless crippling load 


against b’/t. The area of an unequal angle section 


FORMED 


Inasmuch as 
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of flange widths a and 6, and inside bend radius R is 
given by 
A = [(b’/t) — 0.214(R/t) ]2# (17) 


where, as before, 6’/f = (a + b)/2t. 
load of the angle is, from Eq. (15), 


The crippling 


Pee = [Ce V Fe, E/(b'/t)"7] A 





- TABLE V- 


CRIPPLING TESTS OF SQUARE AND 
RECTANGULAR _ TUBES 














+ 








TABLE V 
E F Fey . 

(Kips (Kips (Kips ——— 
Material and 2a 2b t R aa oa om V FE 
Reference (In.) (In.) (In (In.) b’/t Sq.In.) Sq.In.) Sq.In.) 
148-T6!* 6.85 6.85 0.158 0.094 21.6 10,700 31.6 62.2 0.03888 
148-T6!2 6.85 6.85 0.158 0.094 21.6 10,700 31.5 62.2 0.0387 
148-T6!? 4.83 4.83 0.153 0.094 15.8 10,700 43.5 61.4 0.0537 
148-T6"? 1.84 4.84 0.156 0.094 15.5 10,700 43.2 61.4 0.0532 
148-T6"* 4.09 4.09 0.154 0.094 13.3 10,700 53.9 63.2 0.0656 
14S-T6!? 4.09 4.09 0.154 0.094 13.3 10,700 54.7 63.2 0.0665 
148-T6!? 3.58 3.58 0.153 0.094 11.7 10,700 57.6 60.0 0.0720 
148-T6!? 3.58 3.58 0.153 0.094 11.7 10,700 57.9 60.0 0.0723 
148-T6!? 3.34 3.34 0.1538 0.094 10.9 10,700 61.6 63.0 0.0752 
148-T6!2 3.34 3.34 0.158 0.094 10.9 10,700 61.0 63.0 0.0744 
148-T6!? 3.34 3.34 0.158 0.094 10.9 10,700 61.0 63.0 0.0744 
61S-T6!8 2.00 2.00 0.062 0 15.6 10,700 37.8 14.0 0.0552 
61S-T6'3 2.00 2.00 0.062 0 15.6 10,700 36.8 44.0 0.0537 
61S-T6'3 2.50 2.50 0.047 0 26.1 10,700 21.4 14.0 0.0812 
Brass!4 1.00 1.00 0.0065 0 77.0 16,400 11.2 31.4 0.0156 
Brass! 1.00 1.00 0.0056 0 89.2 16 ,000 9.4 28.6 0.0139 
Brass!4 2.00 2.00 0.0065 0 154.0 16,400 6.3 31.4 0.0087 
178-T® +. 00 $.00 0.040 50.0 10,000 13.5 15.0 0.0185 
178-T 3.00 3.00 0.040 37.5 10,000 16.1 16.0 0.0239 
178-T® 2.00 2.00 0.041 25.0 10,000 23.2 16.0 0.0342 
178-T" 1.00 1.00 0.040 12.5 10,000 15.6 16.0 0.0674 
178-T% 0.938 0.938 0.031 15.1 10,000 12.3 44.0 0.0643 
178-T® 0.818 0.813 0.016 25.4 10,000 22.2 14.0 0.0335 
178-T" 0.688 0.688 0.016 21.5 10,000 24.0 44.0 0.0362 
17S8-T 0.563 0.563 0.016 17.6 10,000 32.2 14.0 0.0486 
24S-T* 0.930 1.176 0.070 0.156 7.5 10,500 59.8 15.0 0.0870 
24S-T* 0.934 1.176 0.070 0.156 7.5 10,500 59.3 15.0 0.0863 
248-T" 0.948 1.694 0.052 0.250 12.7 10,500 40.8 15.0 0 0594 
“4S T 0.948 1.696 0.052 0.250 12.7 10,500 10.8 15.0 0.0594 
<48-T* 1.160 1.974 0.033 0.188 23.7 10,500 24.7 15.0 0.0360 
248-T* 1. 160 1.974 0.033 0.188 23.7 10,500 24.8 15.0 0.0361 
“4S T* 1.670 1.670 0.084 0.375 10.0 10,500 55.7 15.0 0.0810 
“4S-T* l 670 1.670 0.084 0.375 10.0 10,500 54.4 15.0 0.0792 
“4S-T* 0.702 0.702 0.050 0.094 7.0 10,500 62.2 15.0 0.0905 
“45-T 0.702 0.702 0.050 0.094 7.0 10,500 62.4 $5.0 0.0908 
“45-T* 1.692 1.692 0.058 0.375 14.6 10,500 46.9 45.0 0.0682 
24S-T* 1.692 1.692 0.058 0.375 14.6 10,500 16.8 45.0 0.0681 


rests conducted by writer. 


Compression yield stress assumed equal to 45,000 Ibs. per sq.in. 








1954 















































































































































226 JOURNAL OF THE AERONAUTICAL SCIENCES—APRIL, 
0.09 T 
0.08 ht CLAD = 
\\ i 
0.07 > - 
}— —}— +4) 245-T3 
0.05 \ 
foc \ \ 
aaa \ 
0.04 . \\ 
A 
NY N FR 
one NS ON GE FREE 
: be TWO EDGES FREE 
~~ 
=e 
0.02 7 
>—_—— 
———— 
0.01 
O 
16) 10 20 30 40 50 60 70 80 
bl a+b 
& ze 


FIG.9 - DIMENSIONLESS CRIPPLING STRESS vs. b/t 


or, from Eq. (17), 


= R/t , . 
= 2C,}1 — 0.214 (~~) | (o’/1)- 
b'/t 


If, the ratio of bend radius to thickness, R /t, is taken as 
three (standard aircraft practice), 
ad 
(1S) 
b’ /t 


bh’ 0,25 
= 2, ( ) [ — 
t 


Eq. (18) tis plotted in Fig. 10 for the three degrees of 
A dimensionless area curve (A //? vs. 


Pree 
?VFVE 


lap 
WV FE 


edge support. 
b’/t) is also included in this plot. 
Crippling of Formed Structural Shapes 


In order to check the accuracy of the proposed 
method of crippling stress prediction, the 24 typical 


structural shapes of Table VI (Fig. 11) were fabri- 
cated from 24S-T3 clad sheet. The dimensions of the 
various shapes of Table VI are mid-line dimensions, 
and the inside bend radii in all cases were '/, in. The 
specimens were fabricated with the grain direction coin- 
ciding with the direction of loading. 

The testing technique employed in this series of 
tests was identical with the angle and channel tests of 
Figs. 2, 3, and 4. The tests were made in a Baldwin- 
Southwark testing machine having a maximum capacity 
of 60,000 Ibs. and intermediate load ranges of 24,000 
and 6,000 Ibs. 

The analysis technique employed by the writer 1s 
shown in Fig. 12. Dimensionless crippling loads were 
obtained from Fig. 10 and were converted to dimen- 
sional quantities using the average values of F,, and 
E of Table IV. 
t?5-(A /t?) or (0.0255)2(180.4) = 


| 


The calculated area of the section 1s 
0.1173  sq.in. The | 
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TABLE VI 
Results of Crippling Tests of 24S-T3 Clad Formed Structural Shapes 
Prediction 
Shape a b ( d e t A Test F, ’ 
Al 1.000 0.434 0.0250 0.0672 26.2 25.0 
A2 1.568 0.437 0.0315 0.1188 22.9 23.0 
Bl 1.205 1.340 1.236 0.447 0.0255 0.1833 20.8 22.3 
B2 1.198 1.610 1.001 0.445 0.0320 0.2040 27.1 27.0 
ei 2.457 0.980 0.449 0.0200 0.0989 7.3 16.6 
C2 2.487 0.987 0.414 0.0315 0.1571 24.1 23.6 
C3 2.412 0.892 0.408 0.0382 0.1892 27.5 28.3 
C4 2.236 1.262 0.428 0.0505 0.2660 31.4 32.1 
C5 1.813 0.862 0.425 0.0258 0.1045 24.0 24.4 
D1 1.23 2.260 0.853 0.448 0.443 0.0205 0.0984 17.4 17.1 
D2 1.255 2.195 0.875 0.429 0.429 0.03882 0.1870 28.2 27.8 
D3 0.847 1.887 0.847 0.4381 0.4381 0.0258 0.1070 23.3 24.2 
D4 1.244 1.950 1.244 0.423 0.425 0.0252 0.1261 19.0 20.0 
D5 1.250 1.904 1.250 0.415 0.415 0.0383 0.1886 Ys ee 26 . 4 
El 1.991 0.653 0.961 0.0252 0.0868 15.8 15.1 
E2 1.082 0.492 0.844 0.0253 0.0558 22.6 20.9 
E3 2.141 0.418 1.028 0.0504 0.1715 24.7 25.4 
E4 1.213 0.408 0.9385 . 0.0255 0.0620 21.0 19.7 
Fl 2.331 0.856 1.241 0.424 0.424 0.0318 0.1571 23.9 23.0 
F2 2.118 0.976 1.266 0.732 0.442 0.0318 0. 1651 23.8 23.1 
F3 1.788 0.448 0.953 0 0.438 0.0500 0.1680 35.0 34.8 
Gl 0.977 1.561 0.881 0.422 0.0318 0.1127 6.2 23.9 
G2 0.936 2.457 0.8638 0.425 0.0255 0.1141 17.8 17.4 
H 3.668 1.528 0.393 0.943 0.0318 0.1933 25.0 24.9 
: ee B63 ~ 
+b I } ' 
- HAPE - ' is 
SHAPE-A | Ss £6 MAT'L. 24S-T3 CLAD 425 (Oo ‘J ] 
; ‘ TEST Re 2035 LB. 5 
a t t AREA 0.1141 SQ.IN. 2.457 
! b t 
I fee 17,700 PSI. .0255 > 
a Q —— 
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~—~Cr 1 
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measured area of the shape (obtained by weighing the 
specimen and using a material density of 0.100 Ib. per 
cu.in.) was 0.1141. The disagreement of area values is 
due to the fact that the design chart was developed for 
sections having an R/¢ ratio of 3, whereas the test speci- 
men R/t was 4.90. A comparison of test and pre- 
dicted crippling tests is shown in the last two columns 
of Table VI. The ratio of predicted-to-test crippling 
stress for the 24 shapes tested shows a scatter between 
0.912 and 1.121, while the average ratio is 0.993. 


Predicting Crippling Stresses at Elevated Temperatures 


Eq. (15) is in a form that is readily adaptable to pre- 
dicting crippling stresses at elevated temperatures. All 
that is required is to have on hand curves giving the 
variation of F,, and / with temperature.* While the 
physical significance of C, in Eq. (15) is not known, an 
analysis of Schuette’s equation [Eqs. (3) and (11)] 
indicates C, to be a function of the effective modulus 
of elasticity and the compression yield stress, both of 
which vary with temperature. The present investi 
gation demonstrates, however, that at room temper- 
atures C, is a constant for a given edge condition. Fur 
thermore, an analysis by the author of the elevated tem- 
perature crippling tests of extruded H-sections con- 
ducted by Heimerl'® indicates that Eq. (15) will yield 
reasonably accurate values of the crippling stress of 
Unfortunately, no test data on ele- 
formed sheet-metal 


these sections. 
vated temperature crippling of 
sections are known to the author; hence, the absolute 
certainty of the accuracy of Eq. (15) at elevated tem- 


peratures cannot yet be established. 


CONCLUSIONS 


From the results of this investigation the following 
conclusions on the proposed method of crippling stress 
prediction have been drawn: 

(1) The crippling stress of equal and unequal angles, 
channels, zees, and rectangular tubes can be predicted 


directly from the equation 


Fi. = CV FE/('/tf* (15) 
where 
‘y = compression yield stress, lbs. per sq.in. 
F yield st It rsq 
Ek = Young's modulus in compression, Ibs. per 
sq.in. 
b’/t = equivalent b/t of section = (a + b)/2t 
C. = coefficient that depends on the degree of 


edge support along the edges of contiguous 
angle elements. Specifically, these values 
are: 0.316 for angles with two edges free, 
0.342 for angles with one edge free, and 
0.366 for angles with no edge free 


*This assumes that time effects are negligible—i.e., that 
creep buckling need not be considered. 
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(2) The crippling stress of structural shapes other 
than those mentioned above can be determined by 
dividing the shape into a series of angle sections and 
computing the crippling loads of the individual angles 
from Eq. (15) and the areas. The weighted crippling 
stress is obtained from 


F > (crippling loads of angles) 


>> (area of angles) 


(3) Although experimental evidence is lacking, the 
proposed crippling stress equation is believed to apply 
to the crippling of formed structural shapes at elevated 
temperatures, provided the variation of compression 
yield stress and Young’s modulus with temperature is 


known. 
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Linearized Theory of the Oscillating Airfoil 
in Compressible Subsonic Flow 


R. TIMMAN? 
Technological Unwersity, Del/t, Holland 


SUMMARY 


Extending previous work, a series expression in Mathieu func 
tions is given for the pressure distribution on an oscillating airfoil 
in two-dimensional subsonic flow. Using an identity between 
different solutions of the corresponding boundary value problem, 
it is shown that the resulting formula satisfies the reciprocity re 


lations in nonstationary flow. 


NOTATION 


x = coordinate in direction of main flow, fixed to the 
fluid 

x’ = coordinate in direction of main flow, fixed to the 
airfoil 

t = time 

: om fay ® = coordinates obtained from x, y, ¢ by a Lorentz 
transform 

V = flight velocity 

( = sound velocity 

M = V/c, Mach Number 

B = (1 — M?*) ‘2, Lorentz (Prandtl) factor 

oy = velocity potential 

rf) = acceleration potential 

v/2r = frequency of airfoil oscillations 

2! = airfoil chord 

w = vl/V, nondimensional frequency parameter 

Q = B*w, compressible frequency parameter 

G(X, } Green’s function of the second kind 

X,, V1) 

sen(n) odd Mathieu function of order 

Ne,,°?(€) associated Mathieu function of the third kind 
of order n 

w(x’) = given normal velocity on the airfoil contour 

Pp = pressure 

p = density 

Il = pressure difference on the airfoil 


(1) INTRODUCTION 


T HIS THESIS,'! the author gave a solution to the 
problem of the determination of the pressure dis- 
tribution on an oscillating airfoil in compressible sub- 
sonic flow. This solution was presented afterwards in 
corrected form by the author and van de Vooren,’ to- 
gether with an extension to the case of an oscillating 
airfoil with flap. The resulting formulas were evalu- 
ated numerically by the Computation Department of 
The. Mathematical Centre at Amsterdam under the 
direction of Van Wijngaarden,* and the results were 


Received August 17, 1953. 

* This research was sponsored by the N.I.V. (Netherlands 
Aircraft Development Board). 

+ Professor of Mathematics. Also, Consultant for Theoretical 
Acrodynamics, National Luchtvaart Laboratorium, Amsterdam 


published in this JOURNAL.‘ The numerical values 
were not in good agreement with the values obtained 
by Dietze by a numerical method for the solution of 
Possio’s integral equation.® A short time after pub- 
lication in this JOURNAL, Fettis® raised severe objections 
to the values given in reference 4, because they did not 
satisfy the reciprocity relations for nonstationary flow, 
relations that at the time of publication were unknown 
to the authors. 

Two possibilities were open: either an analytical error 
in the final formulas given in reference 4, or an error 
in the numerical evaluation by The Mathematical 
Centre. 

In this paper it will be shown that the formula for the 
pressure distribution, as given in reference 4, satisfies 
identically the reciprocity relations in the form given 
by Flax.’ This result was obtained by a re-examina- 
tion of the theory inspired by the work of Hofsommer* 
in incompressible flow, which led to the establishment 
of a new identity which enabled the author to show 
that his pressure distribution satisfied the reciprocity 
relation. Moreover, using this identity, it was possible 
to simplify the expression for the pressure distribution 
considerably, thus getting a complete counterpart to 
the well-known results in incompressible flow. The 
author was informed by private communication that 
Kassner had already obtained similar results. 

During the preparation of this paper the numerical 
error has also been detected, and the author has been 
informed, that the corrected numerical values in a par- 
ticular case do check the reciprocity relations up to four 
decimal places. The author is indebted to Drs. van de 
Vooren and Van Spiegel of the N.L.L. for many inter- 
esting discussions on this subject. 


(2) Bastc EQUATIONS 


In a compressible medium the velocity potential 
(x, y, t) and the acceleration potential g(x, y, ¢) of 
an unsteady two-dimensional disturbance satisfy the 
linearized wave equation 


G'rr + Py (1/c*) gn = 0 (1) 


where x, y are cartesian coordinates fixed with respect 
to the fluid. 
The relation between the two potentials is given by 


e(x, y, t) = B(x, y, t) (2) 
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OSCILLATING AIRFOIL 


We now consider a thin airfoil moving with subsonic 
velocity V, Mach Number 1, in the direction of the 
negative x-axis. The equation for the potentials ¢° and 
@ with respect to a set of coordinates fixed to the air- 
foil will be obtained from Eq. (1) by the transformation 
(Galilei Transformation) 


x’ =xt Vt y=y V=t (3) 
It is, however, well known that the wave equation, 
Eq. (1), is invariant for Lorentz transformations of the 
form 
X=—Px+ VT) Yu=y, T= 
Blt + (iW c)x] (4) 
where 
(5) 


The relation between the Galilei coordinates, Eq. (3), 
and the Lorentz coordinates, Eq. (4), then reads 


X = px’, Y=y’, T = (t'/B) + (8Mx'/c) (6) 


We now consider harmonic oscillations of the airfoil 


with frequency v, and hence we may introduce Lorentz 
disturbance potentials (XY, ¥) and g(X, V’) by 


4 (x, y, t) = ®(X, Ve?" = &(X, Veo" *"-e" 


(7) 
g(x, y, t) = ¢ (X, Ve? = o(X, V)et?o"*- 0 
(S) 
Then (XY, ¥) and ¢(X, VY) satisfy the equation 
¢xx + gry + k’g = 0 (9) 


where 
k = Bv/c 
The relation (2) between the potentials transforms into 
o(X, Y) = ivB*@(X, Y) + BVOx(X, Y) (10) 


We assume the airfoil to extend over the segment 
-l <x’ < +1 of the x’-axis, or in the Lorentz co- 
ordinates over the segment 


—-sl<X<+t+8l 


The boundary conditions for the potentials are given 
by the oscillatory motion, which is assumed to be 
known. For an infinitely thin airfoil, the normal ve- 
locity is equal at opposite points on the upper and lower 
side, and we get the boundary condition for the velocity 
potential 


$,.(x’, + 0) = #,%x’, —0) = 
w(x’)-e” =-—(m ss 42 (11) 


Where w(x’) is the amplitude of the given velocity. 
This imposes on our problem the symmetry condition, 


Hx’, y, t) = —O%(x’, —y, 2) 
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Similarly, we get a boundary condition for g(x’, y, ¢) 
by considering the normal values of the acceleration 


0,w(x’)e™ = O,w(x + Vt)” = (ivw + Vw,-)e” 


which gives 


, 


¢,(x’, + 0, t) = ¢,(x’, —0, 2) = 
(1yw + Vw,-)e = (12) 


These boundary conditions are readily transformed into 
boundary conditions for the Lorentz potentials, 


, e —ivB MX /< ai 


&;(X, 0) = w(x W(X) (13) 


—ivBMX/c_ 


ivB?W(X) + BVWy (14) 


gy (X, 0) = (tvpw + Vw,)e 


By giving the boundary conditions on a closed con- 
tour, solutions of the wave equation are not uniquely 
determined in the outer region. Here Sommerfeld’s 
condition must be added, which states that, at infinity, 
the solution must represent traveling waves, diverging 
from the airfoil. For a closed regular contour, the 
solution is uniquely determined by these conditions. 
The behavior of the solutions near the ends of the seg- 


ment will be investigated separately. 


(3) THE SOLUTION OF THE BOUNDARY VALUE PROBLEM 


Considering the two sides of the airfoil together as 
forming a closed contour, the problem is a boundary 
value problem of the second kind for the wave equation, 
and its solution depends on Green's function of the 
second kind. 

This function, denoted by GLY, V; X1, ¥1) is defined 
as a solution of the equation 


Ag + k?¢ = 0 (1) 
satisfying the following conditions: 
(i) At -—psl < X < + Bl, 
Gy(X, + 0; X1, ¥i) = G(X, — 0; X%1, V1) = 0 


(ii) At Y = X,, Y = Jj, the function must have a 
logarithmic singularity. 

(iii) At infinity, the function must satisfy the Som- 
merfeld condition for outgoing waves. 

In physical language the function represents the po- 
tential of a source in the point X,, ), with the segment 
—Bl < X < Blas a fixed wall. 

Unfortunately it has never been possible to obtain 
an explicit expression for this function. For this reason 
we follow another procedure for the solution of the 
boundary value problem, which then will give us an ex- 
pression for the function G. 

We introduce elliptic coordinates 


X = Bl cosh é cos 7, 
Y Bl sinh & sin 7, 


—Eé>0 
—r<nster 


II 


The derivatives then transform by 
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31d. _ sinh — cos nO: — cosh é sin 70, For large values of 1 it can be shown (see reference 2, mea 

cosh? £ — cos? 7 Appendix 2) that the terms in this series are approxi- 

, mated by the corresponding terms of the ordinary 
cosh & sin 70; + sinh & cos 70, F ier expansion in the inc sssible case ; ; . 
Bldy = — . ourier expansion in the incompressible case and hence whit 

cosh* £ — cos* 7 the convergence properties in the two cases are the 

and the wave equation then passes into ee 

This latter expression is App 
Gee + Gm + Wh*%(cosh 2 — cos 2n)e = 0 (2) | ; poi 
s ~ @ roe 2 ‘ ° _ 
where CX.S BW) =  & sin ”n- / JF(m) sin m mm, dn, I 
1 on rT /0 
pon 
h= BRI 2= Bul ¢ (9 the | 


Particular solutions of this equation having the proper and it is well known that in this case summation and the : 


symmetry can then be found by separation of variables, integration may be interchanged. This can also be tion 


which gives solutions of the type’ done in our case, and we write down 


Ne,‘ (&)s€n(n) " 
g(é, 7) = / f(m)+sin m°G(E, 7; m) dm 9g with 
0 


where se,(n) is the odd Mathieu function of order Eq 
and Ne,‘ (&) is the corresponding associated Mathieu , 
s ; ‘ : Re where 
function such that the particular solution satisfies 
Sommerfeld’s condition. 2 2 Ne, ((z) 
. '/ te ad 4 n Ss 
We now assume a series solution of the form G(s, 7; m) = > Ve, (2)"(Q) Sén(m)sen(m) (10 
T n 1 - on Z 
a, Ne, (E)se, (7) a . ; . ; In o 
p> r . [his series vanishes at the end points of the interval equi 
; es: viz., até = 0,7 = Oandé = 0,7 = =. itd 
In order to determine the coefficients a, we write down : ; , van 
A i pile ; ; In the incompressible case, the series can _ be 
the boundary condition in elliptic coordinates at & = 0 
5 summed 
g: = sin n-f(n) (3) 
i ») e ns 
If f(n) is assumed to be a continuously differentiable, G'(é, 7; m) = > - Sin 2n-Sill nN; Hen 
‘ . ‘ : : YT 21 # 
even function of n, f(n)-sin » can be expanded in a ; 
uniformly convergent series of odd Mathieu functions I cosh £ — cos (n + m) 
: n 
24 cosh € — cos (n m1) 
sin n:f(n) = 2 On Sey (7) (4) 
For the velocity potential the function f(y) 1s 
where equal to @/IV(m) and the corresponding — solution 
9 fF . : : = The 
2, = f(m) sin mse,(m) dm (9) > iti 
T Jo : ; ‘ ° 
&(~, n) = Bl / Wm) sin m-G(é, nim) dn (12 
Differentiating our series solution and putting & = 0, 1d) 
we obtain ; ; ; 
and for the acceleration potential 
o(0, 7) >> a,: Ne, ‘?’(0)+se,(n) (6) ; em is , 
" , bao “ f(n) = Bllive?W — (V/l sin n) W,] 
and equating corresponding coefficients in Eqs. (4) and — and, hence, thec 
(6), we get the solution the | 
" -FO Ne, (é) g(t, n) = BI | (QW sin m — IW, )G(é, 1; m) dm bee 
ey = 24." — "Se, () . 
n=l Ne, ‘’(0) 13) U 
— Ne,‘ (€) $f pres 
x — ray P *S@,(N)* t(m) Sil 1° S@,( Mm) Gm where Kut 
n=1 Ne, ‘’(0) T JS 
(7) Q = Byl/V Act 
shor 
A 
the 
cele: 


(4) THE FUNDAMENTAL IDENTITY 


The two solutions to the same problem obtained in the preceding section are not identical. This is apparent 
we calculate the acceleration potential g*(X, VY) corresponding to the velocity potential, Section (3), Eq. (12) by 
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means of 


e*(X, Y) = (V/1) (Qe + Ble, (1) 

which, on the airfoil, where £ = 0 becomes 
¢*(0, n) = (V1) [s20(0, n) — (1/sin n)-®,(0, n) } (2) 
Apparently this expression is infinite at the end points » = 0 and n = z, but Section (3), Eq. (13) vanishes at these 


points. 

The difference potential must have a zero normal derivative on the airfoil, it must be regular in all internal 
points and singular at the end points, and the kind of the singularity must be the same as that of ¢*(£, 7). Now 
the function G(&, 7; m) given by Section (3), Eq. (10) has these properties, except for the singular behavior. As 


the singularity corresponds to a differentiation, we differentiate with respect to 7. Then G, (&, 7; ) is still a solu- 


tion of the wave equation in £ and yn, and we may put either 7; = O ory, = 7. We now show that 
¢(é, 9) = o*(é, 7) + BVAAG, (&, n; Q) + BV XAG, (é, 9, F (3) 
with suitable coefficients \; and A». These coefficients are determined by the condition, that the right-hand side of 


Eq. (3) must remain finite at 7 = Oand yn = 7, or that 


lim [—(V/1)®,(0, n) + BVA; sin n-G, (0, 7; 0) + BVA» sin n-G, (0, n; 2) ] 0 (4) 
n—O0 
lim [—(V/1)®,(0, n) + BVA, sin n-G, (0, 7; 0) + BVA2 sin 7-G, (0, 7; 7)] = O (5) 


In order to evaluate these limits we remark that in the neighborhood of a singular point the solutions of the wave 
equation can be approximated by the corresponding solutions of the potential equation and hence, up to a function, 
vanishing in the limit 

G, (0, 7; 0) ~ G, (0, 7; 0) sin n/m (1 — cos n) 


G, (0, n; 7) ~ G,'*(0, 7; x) —sin n/a (l + cos7 


Hence we find, calculating the limits in Eqs. (4) and (5), 


7 . m ‘ 
Ay / 11’(m) sin m°G,(0, 0; m) dm (6) 


T r . . —_ 
2 a / W(m) sin m°G,(0O, 7; m) dn (4) 


—e 


The remaining difference function now is finite at the end points; its derivative vanishes on the airfoil and, hence, 


itis identically zero. Thus we have found the important identity, 


Bl Ox / Ii(m) sin m-G(E, 1; m) dm + / W,, (m)-G(é, 03 m) dm + A1-G,,(&, 0; O) 4+ ALG, (E, 3 7) O (S) 


(5) Tue Kutta ConpDITION 


Neither of the two solutions given in the preceding section gives the solution to the actual problem of airfoil 
theory. The regular solution g(Y, JV) for the acceleration potential was obtained from the boundary values of 
the normal acceleration, which were found by differentiating the given values of the normal velocity, and it cannot 
be expected that integration again yields back the proper velocity. 

On the other hand, the velocity potential 6(Y, Y) adapted to the given values of the normal velocity gave a 
pressure distribution (acceleration potential) infinite both at the leading and trailing edge, which is contrary to the 
Kutta condition, applying in the actual problem, which requires that the pressure vanishes at the trailing edge. 
Actually the Kutta condition only requires that the pressure remains finite at the trailing edge, but it can be 
shown that in this case the acceleration potential must vanish. ) 

As we already remarked, ¢(, 7) vanishes both at the leading and trailing edge and also G, (&, n; 7) vanishes at 
the trailing edge — = 0, 7 0 and does not contribute to the normal acceleration. Hence, we assume for the ac- 
celeration potential, 


7 


(én) = BV | / }i2W sin m — W,$ -G(E, 0; m) dm + aoG,,(& 7; | (1) 
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where ad) must be determined by the condition, that the normal velocity, obtained from Eq. (1), must coincide 
with the given normal velocity. As the normal accelerations do already coincide, it is sufficient to consider a single 
point of the airfoil and, following our previous choice, we take the leading edge. 

The derivatives , and g, being related by the same formula as ® and g—viz., 


gy(X, Y) = V/1l (iQ, + Bl’y,) (2) 


#, can be expressed into gy by integrating this first-order differential equation, 


CA. ¥ : J. ~ — x) CGR Ae @& 
®,. r, = "Py. - ) y 3) 
} ) V ane exp | al | a (9) 


if we assume that 4; = Oat XY = —o. 
Substituting, we get the following equation for the determination of dp: 


— pl - ae 
W(x) = / eittl +(x miax’-| ay / GOQW sin m — W,)-G(e’, n’sm) dn | + 
1 J0 : } 0 
e 


X’ = —£l cosh ¢’, Oy = (—1/Bl sinh é’)O, 


we get 


Wr) = —e® / e sea | (QW sin m — W,) G, (&’, 7; m) dm — 
JI 0 


Q 2 cosh &" 7 es , es e 
ae’ / * ie cas Gyy, (E', ©; mF) dé (5 
0 
which is the same as Eq. (7.25) of our previous paper.” 


By means of the fundamental identity, the equation can be considerably simplified. Substituting 


[ W, -G(é’, n’; m) dm = —BlOdxy: [ W sin m-G(é’, 0’; m) dm — MG, (&’, 9’; 0) — WG, (E05 
0 J0 


into Eq. (4), we get 


*—pl = os 
W(r) = ate f oo * as. [2:2% J W sin m-G(é’, ’; m) in| -~ 
x 0 Y=0 


at , fo 
nett f * wade [OvG,, (’, 0’; 0) Jy=0 dX’ + (Az + ao) -e™ / ef X'/PN0,.G, (é’, 0'3 w)]y-0 dX’ + 


—l “x 
ine* f fr jar f W sin m-G(é’, n’; »)| dX’ (6) 
x 0 Yy=0 
and by partial integration of the first term, 


sa 7~— Bl ° x 
W(m) = ploy f W sin m-G(0, +; m) dm —iQe™ / gee ax’ la f W sin m-G(é’, 7’; m) dm|y=o0 + 
0 x 0 


— pl Ms *—£,l i 
die” [ ef X/BUO,.G,(é’, 0’; 0) ]y-0 dX’ + (Ae + ane / ” aaa [vG, (&', 0; m)]y-0 dX’ (7 
Regarding Section (3), Eq. (12), we see that the first term of the right-hand side is equal to W(2). Introducing 


the integrals 


~ 


— pl x 
Dil ate-i J ei®X"/8l VG, (8, 0/3 O)Iy-0 dX" = f ens CO Oe 


7 — pl x 
xX’ / ] r —ift cosh £’ ’ 
R, = — / elt XP! TOG, (€’, 0’; wy -0d X’ = / et Gag lks 0 wee 
— ' 0 


a) = —dz — (AiRi/ Re) (9) 


We find 
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T(Q) = 
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R,/R2 


we finally obtain, for the Lorentz acceleration potential, the expression 


¢(é, 1) = av| f (QW sin m — W,)-G(E, 0; m) dm — (x + MT)-G, (&, 0; | (10) 


It should be remarked here that R, and R, cannot be calculated directly by termwise integrating the series ex- 


pression for Gyn, but must be calculated as described for R, in,* by splitting off the “incompressible” singular part 


of G,,, and integrating termwise the remaining correcting series that is uniformly convergent. The value of the 


contribution of the singular part can be taken from the known result in the incompressible case. 


From the 


Lorentz acceleration potential, we now calculate the physical acceleration potential, related to the pressure differ- 


ence by the formula 


IH(x’)-e' = pt(x’, 0, t) — p-(x’, 0, 1) = 2pe (x’, 0, 2) (11) 


This gives 


Mx’) = 2p-e?!*"/<. (0, ») = 2pBV-e@M? 59 | / @2W sin m — W,)G(0, 0; m) dm — (Ae + AT)G, (0, 0; | 


i./0 


= 2pBV-e'@ a f e MP cosa (iw sin m — Ww, )G(O, 7; m) dm — (2 + AT)G, (0, 7; | (12) 
0 : 


where w = vl/V. 


(6) THE REcrpROcITY RELATION 
We now prove the reciprocity relation in the form 
given by Flax’ and show that the expression in Section 
(5), Eq. (10) satisfies this identity. 
Consider two solutions of the wave equation g and 
$, both antisymmetric in Y and vanishing in infinity. 
Green's theorem enables us to write 


. 


| g(X’, 0) -By(X’, 0) dX’ = 


J | &(X’, O)-gy(X’, 0) dX’ (1) 


Considering now ¢ as the Lorentz acceleration potential 
corresponding to an airfoil moving in the direction of 
the negative x-axis, the corresponding Lorentz velocity 
potential satisfies 


e(X, Y) = 1vB?O(X, Y) + BV Ox(X, VY) 


The second member of Eq. (1) now can be written as 


inst f by dX + av | dy» dX’ 


Fy partial integration of the second term we get 


av [aay | oo av f by Py dX’ 


If we consider @ as the Lorentz velocity potential of an 
airfoil moving with velocity V in the direction of the 
positive Y-axis, the corresponding acceleration poten- 
tial is 

e = 1vB’h — BVO, (2) 


Moreover, at — ©, d, = Oand at +, # = 0, and the 


integrated parts vanish. This gives the result 


/ gb, dX’ = / o- Py dX’ (:3) 


Considering further that in the regions of the x-axis 
before and behind the airfoil both ¢ and ¢ vanish, we 


may write 


> +21 > +4 gl 
/ gb, dX’ = / oe: Py dX’ (4) 
pl Bl 


Returning to the actual pressure distributions and 
vertical velocities, we know that 


I(x’) = 99-6 ™ “.9(X) 

i(x’) = eer *.o(X) (5) 
w(x’) = ef" /<.b(X) 

w(x’) = eW OP Mx BX) 


and we get Flax’s result. 


II(x’)- w(x’) dx’ = | I(x’) w(x’) dx’ (6) 
‘i 


We now show that ¢(X, 0) as given by Section (5), Eq. 
(12) satisfies the identity, Eq. (4). 

Consider a reverse flow with vertical velocity 
Then the normal acceleration is 


o tul 


w(x’ )e 
on the airfoil. 
(iva — V@,-)e'™ 


and the boundary condition for the Lorentz potential 
is 
oy(X, 0) = ivB2?W(X) — BV Wy 
where 
W(X) = w(x’) eM" 
Accordingly the regular part of the acceleration po- 
tential is 
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e(§, 1) = av: f (iQ sin m + W,,)-G(E, 2; m) dm shee 
0 ‘. 
(7) 1 = = W(m) sin mG, (0, 0; m) dn 
For the singular solution we now must take - - ae . (8) 
GG, (E, 7; 0) ~~ = f W(m) sin m-G,(0, viedael 


ae : , ’ we get 
where @ is determined by integration along the seg- 


ment —o < X¥ < —, & = —rA; — dbo TD 


Now substituting the complete expressions for ¢ and 
¢, where 


¢(0, ») = BV / (QW sin ny + W,.)-G(O, nim) dm — (Ar + A2T)G,, (§, n; O v 
into Eq. (4), we get 


/ / (QW sin m — HW, G(0, n;m) dm: W(n) sin ndn — | / (QW sin m + W,,) x 
J0 Jd 0 Jo 


G(O, n; m) dm: W(n)sin n dn — (Ao + AT) / G, (0, n, wr): W(n) sin 7 dn —(Xy + Yel): / G,, (0, n; 0) X 


W(n) sin n dn = — / / W,,-G(0, n; m)W(n) sin » dn dn, — / / W,,-G(0, n; m)W(n) sin n dn dy, + 
2 = es 
(Ae + Ail) deo = (A, 2 i del) +dy 
Tv Tv 


Substituting the identity in the form 


T | "fr 
[ W,,-G(O, 23m) dm = ae / W(m) sin m-G,(0, 1; m) dm — MG,.(0, 730) — AWG, (0, n; x) (10) 
y S sin 0 


this passes into 


-[f W(»)-W(m) sin m°G,(0, 7; m) dm dn+n f W(n)-G,,(0, 0; 0) sin » dn — 
0 0 0 
T "x = < » 7 2 
f | IV’, -G(0, n;m)W(n) sin n dn dny + Xo / W(n)-G,, (0, n; 3) sin n dn + ~ (AcdA2 — Ari) 
0 0 0 . us 


and by partial integration of the first term with respect to 7, we get 


W(0)- Z W(m)-sin m:G(0, 0; m) dm — W(2)- f W(m) sin mG(O, x, m) dm + 
0 0 
| f W,-W(m) sin m:G(0O, 7; m) dn dn — F / W,,-W(n)-G(0, n; m) sin 7 dn dm + 
0 0 0 J 0 
9 


9 


~ (Aid — Aoke) + = (Ave — MA1) = W(O) peg. (O) — Wa) Bey (7) = 0 (11) 


Tv Tv 


since the regular part of &(é, 7) vanishes at the end points of the interval. 
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Calculation of Compressible Cascade Flow 


by the Method of Flux Analysis’ 


SHIGEO UCHIDAt 
Institute of Sctence and Technology, University of Tokyo 


ABSTRACT 


A unified graphical method for analyzing the problems of com- 
pressible flow through a cascade of arbitrary airfoils is presented. 
The method is based on the conception that the partial dif- 
ferential equation expressed in terms of suitably chosen curvi- 
linear coordinates may be replaced by an ordinary differential 
equation in a single variable, the remaining variables being re- 
garded as parameters. The computation reaches, in many cases, 
a convergent solution after two or three iterations of 2 graphical 
procedure. It will be particularly useful for cascades of high 
solidity. 

Comparisons with the exact solutions of incompressible flow 
past cascades of modified Joukowski airfoils and with the ob- 
served results in a cascade wind tunnel are satisfactory consider- 
ing the accuracies of the graphical methods. Applications are 
also made to the study of turbine blade sections, in which effects 
of blade profile, solidity, and Mach Number are investigated. 


NOMENCLATURE 


x,y = cartesian coordinates 

a, B = orthogonal curvilinear coordinates 

5, 0 = arc length of a, 8 axes, respectively 

ha, hg = parameters in curvilinear coordinates 

V = velocity vector 

u, = a, 8 components of V, respectively 

p = density 

p = pressure 

( = sound velocity 

Y = ratio of specific heats 

K = constant 

, ¥ = velocity potential and stream function, respectively 

F(a) = integration constant 

N = value of m at 8B = Bp 

l = chord length 

t = pitch distance along cascade surface 

pi, p = densities of undisturbed flow, up- and down- 
stream, respectively 

pi, pe = pressures of undisturbed flow, up- and down- 
stream, respectively 

Wi, @. = total velocities of undisturbed flow, up- and down- 


stream, respectively 

Wa, W@a2 = normal components of w; and we to the cascade 
surface, respectively 

= tangential components of w, and we to the cascade 


Uy, U = 
surface, respectively 
8, Be = angles of attack of w; and we, respectively 
r = circulation around the unit height of an airfoil 
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i = lift force acting on the unit height of an airfoil 

6/ = summed-up partial torque acting on the blade 
eiement defined by 7 and 7 + 6r 

6Cy = partial torque coefficient 

= radius of the circumferentially arranged cascade 

yy = number of blades in the circumferentially ar- 
ranged cascade 

f,n = cartesian coordinates parallel and perpendicular 
to the chord line, respectively 

a,b = half length of major and minor axes of ellipse 

€ = b/a = ratio of radii of ellipse 

n = order of general ellipse 

f, k = percentage maximum camber and its distance from 
leading edge, respectively 

e, m = percentage maximum thickness and its distance 
from leading edge, respectively 

o = angle of incidence of the airfoils 

M, = w,/c, = Mach Number of the undisturbed inlet 
flow 

h,g,c,d = parametric constants 

Subscripts 

0 = quantities in the state at rest 

a, b = boundary values 

m = orders of the approximate solutions 


(1) INTRODUCTION 


fs OF THE FUNDAMENTAL AERODYNAMIC PROBLEMS 
in turbomachinery is the flow of a compressible 
fluid past cascades. It will be expected, by predicting 
the flow characteristics through such cascades of air- 
foils, to obtain appropriate blade designs for a turbine 
or an axial-compressor of high efficiency. To date a 
number of methods have been proposed to solve the 
idealized flow about airfoils arranged in a cascade, in- 
volving analytical or semianalytical methods using 
conformal transformation, superposition of interfer- 
ence, and series expansion;'~® graphical or semi- 
numerical procedures;’~'* and mechanical or electrical 
techniques.'*—"4 

In the present paper a unified graphical method fos 
analyzing the problems of compressible flow past cas- 
cades of arbitrary airfoils is presented. As the principle 
of the method, which has been briefly stated in a 
previous paper,'! is based on some general ideas, 
similar methods may be applied to differential equations 
of higher order, for example, to those of problems in 
elasticity or plasticity and others. 

The partial differential equation expressed in terms 
of orthogonal curvilinear coordinates may be replaced 
by an ordinary differential equation by assuming one 
of the coordinate axes to be nearly coincident with the 
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curves of solution. In the following method of ‘‘flux 
analysis,’ a stream line is chosen as one of the coordi- 
nate axes. Then approximate solutions are obtained by 
simple integrations, while the stream lines first assumed 
can be adjusted until they agree with the final solution 
resulting from the iterative processes. Close approxi- 
mations may be obtained by the correction of the turn- 
ing angle or the direction of outlet flow by the value of 
circulation calculated from the velocity distribution 
Taking special care with the 
calculations in the neighborhood of the leading and 


along the airfoil surface. 


trailing edges, fairly good results are obtained after 
two or three iterations. Comparisons both with the 
accurate cascade flow of modified Joukowski airfoils 
and with experiments on a cascade of turbine airfoils 
are shown to be in good agreement. Applications are 
also made to the basic investigations of turbine blades 
of high efficiencies. Introducing a series of turbine 
cascades having a mean camber line in the shape of a 
general ellipse, effects of the maximum camber and its 
location, the maximum thickness, the order of ellipse, 
the pitch chord ratio, and the Mach Number on the 
surface distribution of velocity and on the magnitude of 
the circulation are studied and presented, with several 
examples of flow patterns. 


(2) FUNDAMENTAL EQUATIONS REFERRED TO THE 
ORTHOGONAL CURVILINEAR COORDINATES 


The steady two-dimensional isentropic flow of a 
nonviscous compressible fluid is governed by the 
equations of continuity, irrotationality, and the isen- 
tropic relation. Referring to the orthogonal curvilinear 
coordinates, the fundamental equations are expressed 
as follows:" 


— l O(/ ) O(ha v) 
div (pf) = | (Nig pu 4 (Hap |- 0 (1) 


hahglL Oa OB 





a 





Fic. 1. Curvilinear coordinates 
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re l O(hgv) O(hatt) 
rot V => | al i. | = () (2) 
NahglL Oa OB 
p= Kp (3) 


where /ia and /ig are the parameters of orthogonal curyj- 
linear coordinates corresponding to the extension fac- 
tors of the networks. 


ha = lim . = | s) + (22 
5a—0 O@ \ el Oa 


hg = lim 7 = (=) + (2) 
6pB-—-0 6B \ og og 


As a function that identically satisfies the condition 
of continuity, a stream function can be introduced, de- 
fined as follows: 


11 op | 1oy 


— 71 \.) 


phg Op’  phada 


Substituting Eq. (5) into Eq. (2), the differential equa- 
tion of y referred to the orthogonal curvilinear co- 
ordinates is obtained. 


(=) a) 
0 a\phada = Brae 


The variation of density can be most simply expressed 


= 0) (6 


by quantities at the stagnation point, which is given 
by the use of Bernoulli's equation combined with the 
isentropic relationship as follows: 


. _ vw 3 =)" eae 
p = m1 > (« +(° | (4) 


There are many cases where the velocity potential 
is more convenient. The method to be described in the 
following can be applied to the velocity potential by a 
treatment similar to that used for the stream function. 
For this purpose the velocity potential is introduced, 
defining, from Eq. (2), 


1 0 1 Od 


hee hg OB 


Substituting Eq. (8) into Eq. (1), the differential equa- 
tion for ¢ is obtained as follows: 


o ( hg 4) ") ( ha =) , 
p +—is = 0 
Oa\ hada OB\ hg OB 


The variation of density is also given by Eq. (7). 
(3) PRINCIPLE OF THE METHOD 


Considering an ideal flow of compressible fluid be- 
tween two fixed boundaries, the principle and procedure 
of the calculation may be described as follows: 

As an inviscid fluid in contact with rigid boundaries 
flows along the fixed surfaces and has no component of 
normal velocity, the value of the stream function along 
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CALEULATION oF 
the rigid surfaces is constant. Choosing a group of 
curved lines for one component of the orthogonal curvi- 
linear coordinates represented, say, by 8 = constant, 
the curves 8 = 8, and B = 8, are coincident with the 
two rigid boundaries a and 6 as shown in Fig. | and ac- 
cordingly with the stream lines along these surfaces. 
The boundary conditions are expressed, for example, as 


follows: 
y = y. at B = B, (10a) 
y= pat B = B, (10b) 


The remaining curves are assumed of similar form to 
the bounding stream lines. Then more exact shape is to 
be determined further on in the solution. 

On such orthogonal curvilinear coordinates, there- 
fore, OW O@ should be nearly zero along the curves of 
8 = constant, while O¥/08 should be nearly constant 
along @ = constant, and accordingly each of the 
two terms in Eq. (6) will become vanishingly small 
separately. 

Under these conditions a good approximation may be 


obtained by putting 


“beget a ei 
OB\ph,dB/ 


where it must be noticed that Eq. (11) included only 
derivatives with respect to a single variable 8. Based 
on iterative methods, the convergent solution will be 
approached by calculating the correction terms dropped 
in Eq. (11) or by adjusting the coordinate curves that 
were first assumed. 

Rigorously speaking, the basic conception of this 
theory is expressed as follows: 

Referring to suitably chosen orthogonal curvilinear 
coordinates, it will be assumed that a partial differential 
equation can be reduced to an ordinary differential 
equation in a single variable, the remaining variables 
being regarded as parameters. Convergent solutions 
are to be obtained by correcting the parametric vari- 
ables or by correcting the coordinate axes successively. 

For the computational procedures, thus, two sorts of 
methods have been proposed. 
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Fic. 4. Distribution of mass flow. 
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Nomenclature relating to the undisturbed flows. 
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Fic. 8. Adjustment of the inlet stream lines. 
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(a) Method of Fixed Coordinates 


In this method the coordinate axes are left fixed as 
they are first chosen. Separating the terms in the 
parametric variable, Eq. (6) will be converted into an 
ordinary differential equation. 

The equation of mth order approximation, y,,, is 


°( l ove) o ( 1 hg OW,» ) : 
= — ( 
OB\p» hs OB Oa \pm-1he Oa , 


where the right-harid side represents the correction 


9 


term that is calculated from quantities of the former 
By putting Oy%/Oa = 0, we have the 
same expression as Eq. (11) for the first approximation, 


approximation. 


The convergent final solution will be obtained by suc- 
cessive substitutions of y,, through the integrations of 
Eq. (12), a brief description of which has been given in 
the previous paper'' with some general expressions. The 
calculations, however, are more or less complicated, 
since they contain a differential procedure in the course 
of correcting solutions. For practical purposes, the 
second method will be more convenient. 


b) Method of Adjusting Coordinates ‘Method of Flux 

Analysis) 

There has been presented an idea that the stream 
lines calculated from Eq. (11) may be chosen as the 
more appropriate coordinates for approaching to the 
convergent solution.'' Thus, another means of correct- 
ing solutions is by successive adjustment of the curvi- 
linear coordinates. In this case it must be noticed that 
the values of Ha and fg must also be adjusted at every 
stage of approximation. 

The differential equation for the mth order approxi- 
mation is expressed as follows: 


Pay (; Ren) Obm ) “e (13) 
OB», 1\Pm hs m 1 OBm 1 


Neglecting the subscripts for the sake of simplicity, 
integration of Eq. (13) is easily obtained. 


(1/p)(ha/hg)(Ow/O0B) = F(a) (14 
which is rewritten in the form 


hau = Fla) = he.atta (15) 


a-a 


or 


4 / ty = Ne. 2/he (15’) 


where the quantities with subscript a represent the 
appropriate boundary values (cf. Fig. 1). 

Approximate distribution of velocity ratio along the 
curvilinear coordinate a = constant is given by Eq. 
(15’) as the reciprocal ratio of the extension factor ha. 
The constant of integration, F(a) or u,,is not determined 
in this stage, but it is obtained from boundary condi- 
tions applied to the further integration of Eq. (14) as 
shown in the following: 
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CALCULATION OF 


Integrating Eq. (14), y is given by 


y-Yy= J puhgOB = | puon (16) 
Ba 0 


the constant of integration being determined by the 
boundary condition (10a) as y,, which may be set equal 
to zero without fear of losing any generality. 

In many cases, it will be more convenient to express 
Eq. (16) in the nondimensional form 


y [ pu. (") 
= 0 
poCol 0 polo l 


where the density is given by 


ol (y-1 
p -|1 ; (“) | (17) 
3 Co’? J 


Po - 


(16’) 


Substituting Eq. (15’) into Eq. (16’) with the as- 
sumed value of u,/co, W can be calculated as a function 
of u,/co. For the computation it will be convenient to 
use the numerical table of pu/ poco versus %/'Co. 

The value of u, or F(a) is determined from the 
boundary condition (10b) connected with the integra- 
tion of Eq. (16’) from the boundary a to 6, which is 


*\/] 
» it nN 

Yo _ y| e-* (18) 

Pot ol 0 Po Co l 


Generally, u, or F(a) cannot be solved in analytical 


given by 


form except for the case of incompressible flow, in which 


case u, is calculated in the form 


/ =v /I 
Wa _ Yo / Rss 5 n (19) 
Us pUll So ha i 


where LU’ represents a standard velocity. 

In the case of a compressible fluid, calculation of 
y, for several values of assumed u, enable us to deter- 
mine the boundary value “, by numerical or graphical 
interpolation. These procedures are schematically 
shown in Fig. 2. 

Expression of the values of , versus a or s gives the 
velocity distributions along the boundary a as shown in 
Fig. 3. 

The velocity distribution and, accordingly, the dis- 
tribution of stream function or mass flow along the co- 
ordinate of a = constant is also calculated from Eq. 
(15’) and Eq. (16’), respectively, by substituting the 
value of u,. As the latter represents the quantity of 
¥ or mass flow through the divided section lying between 
the rigid boundary and the point concerned, the 
positions where the stream lines should be passed can 
be determined, and therefore, the stream lines can be 
drawn as described in the following procedure: 

Dividing the total mass flow into several equal 
parts values of n/l 
corresponding to the stream lines of '/4, ' 
respectively, will be given as illustrated in Fig. 4. 
Marking these points along the coordinates of a = con- 


for example, four parts—the 
2, and * 4y 
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Fic. 11. Flow through cascade of modified Joukowski Airfoils 


(incompressible ). 


stant and connecting the adjacent positions belonging 
to the same y from point to point, stream lines will be 
drawn as shown in Fig. 5. 

The stream lines thus obtained will generally deviate 
from the coordinate curves of 6 = constant, which were 
first assumed to be nearly coincident with the stream 
lines. 

The second stage of continuing the computation is 
started by taking the calculated stream lines as a com- 
ponent of the orthogonal curvilinear coordinates—.e., 
8 = constant. Another component of a = constant 
will be determined as groups of curves which orthog- 
onally cross the curves of 8 = constant at every inter- 
section, which may be conveniently drawn in graphic- 
ally. 

The differential equation for the second approxima- 
tion referred to the new curvilinear coordinates is also 
expressed by Eq. (13) and is integrated in the form of 
Eq. (14) or Eq. (15’), the latter giving the distribution 
of velocity ratio by the reciprocal ratio of the adjusted 
he. 

It must be noted that in the limit process, as con- 
tained in Eq. (4) in the expressions for ha and hg, the 
accuracy of calculation will be diminished with the re- 
duction of the small division, da or 68. For practical 
purposes, therefore, finite differences will replace the 


derivatives of ha or hg, 


ha = 6s/6a, hg = 6n/6B (4’) 
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Fic. 12. Flow through turbine cascade D5111914 


(incompressible ). 


If, taking an appropriate difference da, the coor- 
(6a/2) and a + (éa/2) constant, 
respectively, are drawn, the coordinate differences be- 


dinate lines of a 


tween these two curves at any value of 6 will be equal to 
ba. (Fig. 6.) 
(15’) therefore will be replaced by the reciprocal ratio 


The velocity ratio expressed by Eq. 


of the are lengths divided by these two coordinate lines 
as follows: 


=" 


U/Uq = 4658,/6s (15") 








BASIC TRIANGLE GENERAL ELLIPSE 
(a) (b) 
Fic. 13 Elliptic camber. 
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The boundary values u, are determined by the proce. 
dure described previously, and thus the successive cal. 
culations are continued and repeated until the cop. 


vergent solution is reached. 


Dividing the stream lines with the same difference of 


6¢ = 
orthogonal curves to the stream lines, and thus flow 
patterns will be completed. 


PROCEDURE THE 
THROUGH A RECTILINEAR CASCADE 


(4) COMPUTATIONAL ON FLow 


(a) Factors Characterizing the Cascade Flow 


In order to get reasonable approximations in the floy 
through cascades, it will be necessary to insert corree- 
tions to the flow angles of the outlet stream. For this 


reason, some _ relationships between characteristic 
parameters in the cascade flow are introduced. 

A diagrammatic representation of the undisturbed 
velocities at points at infinite distances is presented in 
Fig. 7. 

The flow is governed by the continuity, momentum, 


and energy equations with isentropic change, respec- 


tively, 
PiWa = pxWe = (pV) (20) 
L = (pV )t(wy + wy) 2] 
wy" Y pi_ wv? Y pb» 99 
2 : oe 7 2 y— 1p. 7 
p2/ pr = (p2/ pi)’ = (€2/e1)” iia (23 


where L represents the force component in the direction 
of the cascade surface acting on the unit height of an air- 
foil. Denoting the circulation around an airfoil by I, 
the next relationship will be obtained by tracing the 
circuit in Fig. 7. 


lr = t(wy W 2) (24 


where the circulation is also detined by the contour 
integral along the surface of an airfoil as follows: 


lr = fue ds (25 


Substitution of Eq. (24) into Eq. (21) gives 


L = (pV)r (26 

It is obvious that this relationship similar to the 
Kutta-Joukowski theorem holds in the compressible 
cascade flow only in the particular direction tangential 


to the cascade surface. If a circumferentially arranged 


cascade can be approximately replaced by the rectilinear 
cascade, the coefficient of total torque acting on the 
blade height of 6r at the radius of r may be introduced, 
defining 


6Cr = 67/(1/2) pywei?- 7: 2arrér (27) 


where 67° is given by 


uds, equipotential lines can be drawn as the 
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67 = YLrér = 9 (pV )Trér (28) 
Hence, 
— Y T r 
6c 5, = = 2 
T Wail Wart 
= 2 [tan B; + (wy2/Wei) tan Be] (29) 


The corrections of the outlet flow angle are made by 
using the relationship between flow angles and circula- 


tion. 
Transforming Eq. (22) with Eqs. (23) and (20), the 


Wai \’' 
al 
— 

War 


(30) 


jollowing expression is obtained: 
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(31) 


where W,1/Wa2 1S given by Eq. (29) as a function of 


tan Bo, 


ae a tan By» ( . — tan a) (32) 
Wa2 | \Wat 
Under the given values of the angle and Mach Number 
of the inlet uniform flow, the relationship between 6» 
and [—1i.e., outlet flow angle and circulation—will be 
calculated from Eqs. (31) and (32). 

For the case of an incompressible fluid, expressions 
are simplified by putting wa = wW,2, 


tan Bo = (I'/Walt) — tan p; (33) 


b) Calculation of the Inlet Stream Lines 


Stream lines bounded by two fixed surfaces of the 
adjacent airfoils are easily obtained by the successive 
The 
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method, however, must be slightly changed on the 
stream lines at the inlet or outlet portions, where the co- 
ordinate axes of a = constant will be extended to an 
infinite distance. 

In the graphical analyses, as well as the numerical 
ones, it will be impossible to deal directly with the 
problems in the infinite field. It is usual, therefore, to 
introduce some practical assumptions that the flow 
reaches a uniform state corresponding to the point at 
infinity at some definite distance adequately taken from 
the cascade surfaces. 

An example of the inlet stream lines of 6 = constant 
are shown in Fig. 8, with the perpendicular coordinates 
of a = constant extended to the assumed uniform state. 
Since the magnitude and direction of the uniform flow 
is given and the distributions of the velocity ratio are 
computable from Eq. (15’) or (15”), it is evident that 
the velocity distributions along the coordinate a = 
constant can be calculated in this case. In other words, 
the integral constant u, is determined at this stage of 
the approximation. Proceeding with succeeding cal- 
culations, therefore, the corrections of the stream lines 
or coordinate axes can be performed with no difficulty. 
In these processes the zero stream line to the forward 
stagnation point is also determined, though its accuracy 
cannot be expected to be excellent. 

In many cases it is convenient to assume the form of 
mean stream line first for determining the zero stream 
line to start the calculations (cf. Fig. 9). 

It must be noticed, in any case, that the correspond- 
ing stream lines on the adjacent airfoils should have 
the same forms. This property will simplify the calcu- 
lations to a great extent. 


(c) Calculations of the Outlet Stream Line 


In the case of the outlet stream line, the principal 
parts of the calculations are the same as those of the 
inlet flows, but the procedures are more laborious, since 
neither the magnitude nor the directions of the outlet 
uniform flows are known beforehand. In order to de- 
termine these values, corrections are necessary as men- 
tioned previously. 

Assuming the directions of the uniform flows, in the 
first step, outlet stream lines can be calculated and de- 
termined in a similar way as the inlet stream lines. If 
the airfoil has a sharp trailing edge, leaving stream 
lines should be separated at this point correspond- 
ing to Joukowski’s assumption for the isolated air- 
foils. 

The calculations are then carried out through the 
whole field. When it reaches a convergent state, the 
solution of cascade flow under the assumed direction of 
the outlet stream will be obtained. Using the velocity 
distributions along the surface of the airfoil, the circula- 
tion around the airfoil is also calculated from Eq. (25) 


(cf. Fig. 10). 
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Substituting the value of I into Eq. (32), the correc. 
tion of the outlet flow angle will be obtained from Eq. 
(31) by numerical calculation. 

The next approximation will be started by adjusting 
the outlet flow direction, and the whole procedure pre. 
viously described is repeated until a convergent soly. 
tion is practically reached, where the correction to the 
outlet flow direction is regarded as negligibly small, 

As the correction is small in many cases, two or three 
iterations will be sufficient for practical purposes. 


(5) ILLUSTRATIVE EXAMPLES 


For the purpose of inspecting the accuracy of the 
method and of forming an estimate of its validity, com- 
parison is made with some analytical solutions. 

For an example of an exact solution, let us introduce 
the complete solution of the incompressible cascade 
flow,'® '? which was obtained by transforming the iso- 
lated Joukowski airfoils into cascade arrangements. A 
Joukowski transformation connected with the mapping 
function of the cascade of flat plate gives airfoil con- 
tours with no stagger, as shown, for example, in Fig. 
lla. Applying the method of ‘‘flux analysis’’ to these 
cascades of modified Joukowski airfoils, two cases of 
inlet flow angle—i.e., 8, = O° and 7.72 
lated, the velocity distributions and stream lines of 


are calcu- 


which are given in Fig. 11. 

In spite of the low accuracy of the method for the 
flow through a cascade of such a low solidity, the values 
are in fairly good agreement with the analytical solu- 
tions. The angle of the outlet uniform flow is also de- 
termined with sufficient accuracy for practical purposes 
(cf. Fig. 11a). 

Another examination is then made by comparing 
with the experimental results on the cascade of a turbine 
airfoil.'* The experiments were made in a cascade wind 
tunnel having adjustable walls in the downstream sec- 
tion. The outlet flow angle was determined by adjust- 
ing the movable walls so that the lift coefficient calcu- 
lated from the momentum change between undis- 
turbed flows coincided with that obtained by the inte- 
gration of the measured pressure distributions. 

Calculating the stream lines and pressure distribu- 
tions of the flow through this cascade, fairly good agree- 
ment is also shown between the analysis and observation 
as presented in Fig. 12, where the values of both circu- 
lations are nearly coincident. 

From the considerations on computational accuracies, 
it is expected that the presented method of analysis will 
be suited to practical calculations of the flow through 
cascade of arbitrary arrangement. 


(6) TURBINE CASCADE OF IST-SERIES 


In order to get the basic profiles suitable for starting 
the calculation of predicting velocity distributions, 4 
study has been made by the ‘‘flux analysis” of the char- 
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acteristics of flow through a systematic series cascades 


designed at the Institute of Science and Technology. 


(a) Mean Camber Line 


The chord line and the two tangents to the camber 
line at the leading and trailing edges form a triangle 
called the ‘‘basic triangle’ of cascade airfoils, the 
vertical angle of which might be supposed to have some 
relations with the turning angle of cascade flows. For 
the camber line, therefore, a smooth curve that is in- 
scribed in the basic triangle with appropriate vertical 
angle, as shown in Fig. 13a, may be desirable. 

For this purpose a series of general elliptic cambers is 
introduced, the basic form of which is reduced from the 
ellipse of mth order (cf. Fig. 13b). 


(x/a)”" + (y/b)”" = 1 (34) 
a quarter part of which is used for the camber line. 


The coordinates relating to the chord line and its normal 
are calculated by the following formulas: 


&/l = [a(a — x) + by]/(a® + 6°) (35 
n/l = jay — b(a — x)}/(a® + B*)f ‘ 
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Denoting the coordinates for the point of the maxi. 
mum camber Q, by x1, v; or &, m, the following relation. 
ship is easily verified : 


: _ _ o-—l/n 
x,/a = y,/b = 2 (36) 


Substituting Eq. (36) into Eq. (35), the position and 
magnitude of the maximum camber are, respectively, 
given by 


_ 6 §f 7 7 l ea | 
dealt tat‘, (1 (5) hel 
_ B. §, l a t eh. 

.* m = {2(;) ‘Ti+ ¢ | 


where « = b/a. It is convenient for practical computa- 
tions to calculate « and 7 for the designated value of » 
and k. 


The formula is given as follows: 


F 9Q—!1 ‘ 1 ‘ 
e= V(k+2°”" — 1)/(2-”" — k) 38 
If it is necessary, formulas giving m and e for the 


designated k and f are also introduced after brief compu 


tations. 
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The basic cambers thus obtained are then deformed 
by elongating the figures along the normal direction to 
the chord line so that the vertical angle or the maximum 
camber reaches the desired value. 

In some cases combinations of ellipses of different 
order divided at the point of the maximum camber are 
employed to get the effects of form factor in the front 


and rear parts separately. 


b) Thickness Distributions 


Separating the thickness distribution into two parts, 
polynomial expressions introduced in the study of iso- 
lated laminar airfoils * are used. They are given by 
the following expressions ahead of and behind the 
maximum thickness, respectively: 


h off, &\’ é é\?) 
= -<i Zo 7 Fol — (39 
i ‘(ae £) +af+n(2)} , 


for 0 < ¢// S m, 


g/ 


hej. g & \? & \*) 
77 + a(1—!)+af1 -§) +a{1-‘)' 


(40) 
form S$ ¢/l S 1 


Coefficients are determined so that the radius of the 


leading edge circle is equal to g(e//)*/; that, h = e/2 
and dh/d§ = 0 at £/1 = m; and that, h = ce and 
dh/d§ = —d,(e/l) at £/1 = 1. 


In the following examples constants are chosen as: 


m = 0.30, g = 1.1, d,= 1A (41) 


c) Blade Profiles and Their Arrangements 

Examples of camber line and thickness distribution 
are given in Fig. 14. 

The superposition is made by taking the thickness 
height perpendicularly to the mean camber line at the 
corresponding percentage locations along the chord line. 

The cascade arrangements of the airfoils are deter- 
mined by giving the setting angle—for example, the 
angle o between chord line and turbine axis—and pitch 
chord ratio, ¢// (cf. Fig. 7a). 

IST-cascade series are described by three groups of 


numbering—for example, in 


IST 24-2040 
1030-1111 
0 = 35°, t/l = 0.6) 


% 49 
By = Go”, M, 0.24 (sa) 


the first line represents the form of camber line; the 
second, the thickness distribution; and the third, the 
arrangement. Adding the flow conditions in the last 
line, the flows will be completely determined. 

In the numbering of the camber line the first two 
numbers represent the orders of the ellipses ahead of, 
and behind, the point of the maximum camber. re- 
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Fic. 25, Flow pattern (4). 


spectively, so the different numbers indicate a com- 
bined camber line. The two numbers immediately be- 
hind the hyphen represent the amount of the maximum 
camber and the last two the location of the maximum 
camber both as percentages of the chord length. 

The numbering of the thickness distribution con- 
sists of four groups of two numbers, each representing, 
in order, the maximum thickness e//; its position m, 
both as percentages of the chord length; value of 10g 
related to the leading edge radius; and value of 10d,, 
defining the trailing-edge inclination. The coefficient 
of trailing-edge thickness is chosen as c = 0.08 from 
practical considerations. 

TURBINE CASCADES OF IST- 
SERIES 


(7) FLrow THrouGu 


In the graphical, as well as numerical, calculations, 
summary figures will be given for systematic variations 
Examples will be presented of the 
Results 


of parameters. 
effect of various factors on the cascade flow. 
are given in the form of velocity distributions along the 
airfoil surface, variations of the circulation referring 
to each parameter, and some examples of flow pat- 
terns. The velocity distribution is evidently the most 
significant expression of the airfoil characteristics, since 
it gives not only the intensity of circulation but also 
the basic data for the resistance calculations in the 
boundary layers. The conditions that must be satisfied 
by the profile of high efficiency are that the maximum 
velocity be far from the leading edge for the purpose of 
delaying the transition in the boundary layer and that 
the region of high velocity should be restricted as much 
as possible in order to reduce the frictional resistance. 
These conditions might be accomplished to a consider- 
able extent by a profile of uniform velocity distribution. 
As it is difficult to attain these conditions on both sur- 
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faces of cascade airfoils, more attention is usually paid 
to the suction surface, where the higher velocity causes 
more resistance than that on the compression surface. 
Any increase of surface velocity along the suction side 
causes an increase of circulation and surface resistance 
at the same time; therefore, the most effective correc. 
tion may be to approach a uniform velocity distribution, 
On the compression side, low velocities are obviously 
desirable, but the flows are accelerated in the neighbor. 
hood of trailing edge by the interference in cascade ar. 
rangeinents. 


a) Effect of the Location of the Maximum Camber 


The location of the maximum camber controls mainly 
the position of the peak velocity, as shown in Fig. 1), 
In these examples, a profile of k = 30 per cent gives a 
peak point too far forward, and conversely the profile 
of k = 50 per cent gives maximum velocity toward the 


rear. 


(b) Effect of the Magnitude of the Maximum Camber 


The magnitude of the maximum camber seems to 
have an influence on the value of the surface velocity, 
especially over the rear part of the profile (cf. Fig. 16). 

The maximum camber has an influence on the value 
of circulation as presented in Fig. 16b, where I’ is shown 
in a quasilinear relationship with the magnitude of the 
maximum camber. 


(c) Order of the Elliptic Camber 


The ellipses of higher order give a camber line of 
greater curvature in the vicinity of the maximum cam- 
ber point, while those of lower order distribute the 
curvature over the airfoil. Effects on the velocity dis- 
tributions are more or less systematic but are not so 
great as for the other parameters. The case of equal 
order for front and rear parts of the airfoil is shown in 
Fig. 17 and that of different orders in Fig. 1S. 


d) Effect of the Maximum Thickness 


Three cases of e = 10, 15, and 20 per cent are com- 
pared in Fig. 19, where the effect of increasing thickness 
is evident in the increase of velocities along both sur- 
faces. Since the higher velocities cause increased fric- 
tional losses in real fluids, thick profiles will not be de- 
sirable for cascades of high efficiency, even if they have 
Such high 


velocities in the flow passages will reduce the per- 


possibilities of increasing the circulation. 
formance, especially for compressible fluids. 


(e) Influence of the Pitch-Chord Ratio 


The solidity influences the velocity distribution along 
the compression surface in the special form, as shown in 
Fig. 20, where increasing solidity results in increased 
velocities. So the surface friction will be greater in the 
cascade with narrower passages, and therefore the total 
losses obtained by multiplying the number of blades 
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Fic. 26. Flow pattern (5). 


will be increased. In approaching arrangements of 
higher solidity, the torque coefficient is also increased 
but perhaps to a somewhat lower degree than the case 
of total losses. It is supposed, therefore, that there 
may be an appropriate solidity for a high-efficiency 
cascade, which gives a considerable amount of torque 
with comparatively small losses. High efficiencies will 
be expected in the vicinity of ¢// = 0.6 ~ 0.7. 


f) Influence of the Mach Number 


Influence of the Mach Number is studied on the 
standard airfoil. The 
velocity ratio along the profile surface is shown in Fig. 
21. As is well known, the effect of compressibility ap- 
pears in the region of high velocity to the greatest ex- 
tent. The distributions of ratio are not 
greatly different in the low velocity regions of the front 
part. Remarkable variations are presented in the high- 
velocity region near the narrowest part of the flow path, 


calculated variation of the 


velocity 


where, according to the increase of the inlet Mach 
Number, the maximum velocity becomes extraordinarily 
high and, finally, it exceeds the sonic velocity. 

An example of the compressibility effects on the cir- 
culation and turning angles is also presented in Figs. 


21b and 21c, respectively. 


APPENDIX 


Several examples of flow patterns are presented in 
Figs. 22-27. 
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A Vector Method Approach to the Analysis 


of the Dynamic Lateral Stability of Aircraft 


L. STERNFIELD* 


Langley Aeronautical Laboratory, NACA 


SUMMARY 


A description of a vector method approach to the analysis of 
dynamic stability problems is presented. In addition to supply- 
ing information on the damping of the dutch-roll oscillation, the 
phase relation between roll to yaw, and the amplitude ratio of 
roll to yaw, the method makes possible an excellent physical 
visualization of the contribution of the various stability deriva 
tives and mass characteristics to the overall motion of the air 
plane. The application of the vector method to several problems 
is briefly discussed. 


DISCUSSION 


ly THE SPRING OF 1951, two British scientists visited 
the Langley Laboratory of the NACA, and in the 
course of reporting on several phases of research related 
to dynamic stability of aircraft, they mentioned that 
Dr. K. H. Doetsch, a German scientist now working 
in Britain, had developed a vector method analysis of 
dynamic stability problems. Unfortunately, the Brit- 
ish representatives were not familiar enough with the 
method to describe it in detail. Our curiosity about a 
so-called ‘‘vector method”’ stimulated thought on the 
problem, which resulted in the development of a method 
and the design of an analog computer, to be descrited 
Recently, W. O. Breuhaus, of Cornell 
Doetsch 

It ap- 


in this paper. 
Aeronautical Laboratory, visited with Dr. 
and published a paper! describing his method. 
pears that, in principle, the vector method developed 
at Langley is similar to Dr. Doetsch’s method. 

The fundamental principle involved in the method is 
well known from the solution of linear differential equa- 
tions with constant coefficients and will be illustrated 
by considering first a simple spring-mass system and 
then introducing a dashpot into the system. In Fig. | 
we have a sketch of the spring-mass system and the 
differential equation describing the system, mX + kx = 
0, where the two dots signify the second derivative of 
x with respect to time. One complimentary solution 
of the differential equation is x = 
the natural frequency of the system and ¢, is an arbi- 
trary constant. Thus, each term of the differential 
equation could be plotted as a vector, with the vector 
inoment or force due to displacement along the positive 
x axis of magnitude k and the vector moment or force 


C, sin wf, where w 1s 


Presented at the Aerodynamics Session, Annual Summer 
Meeting, IAS, Los Angeles, July 15-17, 1953. 


* Head, Stability Analysis Section. 


due to acceleration 180° out of phase with the displace 
ment of magnitude mw*. As indicated in the equation, 
the sum of the two terms is zero, or the resultant of the 
two vectors is zero. As shown on the plot, the velocity 
is 90° out of phase with the displacement. 

Fig. 2 shows the case of a dashpot added to the sys- 
tem. The differential equation becomes m¥ + dx + 
kx = 0. 
a damped sinusoidal oscillation, and the expressions 


One complimentary solution is x = ce"! sin wf, 


for the velocity and acceleration are given on Fig. 2. 
It should be noted in the analytical expression for the 
velocity that a phase angle is introduced which is a 
function of the damping and frequency of the system, 
« = tan“! (a/w), and the sine of the angle is directly re- 
lated to the damping of the system. 
multiplied by a factor of two is also present in the ex- 
Representing each term 


This phase angle 


pression for the acceleration. 
of the differential equation as a vector, we have once 
again the displacement plotted as a vector along the 
positive x axis of magnitude k, the velocity is 90° + « 
out of phase with the displacement of magnitude dw, 
and the acceleration is 180° + 2 out of phase with the 
displacement of magnitude ma *. As the equation in- 
dicates, the resultant of the three vectors is zero, shown 
by the dashed vectors on the figure. Thus, the im- 
portant property of the vector representation of a 
damped oscillatory system is that the velocity is out of 
phase with the displacement by 90° plus an angle that 
is directly proportional to the damping of the system 
and that the acceleration is out of phase with the dis- 
placement by 180° plus twice the same angle. 

Now the method is extended to two degrees of free- 
dom. Continuing with our example of mass-spring- 
dashpot systems, Fig. 3 shows the differential equations 
of the sketched coupled system with two degrees of 
freedom in x and y. The solutions for x and y are 
also given. Of interest in the two degree of freedom 
case is not only the period and damping of the oscilla- 
tion but also the phase relation between x and y de- 
noted by the phase angle @ in the solution for y and 
the amplitude ratio of y to x, the ratio c./c,, where 
c, and cs are the coefficients of the solutions in x and 
y, respectively. The vector plot shown on Fig. 3 
indicates how the phase relation between x and y 
and the amplitude ratio of y to x can be determined 
from the first equation if the damping, a, and the fre- 
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Fic. 1. Vector representation of mass-spring system. { 
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Fic. 2. Vector representation of mass-spring-dashpot system. 
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in a two degree of freedom system. 
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On the left-hand plot, the terms 
in x are represented. Note k; + ks, the displacement 
vector; (d; + ds)wo, the velocity vector that is 90° + 
the angle e« out of phase with the displacement vector; 


quency, w, are known. 


and mw", the acceleration vector located 180° + 2¢ 
out of phase with the displacement. Adding these 
vectors, we have the resultant vector X. Similarly, 
on the right-hand plot, the terms in y are plotted on the 
assumption that the phase angle, 6, is zero—that is, 
y like x is considered to be positive along the horizontal 
Also, the amplitude ratio ¢2/c; is 
that is, the vectors have the 


axis to the right. 
assumed to be unity 
magnitude of the constants in the equation, kz and dz, 
and are not multiplied by the ratio c./c;. Since the 
right-hand side of the equation is zero, the sum of the 
two resultant vectors shown on the plot as X and Y 
must be zero. It is obvious that the condition required 
for the sum to be zero is that the vectors be diametri- 
cally opposite each other and of the same magnitude— 
that is, ¥ should assume the position of }’ and be of 
the same length as Y’. Therefore, the phase angle 
that was initially assumed to be zero is actually the 
angle through which Y is rotated to assume the position 
of Y’, and the amplitude ratio, initially assumed to be 
unity, is the ratio of the resultant vectors }’/ Y. 

If only the frequency is known, the damping of the 
system, as well as the phase relation and amplitude ratio 
of x to y, can be determined, working with both equa- 
tions simultaneously. Fig. 4 illustrates the procedure. 
The equations of motion of the coupled system are 
shown once again. Initially, we assume that the phase 
angle between x and y is zero, that the amplitude ratio 
of y to x is unity, and that the damping angle, «, is 
zero—that is, the system is assumed to have zero damp- 
ing. The upper plot on Fig. 4 contains the resultant 
vectors of X and Y for Eq. (I) and Eq. (II) when e = 0. 
In both equations, (I) and (II), the sum of the result- 
ants must be zero and the phase and amplitude ratio 
of Y to X must be the same. Considering only the 
phase relation between X and Y, it is apparent from 
the figure that, in order for Y and X to be diametrically 
opposite each other, Y would have to be rotated in the 
counterclockwise direction in Eq. (I), whereas in Eq. 
(II), Y must be rotated in the clockwise direction. 
Thus one may conclude that the assumption of « = 0, 
zero damping for the system, is not correct. For the 
damping angle equal to 2°, « = 2°, the resultant vec- 
tors for Eqs. (I) and (II) are plotted on the middle por- 
tion of the figure. As in the upper plot, here too, the 
phase angle required to make the sum of the resultant 
vectors zero in Eq. (I) is different from the phase 
angle required to make the sum of the resultant vectors 
zero in Eq. (II). does not correspond to 
the correct damping in the system. For e = 4°, the 
resultant vectors of Eqs. (I) and (II) are plotted in the 
lower region of the figure. In this case, rotating Y to 
Y’ about 10° counterclockwise in both equations and 
also multiplying Y by the same constant, about 0.85, 
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VECTOR METHOD 
in both equations would satisfy the condition that the 
sum of the resultants is zero in each equation. The 
damping of the system is therefore determined directly 
from the damping angle, e = 4°, the phase angle be- 
tween and X is about 10°, and the amplitude ratio 
of Y to X is about 0.85. 

In Fig. 4, three discrete values of « were selected to 
show the variation in magnitude and direction of the 
resultant vectors as € varies. It appeared possible to 
design an analog computer or “‘vector machine’ to 
enable one to introduce a continuous variation of € in- 
to the system and note visually the corresponding 
changes in the magnitude and direction of the resultant 
vectors. <A ‘‘vector machine’’ was designed and built at 
Langley,* based on the vector method, and was used 
for the problem of analyzing the lateral stability of an 
airplane. For the sake of simplicity, the problem in- 
volved two degrees of freedom of airplane motion, roll, 
and yaw, with the assumption that y, the angle of 
yaw, is equal to minus £, the sideslip angle. Fig. 5 
shows the nondimensional linearized equations of mo- 
tion used in the analysis. For several airplanes inves- 
tigated, the results obtained from the two degree of 
freedom analysis on the damping of the dutch-roll oscil- 
lation, the phase angle between roll and yaw, and the 
amplitude ratio of roll to yaw were in very good agree- 
ment with the results obtained from a three degree of 
freedom analysis if, when using only two equations of 

* Particular thanks are due F. S. Malvestuto, Jr., who assisted 
in the design of the machine, and P. R. Thomas and I. B. Bloxom, 
of the 7- by 10-ft. Tunnels Shop, for the excellent job they did in 


constructing the machine 





I. m,x+(d,+d,)x+(k+k5)x-d, y-kpy=0 
i. -dox- Kox +Moy+doy +koy =O 


— I EQUATION II 
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Determivation of the damping, amplitude ratio, and 
phase angle in a two degree of freedom system. 
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An approximate form of the nondimensionalized lateral 


equations of motion assuming y = —8 


Fic. 5. 


motion, the values of C,, and C,, inserted in the yaw 
equation are determined from the expressions shown on 
the figure. The undamped natural frequency of the 
system is well approximated by the expression given on 
the figure, where (Cx,)4 and (Cx,)4 are the airplane's 
values. A detailed description of the machine is 
given in the Appendix. 

Application of the 
Problems.—As noted in the 
method makes possible the physical visualization of the 
contribution of the stability derivatives to the damping 
of the dutch-roll oscillation. If, therefore, an airplane 
requires additional damping, which may be introduced 
through artificial damping devices, this property of the 
method enables the selection of the most effective type 
of auxiliary damper where a perfect damper is as- 
sumed—that is, a damper without any lag. 

The vector method may also be applied to the prob- 
lem of determining the stability derivatives of an air- 


Vector Method to Several Other 


Appendix, the vector 
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Fic. 6. Determination of the derivatives C;, and Cig from flight 


tests. 
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VECTOR METHOD 
plane if the mass characteristics, as well as the period, 
damping, phase relation, and amplitude ratio, are avail- 
able from flight tests of the airplane. Fig. 6 illustrates 
the approach to this problem. The roll equation is 
selected as an example. Since the three stability deriv- 
atives, C),, C;,, and Cig are not known and only two 
conditions are known—namely, the phase between 
roll and yaw and the amplitude ratio of roll 
it is impossible to uniquely determine the 
However, if one of the derivatives 


to yaw 
three derivatives. 
is estimated, the other two are uniquely determined. 
Assuming that C,, is estimated, the vectors correspond- 
ing to C,, and A,” are easily plotted. Since y is out 
of phase with ¢ by a known angle, the coordinate sys- 
tem for representation of the y terms is rotated through 
this phase angle as indicated on the figure. The A,, 
vector multiplied by the amplitude ratio ¥/¢ and the 
resultant y vector required to satisfy the equation are 
plotted on the rotated coordinate system. The direc- 
tion of the C;, and Cig vectors 1s known, and therefore 
the magnitudes are uniquely determined so that the 
sum of the three vectors, A,., C;,, and Ci multiplied 


by the amplitude ratio y/¢, gives the dashed resultant 
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tives should be estimated in order to obtain the most 
accurate determination of the other two derivatives 
requires additional research. 

Another application of the vector method arises in 
automatic stabilization problems, where it is generally 
necessary to know the frequency response of the system. 
Through use of the vector machine, data required to 
calculate frequency responses may be obtained with 
ease and rapidity. A change in frequency is intro- 
duced by simply changing the lengths of the vectors. 
The phase angles and amplitude ratios for each degree 
of freedom, read directly from the vector machine, are 
the data needed to compute the frequency response of 


the system. 


APPENDIX 


A schematic sketch of the mechanical analog machine 
designed and constructed at Langley is shown in Fig. 7. 
The analog machine consists of two tables with two 
decks on each one of the tables. On the yaw table, the 
lower deck represents the yawing moment due to an- 
gular displacement, velocity, and acceleration in yaw, 
and the upper deck represents the yawing moment due 
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An exploded view of the roll table 
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to rolling velocity and rolling acceleration. On the 
roll table, the lower deck represents the rolling mo- 
ment due to rolling velocity and rolling acceleration, and 
the upper deck represents the rolling moment due to 
angular displacement, velocity, and acceleration in 
yaw. 

The exploded view of the yaw table is shown in Fig. 
8 and corresponds to the case of zero damping of 
the dutch-roll oscillation—that is, the vector moments 
due to velocity and acceleration are 90° and 180°, re- 
spectively, out of phase with the vector moment due 
to displacement. The vector moments corresponding 
to the stability derivatives and mass characteristics 
are represented by metal rods. On the lower deck, one 
horizontal rod corresponds to the directional stability de- 
rivative, C,,, whereas the other horizontal rod corre- 
sponds to the moment of inertia in yaw, —2uA,*«,*. The 
rod normal to Crug and —2uA "wy" is the damping-in yaw 
derivative, —(1/2)C,,,a). The vector attached to the 
Crs arm through the vertical pin is the resultant of the 
three vectors. This pin moves in a groove cut into the 
resultant vector. Each one of these arms is scaled in 
tenths of inches so that values of the derivatives and 
inertia characteristics may be set accurately through 
thumbscrews. The magnitude of the resultant vector 
is the value shown on its scale at the point of intersec- 
tion of the vertical pin with the resultant vector. 

The vectors — (1/2)C,,,wy) and —2uA "wy" are geared in 
a manner to cause the acceleration vector —2uA-°a," 
to rotate through twice the angle through which the ver- 
tical vector —(1/2)C,,w9 rotates. Rotation of the vec- 
tors is performed by rotating the shaft. The angular 
rotation of the shaft causes the vectors to assume a par- 
ticular position relative to each other corresponding to a 
system with a specific amount of damping of the dutch- 
roll oscillation. The displacement vector Cn, is not 
affected by the shaft rotation but always lies in a hori- 
zontal position. 

The upper deck of the yaw table is made of lucite. 
On this deck, the vertical arm corresponds to the cross- 
coupling derivative —(1/2)C,,wo, and the horizontal arm 
is the product of inertia term —2yuA,,0*. The result- 
ant of these two vectors is easily placed into position 
and rests on the small protruding pin that fits into the 
hollowed groove of the resultant vector and moves 
within the groove as the arms rotate. The angular 
motion of the vectors is geared in a two to one ratio—the 


acceleration vector —2uA;,w)? rotating through twice 
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the angle of the velocity vector —(1/2)C,,o9. One 
set of gears controls the rotation of the vectors on the 
table. As the vectors rotate, the magnitude and direc- 
tion of the resultant vector changes and the angular ro- 
tation of the resultant is read on a fine angular scale 
engraved on the lucite. 

Since the derivative C,, and the product of inertia 
factor A,, may, for presently designed aircraft, have 
either positive or negative signs, this deck was so de- 
signed that the two arms can be located in any one of 
the four quadrants. As they are presently set, they 
correspond to negative values of C,, and A,,. If K,, 
were positive, the horizontal arm would be shifted 180° 
to its present position. Similarly, if C, 
the vertical arm would be shifted 180 

The roll table was designed and constructed similar 


» Were positive, 


An exploded view 
of the roll table is shown in Fig. 9. On this table, the 
lower deck contains the metal rods corresponding to 


to the yaw table just described. 


the damping-in-roll derivative —(1/2)C,,wo and the 
moment of inertia in roll —2uA,°a)? and their resultant 
vector. The metal rods on the upper deck correspond 
to the effective dihedral derivative Cig» the cross- 
coupling derivative —(1/2)C;,a, the product of iner- 
tia term —2yuK,.wy?, and their resultant vector. Here 


again, one set of gears controls the rotation of all the 
vectors on the table. 

The solution of a problem for a specific airplane re- 
quires working with both equations simultaneously 
The undamped natural frequency of the system w 
is approximated by the expression presented in Fig. 5. 
The solution obtained is the damping of the dutch-roll 
oscillation, the phase angle between roll and yaw, and 
the amplitude ratio of roll to yaw. The damping is 
determined from the angle through which either the 
velocity vector must be rotated or half the angle 
through which the acceleration vector must be rotated 
in order that the phase angle between the resultant 
vector on the yaw table is the same as the phase angle 
between the resultant vectors on the roll table. Also, 
the amplitude ratio of the resultant vectors on one 
table must be the same as the amplitude ratio of the 
resultant vectors on the other table. 
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Results of Some Base Pressure Experiments 
at Intermediate Reynolds Numbers with 


A = 


2.84 


L. L. KAVANAUT 
Unwersity of California al Berkeley 


ABSTRACT 


Base pressure data are presented for a cone-cylinder model at 
Mach Number 2.84 and Reynolds Numbers from 45,000 to 400, 
00. These data verify the maximum in base pressure coef 
ficient (po/p) 


tween the completely laminar and turbulent flow régimes 


predicted by Crocco and Lees which occurs be 
Sting 
length requirements for this model are summarized reflecting the 
growth of the critical wake region with decreasing Reynolds 
Number. Results from varying the base pressure orifice loca 
tion and sting diameter are also given 


INTRODUCTION 


i COMPLEX FLOW MECHANISM in the wake which 
determines the base pressure has been the subject 
of analysis by Crocco and Lees! for the case of a two- 
dimensional airfoil with a blunt trailing edge. Their 
analysis qualitatively describes the various effects of 
Reynolds Number (based on total length and noted 
base (ratio of 


as Re,) on the pressure coefficient 


base free-stream static pres- 
noted as p/p). 
maximum of the base pressure coefficient at 
nolds Number intermediate to a completely laminar and 


This maximum results from 


measured pressure to 
The analysis predicts a 


a Rey- 


sure and 


completely turbulent flow. 
transition occurring in the wake region immediately 
behind the base. Assuming that the base pressure 
for a blunt-based axisymmetric model responds to 
changes in Reynolds Number in the same qualitative 
manner as the two-dimensional airfoil considered by 
Crocco and Lees, the maximum base pressure coef- 
ficient should provide the approximate upper bound in 
Reynolds Number for the completely laminar flow 
range. 

No experimental data are available which give direct 
verification of the predicted maximum in base pressure 
coefficient. Bogdonoff? has reported the effect of 
Reynolds Numbers between 0.7 X 10° and 17.6 X 10° 
on the base pressure of a representative body of rev- 
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olution at Mach Number 2.95. His data at the lowest 
Reynolds Number show the base pressure coefficient 
increasing sharply with decreasing Reynolds Num- 
ber. Preliminary experiments at Berkeley with similar 
models and Mach Number at Reynolds Numbers near 
1,000 have shown decreasing coefficients with decreasing 
Reynolds Number. 


In order to obtain experimental data through this 


transition region, a series of tests at intermediate 
Reynolds Numbers was performed utilizing the super 
sonic wind tunnel at the University of California at Los 
Angeles (U.C.L.A.). 


report the results of these tests. 


It is the purpose of this paper to 


TEST OBJECTIVES 


The major experimental objectives of these tests 


were to determine the effect of Reynolds Number over 
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Fic. 1. Base pressure models—dimensional information. 
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Fic. 2. Model installation in the wind tunnel. 
the range 45,000 to 400,000 on the base pressure of a 
particular aerodynamic configuration. 

It was, of course, necessary that final results be void 
of all aerodynamic interferences. For most supersonic 
wind-tunnel experiments upon a fore-body of revolu- 
tion without boat-tail, interferences are eliminated by 
preventing the model bow wave reflection off the tunnel 
boundaries from interacting with the model and by 
supporting the model from the rear by a sting. Ex- 
periments concerned with the base region of the model 
require that further attention be given to the influence 
of the supporting sting on the flow in this region. 
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Studies of some such influences have been reported by 
Perkins* and Chapman?‘ for Reynolds Numbers greater 
than abcut 0.4 * 10°. 
flow were anticipated for the intermediate Reynolds 


Significant changes in the wake 
Number range of these tests. Therefore, it was desijr- 
able to include a limited study of the influence of sting 
length and sting diameter on the base pressure. 

A further test objective was to study the uniformity 
of the base pressure distribution. Preliminary data 
obtained at Berkeley for models in completely laminar 
flow indicate large radial variations of pressure. Ex- 
periments appearing in the literature for models with 
turbulent wakes or boundary layers have not included 
investigations of this problem, and the base pressure 
has been taken to be uniform." * ° 
that appreciable nonuniformities might occur in the 


It was possible 


base pressure distribution for the intermediate Rey- 
nolds Number range of this experiment. 


WIND TUNNEL 


The U.C.L.A. 3- by 3-in. supersonic wind tunnel 
is the open-circuit type and operates continuously with 
atmospheric stagnation pressure. It is powered by 
two 60-hp. reciprocating Chicago Pneumatic vacuum 
pumps in parallel evacuating a reservoir downstream 
of the diffuser section. By introducing a stagnation 
regulating valve, the free-stream Reynolds Number 
was made variable from 78,000 to 205,000 per in. witha 


nominal Mach Number of 2.S4. 
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Fic. 3. 


The variation of base pressure coefficient with Reynolds Number. 
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RESULTS OF SOME BASE 

The nozzle performance over this extended range of 
operation was studied through impact tube surveys. 
These surveys succeeded in identifying a usable test 
rhombus with a maximum Mach Number variation of 
+().02. 
assumed isentropic, and the Mach Numbers and static 


For the purposes of this paper, the nozzle was 


pressures were determined by impact and stagnation 
pressures. It is worth noting here that this assumption 
may not be too inaccurate in view of the comparison 
of present data with Bogdonoff’s results for similar 
models as shown in Fig. 3. The ordinate is extremely 
sensitive to free-stream static pressure p by definition 
and to free-stream Mach Number . through its in- 
fluence on the base pressure p)." * The nozzle calibra- 
tion for atmospheric stagnation pressure is given in 
reference 6. 


MOoDELS 


The models tested were geometrically similar to 
models used by Bogdonoff—i.e., 60° apex angle right 
circular cones with cylindrical afterbodies 2.50 cal. 
long. All models were fabricated from combinations of 
aluminum and brass and were slightly polished to ap- 
proximately a 25-microin. finish. Three basic model 
sizes were used with various radial locations of the 
base pressure sampling orifices and sting diameters. 

A control was provided for adjusting the effective 
sting length. This control was accomplished by mov- 
apex angle cone-cylinder along the 
When placed at the base 


ing a coaxial 60 
sting to the desired location. 
of the model, this disturbance cone provided mechani- 
cal support for the model during starting and stopping 
of the tunnel. 

Fig. 1 is a schematic drawing of the model configura- 
tion and a listing of pertinent dimensions for the three 
2 is a photograph of a typical 


models tested. Fig. 


model and support installation in the tunnel. 


SPECIAL INSTRUMENTATION 


Model base pressures were measured by means of a 
precision oil manometer,’ and the wind-tunnel stagna- 
tion pressure was measured by means of a mercury 
micromanometer. Both manometers were referenced 
to a partial vacuum of approximately 9 microns (0.009 
mm. Hg). A coincident beam schlieren system was 
used for flow visualization. 


iSXPERIMENTAL RESULTS 


Base pressure coefficients for a cone-cylinder model 
at .\J = 2.84 were obtained for Reynolds Numbers from 
$5,000 to 400,000. 
shown in Fig. 3. A curve indicating some of Bogdonoif’s 


The final results for dry air are 


results at higher Reynolds Numbers is incorporated for 
comparison. 

The effect on the base pressure of varying the sting 
length at constant Mach and Reynolds Numbers for a 
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Fic. 4. Typical effect of variable sting length on the base pres- 


sure coefficient at constant Reynolds Number 


typical run is reproduced in Fig. 4. The accumulation 
of all such data is presented in Fig. 5 as the critical 
sting length in calibers plotted versus Reynolds Num 
ber. The critical sting length is defined as the mini 
mum sting length possible for obtaining a base pres- 
sure within 0.5 per cent of that obtained with an “‘in- 
finite’ length sting. 

Results from tests of the two sting diameters and 
various orifice locations indicated no significant in- 
fluence of these variables on the base pressure coef- 
ficient over the Reynolds Number range tested. All 


data were within +0.7 per cent of the reported values. 
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Fic. 5. Effect of Reynolds Number on the critical sting length. 
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WAKE THROAT DISCUSSION 


The existence of a maximum in the base pressure co- 






EXPANSION TRAILING 


efficient as predicted by Crocco and Lees was verified as 


















SHOCK pec eee Fe ; : ; 
MIXING REGION seen in Fig. 3. Their analysis attributes this maximum 
mainly to the variation with Reynolds Number of the 
M,p mixing rate between the outside nearly isentropic 
— — ee stream and the critical wake region (see Fig. 6) im 
= + 4! < a : ie , Thi : ee 
r tS 2 mediately aft of the base. This change in the mixing 
_ we a rate is due to the transition moving into the critical] 
\N wake region from downstream of the wake throat as 
: the Reynolds Number increases from the completely | 
BOUNDARY laminar flow range. ae 
LAYER Crocco and Lees describe the action of the wake suct 
throat or “critical point’”’ in isolating the critical wake lea 
CRITICAL WAKE REGION eg . ¢ ease disr 
ey est : ; region from influences farther downstream. This wake 
Fic. 6. Exaggerated sketch of flow over axisymmetric model. ; : : dist 
throat was present with the sting-supported models of ' 
clos 
these tests,* and its changing position with Reynolds is { 
Numbers significantly influenced the sting length re- upp 


quirements as shown in Fig. 5. It is seen that the 
critical sting length increases with decreasing Reyn- 
olds Numbers. 

Fig. 7 presents two schlieren photographs of one 
model illustrating the critical wake region in a dis- T 
turbed and undisturbed state. 





Consideration has been given to testing two-dimen- wit 
sional models in this intermediate range of Reynolds uns 
Numbers where the critical wake region is so extensive. the 
These models would almost certainly have to incorpo- atti 
rate unique devices such as fences extendi.;; from the The 
model leading edge downstream past the wak 2 throat or in V 
some system of tunnel boundary-layer control in order the 
to keep the turbulent fluctuations in the tunnel-wall nos 
boundary layer from disturbing the critical wake re- son 
gion. Unless such precautions are taken, the critical In 1 
wake region, under conditions of decreasing Reynolds the 
Number, will not become free of turbulence as it would be 
0. OISTURBED WAKE, R/p = 4.0 in the free-flight sense described by Crocco and Lees. sup 
This would result in excessively high base pressures and this 
| would prevent the drop in base pressure coefficient that 
which is characteristic only of the completely laminar flow 
model at these Reynolds Numbers.' T 
i not 
CONCLUSIONS te 
(1) Base pressure data have been obtained for a in re 
cone-cylinder model at Mach Number 2.84 and Reyn- is tl 
olds Numbers from 45,000 to 400,000. A maximum pres 
base pressure coefficient (p,/p) of 0.75 occurs at a Reyn- Proj 
olds Number (Re,) of 140,000. The occurrence of pres 
such a maximum is in accord with the qualitative the- The 
resu 


ory presented by Crocco and Lees.' 
(2) Variation in size of the critical wake region re- dire: 





sults in significant increases of the critical sting length com 
’ cura 

*It should be noted that this phenomenon has been pre- I 
SOLU! 


viously observed by experimenters at higher Reynolds Numbers 


b. UNDISTURBED WAKE 5 R/ D = 7.2 than those of these tests, as well as in the preliminary experi- : 
- . . : 2nts at Berkeley in the laminar flow region 
Fic. 7 Schlieren photographs illustrating the effect of sting ae Sere . 


length on the critical wake region at Rez, = 112,000. (Continued on page 274) 
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The Flow over a Flat Plate with a Small 
Angle of Attack at Mach Numper | 


GOTTFRIED GUDERLEY* 
WADC, Wright-Patterson Atr Force Base 


ABSTRACT 


rhe pressure distribution over a flat plate at Mach Number | 
At the 


suction side of the plate, a cavitation region will arise close to the 


is determined with the usual transonic approximations 


leading edge; the influence of this region on the flow field can be 
disregarded if the angle of attack is small enough. The pressure 
distribution on the pressure side is obtained as a relatively simple 
closed expression; on the suction side the pressure distribution 
characteristics, but it can be 


the method of 


by a 


is found by 


ipproximated similarly simple analytical  expres- 


sion 


INTRODUCTION 


—— PRESENT PAPER determines the pressure dis- 
tribution over a flat plate in two-dimensional flow 
with a free-stream Mach Number 1. The only other 
unsymmetric flow field that has been determined up to 
the present is the flow over a wedge with an angle of 
attack that is small in comparison with the wedge angle.° 
The present problem can be considered as the extreme 
in which the angle of attack is large in comparison with 
the opening angle of the wedge. The flow around the 
nose of the profile causes, in either example, a super- 
sonic region to arise at the suction side of the profile. 
In the case of the wedge with a small angle of attack, 
the effect of this local region on the lift can, however, 
be neglected; in the present case of a flat plate, the 
supersonic region extends over the whole profile. In 
this regard the present problem shows a phenomenon 
that did not occur in previous examples of transonic 


flow patterns. 


The method by which the solution will be found is 
not quite systematic. Certain properties of the solu- 
tion can be seen from the particular solutions discussed 
in reference 1. 
is then found by an attempt to construct a simple ex- 


these 


The expression for the stream function 
pression that has these properties. However, 
properties are not sufficient to determine such an ex- 
pression uniquely, so this solution must be verified. 
The advantage of this procedure is the simplicity of the 
result. Previously, the author has carried out a more 
direct approach, which was successful; however, the 
computations were very cumbersome, were less ac- 
curate, and did not reveal the simple properties of the 
solution. 
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(1) THE BOUNDARY VALUE PROBLEM 


Fig. | shows the flow over an inclined plate at Mach 
1. From the leading and the trailing edges, sonic lines 
extend to infinity. At either edge one has a Meyer ex 
pansion. The Meyer expansion at the trailing edge 
occurs because a supersonic pressure prevails at the 
upper side of the plate. Close to the leading edge one 
finds a stagnation point. Stream lines that lie above 
the stagnation stream line “‘spill over’’ the leading edge 
The flow over the sharp leading edge that thus arises, is 
accomplished by a Meyer expansion. 


Among the waves of the Meyer expansion there is one 
in each case, the so-called limiting Mach wave, which 
separates the Mach waves that reach the sonic line 
from those that do not reach it. For the subsonic part 
of the flow field, only those portions of the supersonic 
region which lie upstream of the limiting Mach wave 
are of importance. 

By the Meyer expansion at the leading edge, an ex 
tremely low pressure (neglecting the boundary laver, 
even a vacuum) is reached; the stream lines separate 
from the surface, even if no boundary layer is admitted. 
Behind this low pressure region one has a rapid recom 
pression and a flow reattachment. 

It is practically impossible to prescribe here a realistic 
Fortun itely, this boundary con 
De 


boundary condition. 
dition is unessential for small angles of attack. 
tails will be found later 


The hodograph formulation of the problem is shown 
in Fig. 2. Corresponding points in Figs. | and 2 have 


the same notation. All stream lines start at point 0, 


which corresponds to the free-stream velocity. The 
lower side of the plate maps into the straight line 
DCB. Point D corresponds to the leading edge; 


point C, to the stagnation point; point B, to the trailing 
edge. The figure shows also a few stream lines, in par 
ticular the one that runs to the stagnation point 
Stream lines that lie in the vicinity of the stagnation 
stream line will, in their further course, run along the 
side of the plate. Those that lie below the stagnation 
stream line will follow the portion CB; those above the 
stagnation stream line will run along CD. Thus they 
will have practically a reversed flow direction, and only 
later will they turn back into the original direction. 
The limiting Mach waves are represented by the char 


acteristics EO and AO. 
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In this boundary value problem, simplifications that 
arise for small angles of attack may be carried out 
For the important part of the flow field, this process 
introduces the so-called transonic simplifications of 
the hodograph equation. 

Let w and @ be polar coordinates in the hodograph 
plane, and let w* denote the sonic velocity. y may be 
the ratio of the specific heats. We introduce a new 
independent variable 7 defined by 
(1) 


n = [(w — w*)/w*] (y + 1) 


Usually the transonic simplifications are obtained, by 
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SCIENCES—APRIL, 1954 


where / is a parameter. The hodograph equation js 
then simplified by retaining only the lowest powers in 
In the 
present example the quantity ¢ has a geometrical sig- 


t. Thus Eq. (4) (see next section) will result. 


nificance—viz., one can set 


t= 0 3 


where 6) is the angle of attack. Small values of / will 


then correspond to small angles of attack. The slope 
of the plate will then be given by 6 = —1. As @& tends 
to zero, the portion DE of the boundary will move to in- 


finity in the 7, 6-plane. The characteristics are then 













introducing semicubic parabola. 
of (2a This hodograph filed is shown in Fig. 3. The origin 
7 * <a) of the n, 6-system is shifted to point B. Then the free- 
6 = 6/t (2b) stream velocity will lie at @ = 1. Along ABC one has 
a ? 7 — 
i, 
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tion is the boundary condition y = 0. At infinity, a condi- 
vers in tion must be derived from the original requirement aa 
In the that Y = Oalong CDE. The reader will suspect that, 
‘al sig- | because of the limiting process that this line moves to CHARACTERISTIC. Lae anncremene 
infinity, the exact shape of this line—i.e., the exact “Ys 
3 boundary condition prescribed in the cavitation re- - 
gion— will be unimportant. This is a significant sim- D \\ 
twill plification brought about by the assumption that the Se ; 
> Slope angle of attack is small. The requirement that at zero / | 
)tends | velocity there must be a branch point of the stream ' Fe 
> to in- lines is, however, essential, and it is not automatically o [ ‘6 1, / CHARACTERISTIC 
> then fulfilled. How this condition expresses itself in the ' Cc 7 
; limiting case of small values of the angle of attack 4 will \ 8: 
origin be shown after suitable particular solutions have been ; stasnanen / 
e free- introduced. STRE AMLINE* 
ne has = ee 
—— (2) PARTICULAR SOLUTIONS FOR y 
Fic. 2 


The limiting process that 4 tends to zero leads to the 
following simplified hodograph equation 

Yi, — m5 = 0 (4) 

Here y is the stream function. 

need certain the particular solutions that have been 

studied in detail in reference 1. Introducing in Eq. (4) 


new variables, 


For the following we 


where A is an arbitrary parameter. The function G 
can be expressed by hypergeometric series (see refer 
Naturally G depends also upon \; therefore the 
Further 


ence 1). 
notation G(£, A) will be used occasionally. 
more one has for each value of \ two linearly inde 
pendent solutions G. They may be denoted by G* and 
G*, the superscripts a and s indicating that the corre- 
sponding solutions are antisymmetric or symmetric 


p= —7° + (9/4)0 (5a) 
eatie i a with respect to the line 6 = 0 (§ = —-« Formulas 
= 7/[(3/2)0]” (5b) ¢ ‘a , % , . . > ; ( . 
. for the G* and G* are found in Eqs. (18) and (19) of 
see Fig. 4), one is led to the differential equation eaacchaaaea l. 

We now proceed to determine the particular solutions 
ir On, ‘a = eu _2 :, (7) which will be utilized for the representation of the 
¢ Ve — 58 — Eve = 90, + > rh (0) solutions y in differe ‘gions of the 7, 6-plane j 
¢ 2 2 erent regions of the 7, #-plane. In Fig. 

3, the line p = 9 4—that is, the line p = const. which 


Particular solutions of this differential equation are 
found by the hypothesis 


2 


Y = p~(l 12)+(11/3) VAG (E) (7) 


In the vicinity of point B, the solution is subject to the boundary condition 
Because of the condition along BC, the solution is antisymmetric with respect to 


runs through point 0-——has. been drawn. The regions 
for which p < 9/4 and p > 9/4 will be denoted as inner 


or outer region, respectively. 


¥Y = 0 along AB and along BC. 
& = —o. In the vicinity of the 


line AB (that is, & = 1), the arising expression is then represented by 
2 2 Bes 1 
jr(5 va sine (CVA +, 
2 A, 1 d 
aT ae 


bo l™ 
A 
| 
fo 


see Eq. (18d) of reference 1; for 1 — — < 1 the hypergeometric functions occurring in this equation can be re- 


placed by 1]. A form more suitable for discussion is obtained if we insert here the definition for p [Eq. (3a) ], 
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r( Vx) sin x ( Va ) 
3 3 4 


nm, 
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ra 
>) 
‘Seal 
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Here V/ \ can assume either sign. The condition y = 0 
along AB requires that the first term of the bracket be 
zero. Besides, the exponent of the term 7/é in front 
must be positive, otherwise point A (the trailing edge) 
would not map into the finite part of the physical plane. 

Thus one finds as possible values of +/ \ in the vicinity 
of point B 


Vr = [(9/4) + 3h] baw @ig... (9) 


i.e., the desired function will be represented by a super- 


position of the particular solutions 


y = p~’” 


As can be seen from Eq. (18) in reference 1, in the vi- 


rm 3h]*} (10) 





-G*é,[(9/4) + 


cinity of the line € = —© these solutions are repre- 


sented by 


y = p®/)t* (—$) (11a) 
and in the vicinity of the sonic line (& = 0) one obtains 
ae: 
| =m 9 I & 
iad, 
”o 
[r ( +h)-0 ( - i) 
5 6 
} 2 l 
rye) | 
Fe 3 ys 
é (11b) 
1 ] 
r ( +h)-T (- —h 
3 6 
Ihe last expression can be written in the form 
(5/6) + A} 
J m , Jo ys (5/6) 
= - ; 2 —)- _ 
4 [* r[(5/3) + h] 
Pays) 
rlG 6) + h| 
_ (7/6) | 
Ga i alee dads (I 1c) 
8 ee) 
(4/3) 


If one studies the behavior of the functions G for high 
real values of /X (using, e.g., the method of Langer*) 
one finds, that everywhere in the range — © <& <0 
the function G is of the order of magnitude of its value 


at the sonic line We mention this for the following 


reason: Later y will be represented along the sonic 


line by a development in powers of p. This develop- 
ment will converge along the sonic line (£ = Q) in the 
We shall then obtain 
from this development a representation of yY for values 
of & 


Pa 


region from point B to point 0. 


~ () by inserting for each power of p, which 
occurs in the development of y along the sonic line, the 
corresponding particular solution for ¥. Because of 
the property of G just mentioned, this representation 
will then converge (even absolutely) for all values of 
p smaller than the value that p assumes at point 0—i.e., 
for p < 9/4. Also, the derivatives will converge abso- 
lutely; therefore within the region p < 9/4 (within the 
inner region), the differential equation, Eq. (4), for y is 
fulfilled. 

The particular solutions for the outer region can be 
found by the following reasoning: The part of the 
hodograph field in which we are interested extends to 
the limiting characteristic EO of Fig. 2. But one can 
extend the solution beyond this line by prescribing the 
boundary condition y = 0 along the continuation of 
the characteristic DE beyond point E 
tinuation occurs physically if one constructs the flow 


Such a con- 


field downstream of the limiting characteristic. We 
consider now the hodograph field up to the character- 
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THE FLOW 
7, -plane. Because of the boundary condition y = 0 
along DEE’, no singularity will propagate along E’A. 
Therefore the second term in Eq. (8) must vanish. 
Furthermore, the solution must vanish as p tends to 
infinity because of the condition along DE in the 7, 0- 
Therefore +/ must be smaller than '/,; [ef. 
These conditions yield the same values of 
A derivation of the same re- 


plane. 
Eq. (S) |. 
\ as in the inner region. 
sult by another point of view is found in references | 
and 8. Clearly, the exact shape of the boundary DEE’ 
does not matter after the transition to small values of 


4, has been made. 


As a result of these considerations, one finds that 
for 9 4 the solution can be represented by a superposi- 


tion of the expressions 


y= p lO". Gre, (9/4) + 3h]?} (12) 

We must still consider the condition that the stagna- 
tion point is a branch point of the stream lines. In 
the 7, #-system the stagnation point will correspond to 
7 = —©. Inspecting the functions G occurring in 
Eq. (12), one finds that only for h = 0 is the function 
G different from zero everywhere throughout the sub- 
sonic region (except naturally along the line & = —@), 
all other solutions have lines — = const. along which G 
and therefore also y vanish. For these other solutions 
infinity represents a branch point of the stream lines. 
At infinity, the terms of the expressions (12) with the 
lowest value of / will prevail over the others and thus 
determine the character of the branch point. Thus the 
stagnation point will be a branch point of the stream 
lines only if the outer representation does not contain 
the term 

p”*-G[E,(9/4)?] 

Thus the outer representation will be given as a super- 


position of the terms 


¥ = p (OVO +h. Galel(21/4) + 3h]?} 


(hk = 01,2...) (18a) 
Along £ = — ~ they are represented by 
y = (—§)—3/2. p-l/ete (13b) 
in the vicinity of the sonic line, one then has 
| r{((11/6) + h] 
jc, es 3 9! M(11/6) 
¥= a 
“ r[(s/3) + A] 
r(8/3) 
r((13/6) + A 
4 aa l'(13/6) 
{ P £- 2 Sr ySrEES (13¢c) 
64 r[ or 
(7/3) 


In the outer region (p > 9/4), the expression y for — # 
0 will be constructed from the expression for y along the 
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sonic line (£ = 0) in a manner similar to that used in 


the inner region. 


Let us now study the behavior of the solutions in the 
vicinity of point 0. One remarks that the differential 
equation (4) is invariant against a shift of the coordi- 
nates in the §-direction. We assume momentarily that 
the origin of the 7,6-system lies at point 0 and that the 
quantities p and & are defined within this system by 
Eqs. (5). They will then be denoted by j and —. We 
try to find particular solutions (7), by which the solu- 
tion of the present boundary value problem can be rep- 
resented in the vicinity of point 0. The conditions 
imposed here are that, along the limiting character- 
istics AO and EO, no singularities propagate and that 
the limiting characteristics lie in the finite part of the 
physical plane. that 
must propagate along these characteristics originates 
from the fact that at point A and E there is no singu 
In the vicinity of 


The condition no singularities 


larity in the boundary conditions. 
point 0, symmetric, as well as antisymmetric, solutions 
may occur. 

To discuss these conditions, we write the particular 
solutions for y in the vicinity of the lines AO or EO 
( = 1) in the form of Eq. (8). (This form corre 
sponds to antisymmetric solutions; the treatment of 
symmetric solutions is similar.) If the second term 
in Eq. (S) is different from zero, then for positive values 
of the exponent a singularity will propagate along the 
characteristics; for negative values of the exponent, the 
line will map to infinity. Therefore one obtains, as gen 
eral requirement, that the second term in Eq. (S 
vanishes. One finds the following values of y/ 
for antisymmetric solutions (11 4) + 3h, (7/4) + 3h, 
(3/4) + 3h 
solutions (1/4) + 3A, (5/4) +. 3h, 
8 i 


for symmetric 
—(3/4) + 3h (h = 


As discussed in reference 1 on page 50 (although with 
a somewhat different notation), not all of the negative 
values of YX listed above but only YX —3 4 and 


VA = —9/4, will occur, since otherwise the vicinity 
of point 0 would map into several sheets of the physical 
plane. The particular solution with y\ = —3 4 has 


no physical significance if it occurs as the only singu 
larity at point 0, but superimposed to the particular 
solution with 4/\ = —9/4 it shows at a far distance 


from the body the influence of its unsymmetry. 


Which terms will occur in a development that re; re- 
sents y along the sonic line? _In the original 7, #-system 
with the origin at point B, point 0 corresponds to 6 
1. Then one finds from Eqs. (7) and (5) that the de 


velopment around point 0 has the form 
y ~ (6 —1) 1/6) +(2/3)- VA 


In the same arrangement as in the previous table one 


finds, as possible exponents of (@ — 1), 
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(5/3) + 2h, 1 + 2h, — (5/3) + 2h 


2h, (2/3) + 2h, — (2/3) + 2h 


Among the corresponding particular solutions for y, 
those with the exponents (2/3) + 2h and (5/3) + 2h 
happen to give y = 0 along the sonic line; but y, is 
different from zero. For particular solutions with the 
exponents 24 and | + 2h, one has y ¥ 0, ¥; = 0. 

The particular solutions valid in the vicinity of point 
0 represent y along the sonic line for 6 < 1, as well as 
for@ > 1. The connection between the expressions in 
these regions is given by the fact that in the 7, 6-plane 
the particular solutions under consideration are either 
symmetric or antisymmetric with respect to the line 
6 = 1. Following this idea, one finds that the ex- 
pression y for 6 > 1 is found from expression y for 0 < 1 
by taking in any case the real root of the fractional 
power of (9 — 1). (One should be aware of the fact that 
this rule does not hold for more general values of ~/X.) 

Summarizing this, one can say that, along the sonic 
line in the vicinity of point 0, Y can be represented in 


the form 


We shall furthermore utilize the particular solutions (7) for imaginary values of YX. 
lie in the question of convergence of a Fourier integral formed with these particular solutions. 


view their asymptotic behavior for high values of | VX 


be derived from the differential equation for G with a standard method (for instance reference 5). 
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y = (6 — 1)~”* P,(6 — 1) + P26 — 1) (14 


where P, and P, (and later P;, Ps, or simply P) denote 


suitable power series. 
From the point of view of Tricomi’s uniqueness and 


existence theorem, two different singularities— namely, 


y=p 


are admitted at point 0. The boundary condition is 
Y = 0 along the entire boundary. 


these singular particular solutions are unknown 


The coefficients of 
One 
of the arbitrary parameters that thus appear in the 
solution is unessential, since a factor, with which the 
hodograph solution is multiplied, changes only the scale 
of the physical plane. The other parameter is deter- 
mined by the condition that the stagnation point is a 
branch point of the stream lines. The proper value of 
this parameter will appear in the future computations 


quite automatically. 


Our main interest will 
From this point of 
is of special interest. Asymptotic representations can 
One finds 


9) 
: 3 ; VIA 
G(&, A) ~ r (—£)~*/* (41 — 6°)" igh v(é) (A < 0, |A} > 1) (15a) 
s s s 9 s 
2 r 5) 
where 
v(t) = 2-tg-"(—#) 7°" (15b) 


~ 


The presentation is valid for & < 0—1. 
expressions for y will occur: 


»., for the subsonic region with the exception of the sonic line. 


The following 


y =p cos [(V/|\!/3) log p]-G(é,A) (16a 
y = p | sin [(V//\1/3) log p]-G2(é,A) (16b) 
Near the line & = — ~, according to the definition of G*(é,A), they are represented by 
y = (—£)7*?-p7'/*- cos [(+/|A|/3) log p] (17a 
y= (—£)7°? p7 sin [(o/A|/3) log p] (17b 
The asymptotic expressions for y are found by inserting into Eqs. (16) the above representations for G: 
3 1/12 settee VIA : VIA 
y~ — . 7s dee | — £") *-c0s ( log 6) sinh | : v (&) (18a 
2V/ |N| 3 3 
3 ~1/12 ial ace , \ 
i ~ —— 2p ; al W4q- e*) °~ -emi ( log p)-sinh [ v(&) (18b) 
2V |A| L 
Introducing a complex variable ¢ by 
¢ = log p+ wp (19) 
one can write the last expressions as 
9 | | 
o = 9 a 1 . NV r 
i op M2 g- 0/4 7. — 89)". Tp sin ( ; +) (20a) 
2v || 3 
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' VIA . 
y~ -p ' al cae} £3) ".Im cos ( = .¢) 20b 


>) j 


(3) THE CONSTRUCTION OF THE SOLUTION OF THE BOUNDARY VALUE PROBLEM 


The derivation of the solution given in the following will contain certain tentative elements and therefore cannot 
claim to be applicable in general. However it may be useful in similar cases. The solution will subsequently be 


verified. 

We try to find a representation of the function along the sonic line. In the previous paragraph we found the fol- 
lowing properties: (a) In the vicinity of point B, y is represented by a power series that proceeds in powers 
po” *'*". Because of relation (5) and since along the sonic line 7 = 0, one has 


y = 6°? P,(6?) 
b) For values of 6 > 1 one finds similarly 
y = 67" P,y6-?) 
c) The development at @ = 1 is given.in Eq. (14) as 


y == (A _— 1) -P\(6 — 1) + P.(6 — 1) 


Since the inner and outer representations proceed in terms of @°, we introduce 


= 2 (21 
and then can write the above conditions as 
y=2 P3(z) (fors < 1) (22 
y=! Pils! (fors> 1 (23 
y = (s — 1) *-Pi(s — 1) + Pals — 1) (in the vicinity of zs | (24) 


In the last expression one must always use the real roots of the fractional powers. From the requirements of con- 
vergence in the 7, @-plane, one can conclude that Eq. (22) converges for z < 1, Eq. (23) converges for z > 1, and Eq. 
(24) converges for |(z — 1)} < 1. 

These statements are certainly not sufficient to determine the function along the sonic line or to characterize it 
as a hypergeometric function. So it is only tentative if we assume that at the origin the function y is represented 


by 


yy = z’'*F(a,),c,2) 25) 

In the vicinity of the point z = 1, the solution will then be represented by two other hypergeometric series that 
: I I : yperg 

proceed in powers of (s — 1). Because of the conditions that one obtains for the transition through point 1, one 


finds that the outer solution must be represented by a superposition of the same hypergeometric functions. Thus 
it follows that the hypergeometric differential equation fulfilled by the function in the outer region is the same as 
that fulfilled by the inner representation, for a second-order differential equation is fully determined by two 


linearly independent solutions. Now, if for z = 0 one solution has the form of Eq. (25), then because of the rela- 
tion 
T(a)-T(d) I'(a)-T(b — a) , 

: - F(a,b,c,z) = : (—s)-*-F(a,1+a—c,lt+a-—b),2')+ 

l(c) I'(c — a) 


r'(b)-T(a — b) 


—z)~-F(b,1+6-—c,1+6-—a,27') (26) 
I'(c — dD) 


(cf. reference 6, page 257), the linearly independent solutions at infinity have the forms 
«7? 3.6 @. P(z-!) 
52/3 6. P(s-1) 

One finds from the first expression, by a comparison with Eq. (23), 


5/9 27) 


af 


The other linearly independent solution of the hypergeometric equation will not be needed. To investigate the 
behavior at point zs = 1, we use the formula 
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I'(c)-T(e — a — Db) 


F(a,b,c,z) = 2 : -F(a,b,1+a+6-—c¢,1-—2)+ 
l(c — a)-T(c — BD) 
(c)-T(a + 6b —c) — 
: ; — z)°*”.F(c —a,c —b,1+c¢—-—a-—b,1—s8) (28 
l(a) -T() 
(reference 6, page 252). This equation shows that assumption (25) yields, at z = 1, an expression of the desired 


form if the exponent of the second hypergeometric function has the right value—i.e., if 
c-—a-—-b= —5/3 
Thus, with the value of a found previously [Eq. (27) ], 
c—6= 5/6 29 


It is not necessary to make the same consideration for the outer representation since the outer solution fulfills the 
same differential equation and therefore has, at point z = 1, the same linearly independent solutions. 
Up to now we obtained as a likely expression for the inner region 


y = 2”F[5/2, b, b + (5/6), g] 
Then the particular solution, with the exponent a at z = ©, will, according to Eq. (26), be given by 


yp = Cz 11 oR 





3/2, (8/3) — 6, (7/2) — 6b, s-"] (30a 


~) - 1° 


Here Cisa constant that will be determined subsequently. 


Possible values of b are obtained if one examines how the inner and outer solutions can be joined at point z = 1, 


By applying Eq. (28) one obtains, for the inner representation, 


jI'[d + (5 6)]P(—5 : Clb + (5/6)]-T5/3) 
y : 2 -F(5/2), 6, 8/3, 1 — s) + = x 
(rb — (5/6)]-T5/6) (5 /2)-1T(d) 
Be i ie { 
(1 — 2) -Flb (3/3), 2/6 2/d, | Z| sla 
f 
and, for the outer representation, 
.(((7/2) — b|l(—5/3) 
y=C:s 16 J | | ; -F[5/2, (8/3) — b, 8/3, 1 e") 
( (1 — b)-T(5/6) 
>r[(7/2 b]-1T(5/3) see 7 E } 
(1 — g~!) ". Fl — 5, 5/6, —2/3,1 — 3s 31b 
Pta/2) ri(S 3) — dl { 
Here the second term may be written in a slightly diferent form, 
r{(7/2) — 6]-T(-5/3 
y C-3 6 Jo : F[5/2, (8/3) — b, 8/3, 1 — s—'] - 
{ r(1 — 6)-(5/6) 
r{(7/2) — b]-T(/3) ss 
: : ed —-z FI b, 5/6, —2/3, 1 g—!) (31¢ 
r(5/2)-T[(8/3) — d] f 
Since all hypergeometric functions (with the powers of s and z — 1 in front included), fulfill the same differential 


equation, it 1s self-evident that the power series with the exponents 0 and —5/3 in Eqs. (8la) and (3lc) agree, 


although this is not obvious from the form indicated. Therefore one needs only to equate the first coefficients 


I'[d f+ (9 (4) | - 1 — 9/3) Be all (4 2) — b| P(—5/3) si 
(92d 

r(é — (5/3)] 15/6) r(1 — 6) (5/6) 

rb t (5/6) ]-T(5/3) r((7/2) — b]-T5 3 
: —- —(C 32b 

'(5/2)-T(6) r(s/2)-T[(8/3) b], 
By a division, C is eliminated 
l(b) l(b (3/3 —T(8/3 b| (1 h) 


an one finds, because of the relation 


(zs) (1 — sg) = w/sin wz 
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that 
sin br = —sin [b — (5/3) |x 
Hence, b = (5/6) + k, where kis an integral number. Actually, only b = 5/6 will be used in the solution (this is a 
further assumption). C is then found from one of the Eqs. (32) as C = —1/2. Thus one has, as likely repre- 
sentation for y along the sonic line, 
y = 2"'F (5/2, 5/6, 5/3, s 33a) 
forz < 1 and 
y = —(1/2)-s7''’° F(5/2, 11/6, 8/3, s-) 33b) 


forz>1. 


(4) VERIFICATION OF THE SOLUTION 


From the expression for y along the sonic line given 
in Eqs. (33), one obtains solutions in the 7, 6-plane 
simply by replacing the various powers of z (or of 6) 
which occur in these series representations by the cor- 
responding particular solution (11) and (13) for @ < 1 
and 6 > 1, respectively. 

As mentioned before, these series representations will 
converge automatically inside and outside of the line 
p = 9/4. Thus it follows immediately that discon- 
tinuities can occur only along the line p = 9/4. 

It will be the object of this section to show that no 
discontinuities are present along this line except at 
point 0. This is done by deriving from Eqs. (33) an- 
other representation of y whose region of convergence 
includes the line p = 9/4 and thus shows that along 
this line no singularities are encountered. 


This solution will be verified in the next section. 


> 


From Eq. (33a) one obtains for the general term of 
the inner solution along the sonic line 


r[(5/2) + kA] T[(5/6) + A] 
_ (5/2) (5/6) i 
a 
r(5/3) +h] . 
one ‘Tih + 1) 
(5/3) 


Then one obtains by means of Eq. (11) for the general 


term in the development of y 


4 (2/3) +h g - 
g p <a “7 x 
‘ o 
9 : 
G2 E (; re 3h) | 


Thus, because of Eq. (lla), ¥ is represented in the vi- 


r[(5/2) + h] 
T'(5/2) 


(hk + 1) 
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cinity of & = — © by Hence, 
r(5/2) +h 8 , [4 j /2 
’ | ) | y= 28 (9 *(o 9) (1+ 34 
3/2 8 1/3 (4 - I'(5/2) 3 9 9 
ae hea ee Ore oe ‘ ; 
re 3° 9°) 4 r(h + 1) 
= , h=0 ’ The treatment of the solutions for p > 9/4 is similar 
h ° pan . — - 
(—) (| p) and yields the same result. The simplicity of the last 
9 expression is remarkable. 


From this result one can derive a representation of the solution, by means of the functions G" for negative values 
of A, which holds through the entire subsonic region. Assume for this purpose that in the following Fourier in- 


tegral the functions g,; and g. are known: 


a I ai wa I" VA VA VA 
: : 0 5) >) >) 


Then because of Eq. (16), Y can be expressed by 


S : i d |d ; d 
y= --2 ast 2) (¥. )-cos (¥, log p) Gre — |Al)-d ~ + 
vo Oo »o *) oO 
" WN d In 
/ a(™, ) sin (~ log p) Gr (&, — |A})d ~ (36) 
0 ” 0 >) 


The important question is the range of values £ in which this expression converges. It was mentioned previously 


that the convergence of such expressions is the same as the convergence of the corresponding asymptotic repre- 
Using Eqs. (16), (17), (18), (19), and (20), one obtains approximately 


e NY ) 

_ cos Og p 

4 ‘ - 2 /\Xr 3 : ap, 

¥y~ —2 aie, inked! Med ee “i 21 € ). : - ‘sinh |* ¥@) |-d0 AI/3) + 
0 5) V IAI /o p 


5 sin a og p -_ 
VIA 3) P VIA] An . 
29 : . = ¥ sinh v(E) |-d(V|Al/3) ¢ (37) 
0 3 VIX 3 3 ' 


With the definition of complex variable ¢ given in Eq. (19), the last 


sentations. 


where v(&) has been defined in Eq. (15b). 
expression can be written as 


VIA 


4 ‘ sin ( 3 J 
pm Sgt pote, = em) f a (© ). : -d(v/|x1/3) — 
| 0 5 / Vv A 
3 


. 3! ) 
£2 ; -d(v/ |A|/3) ¢ 
y) ( 3 ) VIN f 


or 


4 uy » 2 = ‘IN rn 4 / rN 
y~ ~.2 1/3, 5 l Sg Wty ¢3)! "am f 2 ¥ )-4 Ys f cos (™. ‘) dé + 
0 o oO ) 
[ ( _ ; { | (* ; +) ; 
£2 ) : ) sin 9 s lia 
/0 — ” ”o 


v 


In the region of uniform convergence, the sequence of the integrations can be interchanged, 
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The second integrations can be « 


The particular solutions occurring in Eqs. (36) and (37) 
The 
limit of convergence will be given by a line £ = const. or 
»(é In the complex ¢-plane, at least one 
singularity will lie at the line of convergence. The 
singularities of the integrand in Eq. (28) are the same as 
The line £ = — © cor- 


have an exponential behavior in the é-direction. 


v(t = const. 


the singularities of the integral. 
responds to vp(&) = 0. Obviously the singularities of 


the integrand in Eq. (38) lie at 


c log 4/9 + t2r(2n + 1 


where 7 is an integral number. Hence it follows that 
the development of Eq. (37), and therefore also Eq. 
36), converges for 


—mr< vt) < 


By Eq. (15b) we then find that in terms of £ the expres- 
sion converges for § < 0-—1.e., in the entire subsonic part 
of the hodograph plane. This shows that the inner 
and outer representation match along the entire line of 
convergence p = 9/4. 

As a check one can determine, from the original inner 
and outer representations, y; along the sonic line (it 
will also appear as a hypergeometric series) and see 
whether it fulfills the matching conditions at point 0. 
From the point of view of the construction of the solu- 
tion this check is unnecessary. 

PRESSURE DISTRIBUTION ON THE PRESSURE 
SIDE OF THE PLATE 


5) THE 


Putting the stream density at sonic speed equal to 1, 


one has the relation 


The x coordination is determined by 
x, = (y+ 1)? 
xo = (y+ 1)", 


From the definitions for p and & [Egs. (5)], one finds, 
for # = 0, 


0p/06 = 0 


(—9) 


‘arried out in a closed form because of Eqs. (35 


] 


a [ dc (<. ,3 


Remembering the definitions (2) and (3) for 7 and 6, one 


thus finds, from Eq. (34), 


no a(t) 


) x 
: Fhe, n 
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3 5 

7) Ay 
9) » 


“( 
J 
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Putting the origin of the x,y-system into the stagnation 
point and introducing a variable of integration 


T -n [(3/2)05] 


one obtains 


r= 1(? 
2 


The overall length of the plate is given by 


L=4(*=-) (° 4) / r'5 (1 


Here x and Z appear with a negative sign. 
to obtain the positive sign (which corresponds to Figs. 


dr 


~ 


In order 


1, 2, and 3), the expression must be multiplied with 


—1. By the substitution 


the last integral assumes the form 


-0 ; | 71 ; 
— se : 6 2/3 
/ ri) (1 + 1%) de = — = v (1 — 9) dv 
vo. oO 


This is a B-Integral (reference 6, page 212). One 


finds 


was ; 1 1(5/6)-T(5/3) 
7i5(] + 73)-*5 dr — = 
x 3 (5 2 
thus 
4 y+ 1 , Ft ) 71(5/6)-T(3/3) 
_ R = A, ; 
. 3 ( 2 2 (5/2) 


The lower part of Fig. [(3/2)0)]"’” versus 
x/L. Since the pressure coefficient is, in the present 


approximation, expressed by 


5 shows 7 
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Fic. 5. Pressure distributions over plate. 
— 1 —1/3 , 
oO = 2(y + 1) U] 
this figure can be interpreted as a representation of 
= 6,/2(y + 1)". [(8/2)60]"”* 


Occasionally one may be interested in the behavior of 
the flow at a far distance from the body. It can be char- 
acterized by y along the sonic line. The portion of the 
solution which is singular at point 0 is represented by 
the second term of Eq. (31a). One has, with b = 5/6, 
1'(5/3)? 


— 1 pom z)~° 3. FB = 6, 
P(5/2)0- (5/6) iad 


y = 2/3 . 


5/6, —2/3, 1 — 2) 


Remembering that 


one obtains as the behavior of y along the sonic line 
in the vicinity of point 0 


—5/3 1'(5/3)? I( ay ‘3 
y = 2 ik ee ee Eee a sae seh | == + 
1'(5/2) -T(5/6) 6 


e-97) 
12 0 
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This expression corresponds to a plate with a length 
given by Eq. (40). With the assumption made pre- 
viously that the stream density at sonic speed is 1, one 
then obtains for a plate of length Z as the behavior of 
y along the sonic line in the vicinity of point 0 


2237 t+ 1\~“ r/3) 
“a «i.g-'. (% ) ; 
v 4 2 1'(5/6)? 


3 | 4 — o\-7" 
4 j ~- ~ = 
( : ( A T 19 (: 4 


(6) THE PRESSURE DISTRIBUTION ON THE SUCTION 
SIDE OF THE PLATE 


At a first glance the task of determining the pressure 
distribution at the suction side of the plate appears 
rather hopeless, because of the presence of the cavi- 
tation region. However, the size of the cavitation 
region depends upon the angle of attack. If the angle 
of attack is made smaller, then the region of the plate 
which is influenced by the cavitation region becomes 
smaller, so that, for a sufficiently small angle of attack, 
the boundary condition in the cavitation region need 
to be introduced only in avery rough form. We have 
here the same reasoning as in the boundary condition 
for the pressure side of the plate. Mathematically, the 
diminishing importance of the cavitation region ex- 
presses itself by the fact that the map of its boundary 
moves to infinity in the 7, 6-plane as @ tends to 
zero. 


Moreover, the region of the flow field in which the 
transonic approximation cannot be applied shrinks with 
decreasing angle of attack; thus it is even possible to 
apply transonic approximations for the determination 
of the pressure distribution. 


Thus one obtains the following problem: By the 
previous investigations the flow field up to the limiting 
Mach wave AE has been determined ; one must find the 
flow downstream of it. One possible method of ap- 
proach is the method of characteristics. 


To reduce the amount of numerical work the solution 
is continued beyond the limiting Mach wave OE by 
analytical means up to the line 6 = 0. For 7 > 0 the 
location of this line in the physical plane can be found 
in a general manner, since the representation of the solu- 
tion along the line 6 = 0 is known for 7 < 0. One can 
use for this purpose Eqs. (3) and (43) of reference 4. 
In the continuation of the solution beyond the limiting 
Mach wave EO, the line 6 = 0 does not represent the 
surface of the flat plate (as one might suppose because 
of the value for 6 which prevails), but it represents a line 
in the flow field where the velocity vector is parallel to 
the plate surface. 


One difficulty must still be overcome: As one goes 
along the line 6 = 0 or along the limiting characteristic 
EO toward the leading edge, the value of 7 tends to 
infinity. Thus one does not have a starting condition 
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THE FLOW OVER A FLAT 


for the method characteristics. Moreover, the com- 
pression waves coming from the sonic line are reflected 
at the plate as compression waves and then coalesce to 
forma shock. It appears desirable to take into account 
the presence of this shock. Therefore the flow field in 
the vicinity of the leading edge has been determined 
analytically. For large values of p the term p 1s. 
[t, (21/4)?] predominates in the subsonic region, as 
well in its continuation beyond the line EO. If the 
preshock solution has this form, then the shock can be 
determined by the method shown in reference 9. The 
only difference between the problem of reference 9 and 
the present one is the value of X. 

The solution for the flow thus obtained has been used 
for the flow field in the vicinity of the leading edge, in a 
region where the neglect of the other terms contributes 
less than about 10 per cent of the entire solution. This 
region covers about 5 per cent of the length of the 
plate. 

The rest of the pressure distribution is then obtained 
by the method of characteristics. The result is shown 
in Fig. 5. 

The pressure distribution thus obtained is close to 
the expression which the formula for x valid in the sub- 
sonic region would yield if it is also applied in the super- 
sonic region—.e., if one expresses x as 


4 1 1/3 3 1/3 r= n/[ (3/2)60)*/s 
+ = 1(7 tt) (5% I x 


rid (1 — r*) 2.5 dr 


This result would have been obtained if one had ad- 
mitted in the supersonic region a multiply-covered 
region of the physical plane (a fold) instead of the shock. 
If such a fold and the change of state occurring in it is 
rather small, the replacement of a shock by the fold is 
quite harmless. In the present case, however, the 
value of 7 along a stream line which passes the fold 
would go to infinity (at the last Mach wave of the 
Meyer expansion) and then return to finite values. 
Under these circumstances a more careful computation 
appeared desirable. The result, however, shows that 
even here the effect of the shock on the pressure dis- 
tribution over the plate is rather small. 


(7) DISCUSSION OF THE RESULTS 


The pressure obtained by a Meyer expansion which 
deflects a sonic parallel flow by an angle 4 gives a quan- 
tity suitable for comparison with the pressures that 
arise in the present case. One obtains then 7/[(3/2) 
|"? = 1. This is also the value that would be ob- 
tained at the suction side of the plate if the flow were 
extended to infinity beyond the trailing edge (or, in 
other words, this is the asymptot of the curve 7/ 
[(3/2)0]*/ versus x on the suction side). The aver- 
age value of 7/[(3/2)@]* * on the suction side is 1.37. 
The same average on the pressure side is 0.81. The 
center of pressure of the upper side lies at 46 per cent; 
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Fic. 6. Comparison of the pressure distribution on the pres- 
sure side of the plate with pressure distributions at a wedge with 
the same slope of the side. 


the center of pressure at the pressure side, at 35 per 
cent; the center of pressure of the entire wing, at 42 
per cent of the chord. 

That the values of 7/[(3/2)0]** at the suction side 
are larger than | (than those obtained by a Meyer ex- 
pansion) is caused by the centrifugal forces of the air 
which ‘“‘spills’’ over. At the leading edge the pressure 
tends toinfinity. This is due to the approximation used. 
n at the leading edge is proportional to x~'°. The be- 
havior is the same for both sides of the plate. The 
power of x which occurs here is dictated by the lowest 
negative power of p which occurs in the outer solution 
i.e., it is caused by the expression p~ ''/°G*[(21/4)?]. 
We remember that originally the outer solution con- 
tained also a term with p °° and that this term was 
omitted because of the condition that the stagnation 
point is also a branch point of the stream lines. At 
first it may appear surprising that the behavior at the 
leading edge is connected to a condition at the stagna- 
tion point, however in the present case the stagnation 
point lies close to the leading edge and its distance from 
the leading edge tends to zero with the angle of attack. 
This explains its strong influence on the pressure dis- 
tribution at the leading edge. One sees that the pres- 
sure recovery in the vicinity of the leading edge need not 
be typical for the flow over a sharp edge in cases where 
the stagnation point lies at a distance from the 
edge. 

In Fig. 6 the pressure distribution on the pressure 
side has been compared with the pressure distribution 
at the side of a wedge with a half nose angle of %. The 
boundary conditions at the trailing edge are the same. 
At the leading edge a difference exists due to the pressure 
build up caused by the centrifugal forces of the air 
which spills over. The general character of the two 
curves is similar; the order of magnitude of the pressure 
is the same. 
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An Exact Solution for Two-Dimensional 
Linear Panel Flutter at Supersonic Speeds 
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INTRODUCTION 


S™ RAL WRITERS HAVE RECENTLY GIVEN attention to the cal 
culation of the flutter behavior of a thin panel, such 


In this type flutter, 


ads 
occurs in wing, tail, or fuselage coverings 
the basic structure supporting the panel acts essentially as if 
In 


ilyzing the phenomenon, the attempts to date have employed 


rigid, while the panel vibrates normal to its own plane an 
the generalized coordinate approach, using the lower normal 
modes of the panel as coordinates. 

In the present note, a theoretically exact solution is deduced 
for the supersonic flutter of two-dimensional panels of uniform 
uniform plates of infinite width—supported on 


thickness—1.e., 


ams running in the width direction. The supersonic flow is 
directed normally to the supporting beams and is past one surface 
of the plate. The fluid on the other side of the plate is stagnant, 
ind it is assumed that its pressure remains sensibly constant 
The flutter of a ‘‘flag,’’ with flow past both surfaces of the flag, 
in be derived from the present results by multiplying the pres 
sure P(x) by the factor 2 

The present study is restricted to a simply supported panel 
resting on two rigid beams. However, a straightforward ex- 
tension of the theory permits taking into account other types of 
panel end support, beams capable of elastic deflection, and the 
case of a panel resting on many supports 

In practice, the problem of interest is the panel of finite span, 
rather than the two-dimensional panel. However, it is hoped 
that the present analysis and its results will throw light on the 
mechanisms of this new type of flutter and will thus enable an 
evaluation to be made of the validity of using a Lagrangian ap 
proach to the finite-span panel problem, in which only the first 


few normal panel modes are used as coordinates 


EQUATIONS OF MOTION 


Consider a rectangular panel of great width and of length 2d, 


the panel being supported on rigid beams along the two wide 
edges (i.e., the leading and trailing edges). A supersonic stream 
with undisturbed velocity V and density p passes over one side 
of the panel, directed normally to the supported edges. The 
panel is of uniform thickness and composition everywhere 

The elastic behavior of the panel, depending on its rigidity, 
may fall within the elastic membrane province, may be described 
by the Bernoulli-Euler bending mechanism, or (in actuality) is a 


ombination of the two. For the present analysis, the integro- 


differential equations for pure membrane or for pure bending ac- 


tion only are deduced. It will be seen that the technique of 
solution is essentially the same in both instances. A combi 
nation of membrane and bending behavior can also be dealt with 


without difficulty 


Now let x be a dimensionless coordinate with origin at the lead 
ing edge of the panel and running in the stream direction to the 
trailing edge, where x has the value 2. Since a linear analysis is 
contemplated, during flutter the panel motion can be described 


x Jwl 


flutter frequency, and ¢ is the time 


as (x,t) = 2(x Here > is the panel semichord, w is the 


The dimensionless function 
(x) describes the shape of the panel during the flutter motion 
In the further work, since the exponential e7” occurs throughout 
in the usual fashion, it will be omitted from the mathematical 


development. 


Then, for a membrane of mass m per unit area and with a ten 
sion 7 per unit width in the x-direction, it is readily shown that 


the equation of motion is 


(T'/b) (d2s/dx?) + mw* (1 


Here P P(x) is the excess of pressure caused by the panel de 


flection, and 


<(x) is taken positive when the panel distorts 


into the stream. The factors 1/b and 6 in the first and second 


terms arise from the fact that x and s are nondimensional 


The fundamental natural frequency of the uniform membrane 


is given by w, = (2/2b)(7/m) ‘*, and so Eq. (1) can be put in the 
alternate form 
(d*z/dx*) + (3?/4) (R?/k,?)s = (6/T)P (2 
where k = wh/V and k, = w,b/V. The boundary conditions for 
the membrane problem are 
atx = Qandx =2;: s=0 3 


For the plate in bending, a similar derivation of the equation of 


motion leads to 


= (h3/D)P 4) 


(diz /dx*) — (mw*h*/D 


where D is the flexural rigidity of the plate and the other symbols 


have the same significance as before. The fundamental natural 


frequency for the simply supported plate is w, = (*/4b*) X 
(D/m) *, and so an alternate form for Eq. (4) is 
(d*z/dx*) — (91/16) (k?/k,»?)s = —(b3/D)P (5) 


where & and k, are as previously defined. The boundary condi 


tions for the plate problem are 


atx = Oandx=2; 32 d?:/dx? = 0 (6) 


With regard to the pressure P(x), an elegant expression has been 
given by Miles.' Converting 
» 


from —1 to +1, instead of the present 0 to 2) to the present defi 


from his x-coordinate (running 


nition for x, the expression for P is 
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pU x dw 4 = 
P= w(O)-G(x) + G(x — &)- + jkw }-dé (7) 
B 0 dt 


where p is the air density, U/ the speed, JJ the Mach number, and 


6? = M?—1. Also, 
G(x) =e KnMX. To xx) (8) 
w = U[(dz/dx) + jks] (9) 
where x = mk/8?, and w(0) is the value of the downwash w at 


x = 0. (Jo represents the Bessel function of the first kind.) 
Combining Eqs. (2) and (5) with Eq. (7) yields the integro- 
differential equations governing the membrane and plate flutter 


problems 


SOLUTION BY LAPLACE TRANSFORMS 


The solution for the membrane problem is now given. While 
the plate problem is not solved explicitly, it will be seen that the 
procedure is identical with that used for the membrane, both 
problems being amenable to solution by straightforward Laplace 
transform techniques. 

Combining Eqs. (2) and (7) and noting that s = 0 at x = 0, 
the integrodifferential equation to be solved is 
d*s 
dx? 

ion d*s ., a p 

G(x — &)- ~ ae Se — k*s}-dte (10) 
0 dt? dé f 
where 270 is the value of dz/dx at x = 0, a? = x?/4k,2, and after 
some manipulation it can be shown that 


(11) 


Now let 2(p) be the Laplace transform of 2(x) in the Carslaw 
Then, using the product theorem, the 


y = (1?/48)-(pb/m)-(1/k,?) 


and Jaeger notation. 
transformed Eq. (10) becomes, subject to the Eq. (3) boundary 
conditions, 

{p? + ak? — yG(p)-(p + jk)?} = = Zr0 (12) 


where G(p) is the Laplace transform of G(x). It can be shown 


that 
G (p) = ((p + jxM)? + x7] 7! (13) 
and so, after some algebraic manipulation, 
2(p)/zx0 = [A(p)/D(p)] + [B(p)/D(p)] -G(p) (14) 
where 
D(p) = (pb? + atk?)? [(p + jx Af)? + x2] — yp + jk)* (15) 
A(p) = (p? + ak?) [(p + jx M1)? + x?] (16) 
B(p) = yv(p + jk)? [(p + jx Al)? + x?] (17) 
Now let F(x) be the inverse transform of A(p)/D(p). Since 


A(p)/D(p) is the ratio of polynomials, the inverse transform is 
Similarly, let F:(x) be the inverse transform of 
It then 


readily deduced. 
B(p)/D(p); the latter is again a ratio of polynomials. 
follows that the solution for s is 


2(x)/t00 = F(x) + f, Fy(t)-e~# M(x 8) Jola(x — £)] dé 


x 
= F(x) + f, Fox — &)-e7~ IMME To(xt) dé (18) 
subject to the boundary condition s = 0 at x = 2. 


If a flutter solution exists, the manner of arriving at the flutter 
condition is now evident. For a given set of values for M/ and 
p(i.e., for a given altitude and Mach Number, so that the U’ value 
is also determined), flutter will occur for those k-values that re- 
sult ins = O atx = 2. The mode shape at flutter is then given 
by Eq. (18). 

If a flutter solution does not exist, there will be no finite k-values 
that lead to satisfaction of the x = 2 boundary condition. 


AERONAUT 


flutter mode is nonexistent. 
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The solution for the beam proceeds in analogous fashion, save 
that ‘wo boundary conditions must be satisfied at x = 2 in order 


for flutter to occur. 


ANALYSIS OF THE RESULTS 


Unfortunately, the analytic solutions to the membrane and 
plate problems are too complicated to permit many general con 
Infor 


mation for practical configurations can be gained only by numeri 


clusions to be drawn concerning the flutter mechanism 


cal calculations for specific cases; this is a drawback that evidently 
characterizes all phases of the flutter problem. 
Furthermore, the numerical calculations, even for a single 
panel, are lengthy. For each combination of panel parameters 
and p, M, and &, the roots of the Eq. (15) polynomial must be 
determined. The value of Eq. (18) must then be calculated at 
x = 2 to see whether a zero value for 2(2) results. If not, further 


k trials must be attempted. The evaluations of Eq. (18) entail 
considerable labor because of the convolution integrals involved 

Because of the lengthy character of the numerics, the present 
authors have been unable to immediately plan any detailed nu- 
merical studies. Hence, the reason for this partially complete 
presentation of the results. 

It is of interest to note certain conclusions applicable to the 
cases where y is either very small or very large. For y very 
small and approaching zero, the solution evidently approaches 
the case of simple panel vibrations in a vacuum—1.e., vibration 
in the panel normal modes. 

For y very large, it can be shown that flutter of both the mem- 
brane and the plate cannot occur. It can be shown by residue 
theory that, for y large and for the membrane, the Fq. (18) ex 
pression for s(x) reduces to 

of VS x 
= _ : ¢ 4 1 cV(x—8) o—jeME Tol xt) dé 19) 
=20 « ¥ “= J0 
which is of order e?*/y. It is evident from Eq. (19) that 2(2) 
cannot be zero for any finite k—i.e., for large y a physically feasible 
For the plate, a similar result is ob 
tained. 

It can also be shown without difficulty that flutter cannot occur 
for the membrane if the Mach Number ./ is sufficiently large. 
However, the lower boundary on MV has not been determined. 

It is thus seen that panel flutter does not exist for all configura- 
tions. The precise flutter determined 
through detailed numerical study. 


boundaries must be 
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INTRODUCTION 


N° SHAPES THEORETICALLY OPTIMIZED for minimum pressure 
drag for given fineness ratio have been derived by New- 
ton' 2? and von Karman.’ Because of the assumptions made in 
the derivations, one might expect the Newtonian shape to have 
minimum drag for hypersonic Mach Numbers and the von Kar- 


man shape to have minimum drag for high fineness ratios and 
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non, Save low supersonic Mach Numbers. However, for moderate fineness et Aa ae T T T 
: in = . ? . 
Mn order ratios and moderate supersonic Mach Numbers, one might expect | 

that other shapes are possible which have less drag than either kt = | 1 t 

the Newtonian or the von Karman. It is the purpose of this " ————- Second-order theory 

note to present the derivation of several such nose shapes and to 12 t t a anes £9s. (1) and (1a) 

| : s . é ah . . —-——-—£s. (1) and (1b) 
ane and compare their profiles and drags with those for the Newtonian ' 
eral con- and von Karman shapes. The derivation is accomplished by the 1 t t — 
Infor- yse of an equation for calculating pressure distributions devel- - 
numeri- oped by Bolton-Shaw and Zienkiewicz.* 10}—*-4 + +———++ + 
vidently ‘ 
. | 
SYMBOLS 9 1 4 | a a | } it 
i single a 7 
umeter B = aconstant, 1.117Mo 8 | = — a 
am 6 P . i 
é eens — " p 
must be ¢ = defined by Eq. (10) : . a 24 Y Fineness ratio, oy «3 
; Cp = pressure-drag coefficient, D/gory-" ° \ ae 
lated at P + cn | + + 
further D = pressure drag 
) entail Ip = defined by Eg. (4) 
volved M = Mach Number 
present p = static pressure 
led nu- pr = total pressure 
ymplete q = dynamic pressure ae 
x,y = coordinates with origin at vertex 
to the y’, y” = dy/dx and d*y/dx*, respectively 
(IN = one ¢ ’ 
y very ¥2/22 Suomen ratio 
roaches Z = By 
bration Subscripts 
0 = free-stream conditions 
* mem- 
; 1 = values at vertex 
residue 
2 = values at base 
18) ex 
. . "2 = 
DERIVATION OF NOSE SHAPES (oe) Approximete Newtonion, rong” 
(19) The pressure-distribution equation of Bolton-Shaw and Fig. .- Comporison of pressure distributions computed by various methods 

Zienkiewicz may be expressed as 
at 2(2) p/pPo = 1.014e" (1) 
easible where “3 yr ‘ae r — 7 . 
is ob 

—— 0.157y’ 

n = 1.405 Moy’ — 0.288 Mo Vy’? cet: leet cad 0 | | | F ‘ ; ; T 
occur V y"* — 2yy’ Second-order theory 
large. (1a) WHA Eqs. (1) ond (10) — 
e —-——.— £ gs. (1) ond (ib) 
igura- By derivation, this equation is valid for computing pressure dis- 10} + 4 
mined tributions over shapes (not too dissimilar from circular-arc ogives) 

for which the Mach Number divided by the fineness ratio is be- 9 + + + + + 

tween about 0.5 and 1. For fineness ratio 3 and Mach Numbers 

between 1.5 and 3, pressure distributions computed by Eqs. (1) 4 T if T T 
nis and (la) for an approximate Newtonian shape and for a von | 
Super- ee . , . } } 4 4 
. 351- Karman shape agree well with those computed by Van Dyke’s " 

second-order theory’ (see Fig. 1). (The approximate New- , | | | 

a hac . . , = wJlele.)/* wac te - T T 
tonian shape, defined by y y2(x/ 2) /*, was suggested by A. J. oa Fineness retie, gies 

Eggers, of the Ames Aeronautical Laboratory.) It is also ob- e z Ye i 

served in Fig. 1 that there is little or no loss in the accuracy of the 

pressure distributions computed by Eq. (1) if terms containing 4 | | | | | | 
} at y” in Eq. (1a) are omitted—that is, if Eq. (la) is reduced to | ge 

n = 1.117Moy’ — 0.157 (1b) - | 

This simplification greatly aids the derivation procedure and is 2 

justified on the basis of the pressure-distribution comparisons in 
oid Fig. 1. / 

The pressure drag of a given nose shape may be expressed in J 
coefficient form by a 
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e 
J 4 + + 4 4 Ee Eee * 
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where B = 1.117My Fig. |- Concluded 
5 








278 JOURNAL OF 


THE 





AERONAUTICAL 


SCIENCES APRIL, 


























SS bes 





Fineness ratio,  - =J3 
2 


—— von Korman, y= a Jb-zsin2g 



































. -14,- 4% 
g=cos a 
Moy’ 
——-—— Newtonian, y = ihe 
4y’ 
4, / ‘a 4 
x= —(Iny'+—+ —--+ 
=e a 4 UNy'* at ayaa! 
rr 7 = x 
a | | Approximate Newtonian, y = & (yi | 
| f- ' ny - i in on ———--——Wew shape, M,2 1.5 | 
( | |_| | ——---— New shape, Me= 3.0 | 
O a 2 | 4 ] 6 a 8 BY 1.0 


Fig. 2.- Comparisons of profiles of minimum drag noses for given length and base diameter. 


As becomes apparent later, it is convenient at this point to 
express the total pressure drag as the sum of the pressure drag on 
a finite area of infinite slope at the vertex: plus the pressure drag 
over the remaining profile. Taking the pressure drag at the ver- 
texas D, = pi’, the equation to be minimized may be written 


as 


qo ee 187 ; 
In = (Cyo— +1 y2? = 
Po 2.028 
0-187 Pr, "x2 
+1) 2 + ebu’ yy’ dx (4) 
2.028\ po 4 . 


Since the integrand function, e%4’yy’, is independent of x, the 
first integral of Euler’s equation yields 


Byy’? e®v’ = const. (5) 


If one of the extreme points (y; or y2) is on the x-axis and y’*e®u’ 


is not infinite, then no usable solution is obtained. This same 


condition arises in the Newtonian analysis. With ye fixed, how- 


ever, y; may be obtained by satisfying the terminal condition 


that® 
po... d Feo187 (Pr, 7 
v (eBy yy’) - a + ] > i > = V0 (6) 
loy dy | 2.028\ po loam 


By letting Z = By’, Z; = By,’, and Z, = By»’, Eq. (6) reduces to 


P ; e9- 157 Pr, 
e~'(Z, + 1) = + ] (7) 
1.014 \ po 


Pr, a =) ‘(Fe —1\7*? a 
bo = 5 6 ) (éa) 


Eq. (5) may be written as 


7 Zo\7—-9 Z > 
y = (Ze%y2e72)Z-2 e~7 (8) 


where 


Cie 
wince 


then, 
i stone : : Z? Aa 
x = — (LZ? y2e%?) InZ—Z+ - + 
: 6 4 18 
ee ae eer) ee , 
96 4 Z? Z* 
When x =x =0,9 = yp #0,Z = Z, # O, therefore, 
: ; ; Z;? Z3 Z\4 z l | 4 | 
C= -n7A+2Z- + — ee TE in ame eee ae 
} 18 96 Z, a: 2;' 
(10) 


Minimum-drag nose shapes can be computed by the following 
steps: 

(1) For given Wo, solve for Z,; by Eqs. (7) and (7a). 

(2) Solve for C by Eq. (10). 

(3) For given fineness ratio x2/2y2, solve for Z2 using Eq. (9). 

(4) For values of Z from Z, to Z2, solve for corresponding values 
of y/y2 using Eq. (8) and values of x/y2 or x/x2, using Eq. (9 


PROFILE AND PRESSURE-DRAG COMPARISONS 


Nose shapes resulting from the use of the previously derived 
equations are a function of both Mach Number and fineness 
ratio. Although the condition of vertex bluntness is imposed, 
this bluntness is generally small, y:/y2 being less than about 2 
per cent for fineness ratios of 3 and greater. For most pressure- 
drag calculations, the drag contribution of the blunt vertex is 
negligible. 

Nose profiles, computed for a fineness ratio of 3 and for Mach 
Numbers of 1.5 and 3, are compared with the Newtonian and von 
Karman profiles in Fig. 2. (Also shown for comparison is the 
three-quarter-power profile which closely 
Newtonian.) It is evident that the new shapes, although lying 
between the Newtonian and von Karman, are in best agreement 
In fact, there is little difference be- 


approximates the 


with the Newtonian shape. 
tween the new shape for My = 3 and the Newtonian shape. All 
of the shapes have infinite slopes at their vertices, and all, except 
the von Karman, have finite base slopes 
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Pressure-drag coefficients computed by second-order theory 


for the various shapes are compared in Table 1. 


TABLE 1 
Fineness Ratio 3 Cp 
Shapes VU 1.5 Mo = 2 Vo 3 
New shapes 0.0845 0.0773 0 0686 
Newtonian 0.0877 0.0778 0.0689 
von Karman 0.0895 0.0857 


It is interesting to note that the pressure-drag differences be- 
tween the shapes at a given Mach Number are small Thus, for 
most design purposes, any one of these shapes may be considered 


fora minimum-drag nose. 
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Streamline Curvature Effects for Finite Wings 
in Wind Tunnels 


Alan Pope 

Supervisor, Experimental Aerodynamics 
Corporation, Albuquerque, N.M. 

September, 1953 


Division, The Sandia 


SUMMARY 


A method of applying the Lotz r correction to wings is outlined, suitable 


for the case of finite span wings in wind tunnels 


INTRODUCTION 


ee CASE OF STREAMLINE CURVATURE EFFECTS for two- 
dimensional models has been discussed by Glauert! and Allen 
and Vincenti.2 The case for large models of finite span has 
been covered by Sivells and Deters? and Swanson and Toll. 
However, the methods presented in the last two references re- 
quire several hours to solve, and a simpler correction of adequate 


accuracy seems worth while. 


GENERAL THEORY 


The basis for a streamline curvature effect is found in the vari- 
ation (along the chord) of the upwash induced by the tunnel 
boundaries. Numerical values of this variation for reasonable 
wings in reasonable wind tunnels indicates that the angle induced 
by the upwash varies almost linearly along the chord, so that to 
all intents and purposes the wing in a wind tunnel is in a circu- 
larly bent flow, or, conversely, the wing has had added to it the 
loading on a circular arc airfoil whose curvature and angle of at- 
tack corresponds to the boundary induced curvature of the flow 

The circular are airfoil’ has the following characteristics of in- 
terest to this specific problem: First, its loading is composed of a 
flat plate loading and an elliptic loading. The portion that is a 
flat plate loading is proportional to the angle of attack between 
The elliptic 


the ‘‘horizontal” chord line and the relative wind. 
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loading is proportional to the amount of camber, which in turn 
defines the angle of zero lift. Since the angle of zero lift is essen- 
tially the slope of the circular arc airfoil at its three-quarter point, 
the elliptic correction is defined by the difference between the 
boundary induced upwash at the half-chord (where the slope is 
parallel to the chord) and at the three-quarter chord. The flat 
plate correction, which has no moment about the quarter-chord 
point, is defined by the difference in boundary induced upwash 
at the quarter-chord and that at the half-chord. The upwash af 
the quarter-chord is covered in the customarv wall connections 
written in the form 


Aa = 6 S/C) Cr, (57.3 
ACp z &S/C) C 


where 
S = wing area 
C = tunnel cross-sectional area 
6 = factor based on model span and tunnel shape 


Since the variation of boundary induced angle is linear, the 
variation between quarter- and half-chord is the same as that 
between half- and three-quarter chord. Hence, once we find the 
(flat plate) angle of attack correction for streamline curvature 


Aa,-, the additive (elliptical loading) lift correction will be 


ACL a = —Aa,,a 
where a = slope of the wing lift curve (either measured or theo 
retical is satisfactory), and the additive moment correction 


(since the moment of an elliptic load is zero about the half-chord 


1s 
ACmic/4)sc = —0.25AC, 


From the above corrections, we see that the finite wing in a closed 
tunnel acts like a wing with excessive camber at an excessive 
angle of attack. The additive corrections increase the measured 
angle of attack to account for the angle change and reduce the 
lift and negative moment 

The method for obtaining the numerical values for the stream- 
line curvature effect is that given by Lotz,® who expresses the 
downwash aft of the wing lifting line (assumed at the quarter 


chord ) according to 
additional boundary induced upwash aft of quarter-chord line 


boundary induced upwash at quarter-chord line 
The value of the additive angular correction is then 
Aaw. = tAa 


Values of 7 are given in Figs. 1 and 2 (replotted from references 
6 and 9) for a wide number of tunnel cross sections and model 


spans. The “tail length’ used for wings should be one-fourth of 


the M.A.C 
The size of Ae 


(Additional values in reference 9 

is usually negligible (around 0.03C, 
For reflection plane models or low aspect 
becomes | irge 


for wings 


of high aspect ratio. 
ratio models, it can rise to 0.8C, or more and 


enough for inclusion 
IFFECT OF COMPRESSIBILITY 
The effect of compressibility on the wing streamline curvature 


to increase the effective chord by 1/yV1 — M?. 
In many 


correction is 
This increased length should be used in determining rt. 


cases the streamline curvature correction at near sonic speeds 


can exceed 0.5 
Example: 
= ().85 for a wall mounted semispan panel model 


speed and VE = 
having the following dimensions: Panel span 7.00 ft., wing area 


32.67 sq.ft., M.A.C. = 2.46 ft., taper ratio = 0.5, wing aspect 
7- by 10-ft. tunnel 


Find the streamline curvature correction at low 


ratio = 6.0, airfoil = 65-006. 
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1.2 \ | (1) From page 22 of reference 4 the boundary induced upwash v= 
Elliptic test section, * orn at the wing results in an angular correction of 
>= ——_] Sa = 1.078CL 
10 | + (2) From Fig. 2, using an extrapolated value for a tunnel 
| | height-width ratio of 0.35 and a ‘“‘tail length” of 2.46/4 = 0.615, 
; , = 0,615/20 = 307, = 0.25 
Circular test section L/B 0.615/2 0.0307, and rt 0.25, 
| Aas.e = (0.25) (1.078C,) = 0.269C, 
| (3) A similar example using the unit vortex approach (page 
| 22 of reference 4) yields a value of 0.254C,. The difference of 
0.015° (at C_ = 1.0) is not surprising in view of the uncertainties 
| ° 
—+—_—+- pea sales SS of the extrapolation and the fact that reference 4 uses an average 
" value rather than that at the centerline. Better agreement 
| would reasonably be expected for a floor mounted model, sine 
| the tunnel shape would then be more reasonable. 
(4) At M = 0.85, 
+ 
' 1,/BQ. — M?*)'7? = 0.0307/[1 — (0.85)?]'/? = 0.0583 mee 
/ % Tunnel height 
| " “Tunnel width — and r = 0.475. Aa, then equals (0.475) (1.078C_) = 0.513C;. 
= +— a Se a 
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Toil length Ry . ‘ 
PABA ae a. Corrections for Partial-Span Models with Reflection Plan, End Plate or No E 
Tunnel width 8 End Plate in a Closed Circular Wind Tunnel, NACA TR No. 843, 1946 Oo 
4 Swanson, Robert S., and Toll, Thomas A., Jet-Boundary Corrections for 
Reflection Plane Models in Rectangular Wind Tunnels, NACA TR No. 770, 
1943 
5 Pope, Alan, Basic Wing and Airfoil Theory, p. 110; McGraw-Hill Book 
FIGURE 1 Company, Inc., New York, 1951. 
6 Lotz, I., Corrections of Downwash in Wing Tunnels of Circular and Elli phic 
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1.2 7 Silverstein, Abe, and Katzoff, S., Experimental Investigation of Wind 
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1.0 Off-Center Positions of the Wing and to the Downwash at the Tail, NACA TR 
No. 547, 1935, : 
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Alan Pope a 
ec 
Supervisor, Experimental Aerodynamics Division, The Sandia 
. Corporation, Albuquerque, N.M. the 
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T APPEARS THAT KEEPING a wind-tunnel model exactly on the 
tunnel centerline is much more important for transonic testing 
Tunnel height than for subsonic or supersonic testing, since unsymmetric wave 
oe _ ss e ‘ . P 
0.2 Tunnel width reflections may cause serious errors in pitching moment measure- 
ments. The transonic case differs from the supersonic one in 
Span that the structure of the detached shock and its reflection cannot Si 
. be computed, and, hence, the Mach Number range in which 
Tunnel width : ie 
| trouble may be expected can only be obtained empirically. 
re) The <j al “oe a 
0 Ol 02 03 0.4 The situation is as follows: SS Cl 
Consider a body in a transonic wind tunnel at a very low super- Ae 
Teil_tength . fa sonic Mach Number (Fig. 1), and let its location be closer to the | 
Tunnel width 8 o1e , , way 5 t J 
floor than the ceiling. The bow shock wave reflected from the Jai 
floor will strike the body further forward than will the shock 
reflected from the ceiling, and a nose-down pitching moment will I 


FI appear that would not occur in free flight. This extraneous 
FIGURE 2. pitching moment becomes larger when one wave completely ty 
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Fic. 1. Bow shock pattern for model low in wind tunnel. 
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Fic. 2. Distortion of pitching motion curve by asymmetrical 


transonic reflection, model low in tunnel. The strength, shape, 
and extent of the distortion vary for every model-tunnel com- 
bination. 


misses the body (see dotted lines in Fig. 1) largely because the 
asymmetric pressure is then farthest from the center of gravity. 

A moment curve for a finned body of revolution is shown in 
Fig. 2 for one case that resulted in a maximum erroneous trim of 
17° angle of attack. Visual observation methods were finally 
used to determine the region of reflection free data. 

While the region and magnitude of the detached shock reflec- 
tion phenomenon vary for each model-tunnel combination, it 
appears safe to conclude that transonic detached shock re- 
flection errors will not be serious with regard to lift and drag, 
that their pitching moment contribution may be minimized by 
keeping the model on the tunnel centerline, and that the Mach 
range in which they exist must be determined by visual or other 


empirical methods. 


Similar Solutions of Compressible Laminar 
Boundary-Layer Equations 


Clarence B. Cohen 

Aeronautical Research Scientist, Lewis Flight Propulsion Laboratory, 
NACA, Cleveland, Ohio 

January 11, 1954 


7 A RECENT NOTE,! Drs. Li and Nagamatsu treated the com- 
pressible laminar boundary-layer equations having similarity- 


type solutions. It was indicated that four classes of solutions are 
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possible. In the present note transformations are given which 
reduce three of these classes to only one, and the remaining class 
is also discussed. 

In reference 1, the boundary-layer equations resulting from 
application of Stewartson’s transformation? were subjected to the 


similarity requirements: 


vy = Z°V,Pf(n) 
S = S.g(n) (1) 
Y = 92°V,9 
where 

Vi = transformed free-stream velocity 

S = (H/H,) -— 1 

H = stagnation enthalpy 

(Z, Y) = transformed longitudinal and normal coordinates, 


respectively 
This led to the system of ordinary differential equations, 


f’’’ = Kif’? — Koff” — Kil + Sug)) 


2) 
g" + Kefe’ =0 j 
with the boundary conditions 
f(0) = f(0) = 0, g(0) = 1 
n—>o, f'’—-Kse g—>0 3) 


Four classes of similarity resulted, corresponding to four com- 
binations of a, b, p, and g. These classes were distinguished by 
their external-velocity distributions, differential equations, and 
boundary conditions, all involving various parameters, A, 
[For class (3) the given definitions of AK; and A¢ were incomplete. | 


The velocity distributions were as follows: 


Vi = [((1 + 2g) (K3Z + cy 


Class (1): as 

Class (2): Vi = Ce®2 , 
Class (3): Vi; = (Ks/m) Z™ | 
Class (4): V; = CZ* 


The differential equations and boundary conditions of classes 
(1) and (3) can, however, be reduced to those of class (4) by the 


following transformations: 


fs = fi V1 + 2q)K1,1 

Class (1) < = Vl + 29:)K1,1 
K34 = (1 + 2q:)~! } 5) 
\fs = fr( Kes )(273— 1)/2@— pat 1) 


Class (3) he, a n( Ke ,) (24 +1)/2(qs— pst1) 


where the last subscript on any quantity identifies the class in 
which that quantity is defined. 

Class (4) itself can be reduced further through the following 
transformation: 


i =f, V (Kas + 1) 2) 


7 =m V (Asa + aia | (6) 
S = Swks 
and the corresponding external velocity can be written 
V, = C(Z — Zo)”™ (6a) 
The resulting system of equations is 
jr + If" = af? 1-5) 
Ss” + fs’ =0 | (7) 
fo) = f'(0) = 0, S(O) = Se 
f(-) =1, S(o) = 0 


. . +4 . > ~~ ) 
which is identical with the system developed in reference 2. 


Thus, class (2) remains the only distinct similarity type. 

Class (2) is a limiting case of class (4) which, however, cannot 
be obtained by transformation alone from class (1), (3), or (4), 
but arises from considering the special combination: @ = b= 








282 JOURNAL OF THE AERONAUT 
0, p = '/2,q = —'/2 in Eqs. (1). This corresponds to the case 
K3.4— of case (4) and can be reduced to the form of Eqs. (7) 
with 6 = 2.0 by the transformation 
f =f:V Ki/2| 
7 =m VK,/2 (8) 
S = Syg \ 
where A, > O (favorable pressure gradient ) 
For A, O (unfavorable gradient) the same transformation, 
but with (—K,) replacing (A,), gives the same boundary condi- 


tions, but applied to the new system 


fp’ = ff’ — Wf’? —1—-S) (9a) 

SS" + fS’ =0 (9b) 

From an examination of the right member of Eq. (9a) for S = 0 
(the insulated-wall case) in the neighborhood of f’ = 1, it can be 


seen that no solution that does not indefinitely oscillate about 


f’ = 1 is possible. 
duced for the case with heat transfer, which thus may be a dis- 


No such simple conclusions have been de- 


tinct case. 
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The Method of Superposition of Planar Wave 
Systems 


Wallace D. Hayes* 
Office of Naval Research, London, England 
January 7, 1954 


— PURPOSE OF THIS NOTE is to show the solution of a linear- 
ized supersonic flow problem by means of the method of the 
superposition of planar wave systems. The equation for the 


velocity potential is taken in the form 


brz + byy — b = 0 (1) 


after the application of the Prandtl-Glauert transformation and 
z is taken in the downstream direction. The problem of interest 
is to determine the potential at a point in the upper half space 
(for which y > 0), with no incoming waves (null initial condi- 
tions), and with ¢, given on the upper side of the plane y = 0. 
The solution to this problem has been obtained by Clarkson! 


and may be expressed 


I dx dz 
P(X, Yo. 50) = — ,(x, O, y) (vy >0) (2) 
rs Jr R 


where R is given by 


(3) 





R= V(x — 2)? — (% — x)? — 3 
and I is the region on the plane y = O for which (s) — =) > O 
and Risreal. Clarkson used the mathematically elegant method 
of Riesz; the result may also be obtained by other methods. The 
method of this note is based on the concept that the solution in 
a singularity-free region, under certain reasonable conditions, 
may be represented as the superposition of planar wave systems, 
the superposition consisting of an integration over the wave 
orientation angle 6 as a parameter. In the present example the 
limitation to the upper half space with no incoming waves limits 
6 toa range of 7, and the general solution may be expressed 


* Scientific Liaison Officer. The opinions expressed here are those of the 


writer and do not necessarily reflect those of the U.S. Navy. 
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9(x, y, 5) = f(s — x cos @ — ysin 8, 6) dé 4 


J0 


In terms of the general solution of Eq. (4) the normal derivative 


on the plane y = 0 takes the form 
Tr 
g(x, 0,5) = — f'(s — x cos 6, 8) sin 6 dé 5 
0 
and the integral appearing in Eq. (2) may be expressed 


i [{ [{. axa: ae 
_ do, = F(x0, Ye, 2, 0) dO 6 
T I R a 0 ; 
i : a sin 6 
F= f'(s — x cos 0,0 dx dz 7 
T e I R 


The transformation is now made 


where 


cos 6 (Sa 


—o = (% — 3) — (% — x 
V £7 — yo? sin? 6 cos Y = (x) — x) — (39 — 2) cos 6 (8b 


for which the functional determinate may be shown to be 


. ' , 
6 sin @ 

J ( ) - 9 
x, 2 R 


The region I corresponds to the inequalities 
OSyYsSrm —-e7 <o Ss —ypsind (10 


The function F may now be evaluated 


1 *—y sin 6 * 
F= | f( — x, cos 0 + ¢,0) dy dt 
T 
- 0 
F f 


= f(t» — xX) cos 8 — yo sin 6, 6 (11) 
whence the desired result Eq. (2) follows immediately from 
Eqs. (4) and (6). 

In case yy = 0, a much simpler transformation may be used, in 


which the variable z is retained and y is defined by 
cos Y = (xy — x)/(29 — 2) (12) 


In the ensuing analysis the relation 


7 


7 dy 7 


0 1 — cos @cos y¥ sin 6 


is used. 
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Comments on ‘‘Ripple-Type Buckling of 
Sandwich Columns’’! 


Paul Seide 
Research Engineer, Structural Research Group, Northrop Aircraft, 
Inc., Hawthorne, Calif. 


December 21, 1953 


N A RECENT PAPER,! Eringen concludes that wrinkling of sand- 
wich columns will take place in neither the symmetric nor the 
antisymmetric modes found by other writers but in an asymmet- 
ric mode. The author also emphasizes that the solution of a 
two-dimensional mixed boundary value problem of elasticity for 


a rectangular orthotropic medium, as reported in his Doctoral 
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Dissertation,” is necessary. The writer has only recently had 
the opportunity to read this paper in detail and has come to the 
conclusion that neither of the author's statements is correct for 
the problem solved. In addition, an effect that may be impor- 
tant has been neglected. The following derivation will indicate 


the defects of the author’s presentation. 


In order for the problem under consideration to be in the form 
of an eigenvalue problem and to have a reasonably simple solu- 
tion, it must be stated as follows: 


A sandwich column is compressed into a state of uniform axial 
strain (Fig. 1). Prior to buckling, there is no friction between 
the load applying mechanism and the ends of the core so that 
movement of the ends perpendicular to the middle surface of the 
core is permitted. After buckling, no further movement of th« 
ends perpendicular to the middle surface of the core is permitted, 
and the axial load applied to the structure remains constant. A 
similar device was used in reference 3 to simplify the problem 
Other formulations of the problem will result in 
The faces of the 
sandwich are of equal thickness and of the same material. It is 
assumed that they may be analyzed by the usual methods of ele- 
The core is an orthotropic elastic me- 


treated there. 
formidable complications of the final results. 


mentary beam theory. 
dium. 

If the destabilizing effect of axial forces in the core is neglected, 
stresses and displacements in the core are subject to the following 


boundary conditions: 


(0, y) = o(L, y) = 0 (1) 
o(0, y) = oA L, y) = 0 (2) 
l . . Ox, &/2) 0*0(x, h/2) t Orz,(x, h/2) 
E yt? + Gert : - + 
12 ox4 Ox? 2 Ox 
h 
a,\ x, = 0 (3) 
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Sandwich column prior to buckling. 
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Fic. 2. Positive directions of core stresses 
1 . O4v(x, — h/2) O7n(x, — h/2 t Or,,(x, —h/2) 
Egt tT Gert — _— 
12 ox‘ Ox? 2 Ox 
h 
oy| x, — = () (4) 
9 
t O°x(x, h/2) l : h Ou(x, h/2 e 
2 a. S Try X, dx — = f) 5) 
2 Ox? Ey Jo 2 Ox 
2 » 9 wa 
t O*%x(x, —h/2) 1 
- + — Tr(x, —h/2)dx 4 
2 Ox? Est Jo 
Ou(x, —h/2 
0 6 
Ox 


where u, v, o:, ¢,, and 7r,, are displacements and stresses in the 
core, in addition to those existing prior to buckling (Figs. 1 and 
2) 

Conditions (1) and (2) for the ends of the sandwich column 
indicate the stated assumptions that no further displacements 
perpendicular to the middle surface of the core are permitted to 
take place and that the axial load.applied to the column remains 
constant. These are identical with Eringen’s boundary condi 
tions for the ends. 

Conditions (3) and (4) imply continuity of the faces and core 
in the direction perpendicular to the middle surface of the core 
An error has been made in reference 1 in the formulation of these 
conditions, as will be shown, and leads to the deduction of an 
asymmetric mode of wrinkling. In Eqs. (6) of reference 1, which 
correspond: to Eqs. (4) and (5) of the present note, both o,(x, 0 
[= o,(x, —h/2) of the present note] and o,(x, t) [= o,(x, 4/2) of 
the present note] have the same sign. However a positive value 
of o, acts downward on the upper face and upward on the lower face 
so that the faces are bent in opposite directions (see Figs. 2 and 3) 
With v(x, h/2) and v(x, —h/2) both having the same positive 
direction, it is readily seen that o,(x, h/2) and o,(x, —h/2) [or 
o,(x, 0) and o,(x, t)] should have opposite signs. Positive values 
of the shear stress, on the other hand, act in opposite directions, 
but one is applied to the bottom of one face while the other acts on 
the top of the other face so that their effect is the same. 

Conditions (5) and (6) imply continuity of the faces and core in 
the direction parallel to the middle surface of the core and are 
similar in concept to Eq. (A10) of reference 3. The first term in 
each represents the face strain due to bending, at the juncture 
of the face and core; the second term represents the face middle 
surface strain due to the shear stresses transmitted between the 


core and faces: the third term represents the core strain at the 
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The stress-strain relations in the orthotropic core are given by 

} st, Ou/OX = Aynor — Ayo, / 5 

~~ A” Ov/Oy = —Ay0r + Annoy > 10) 
(du/Ody) + (dv/Ox) = (1/pry) rey) 


if vis hn where da, @i2, and dz. are constants depending on the materia] 
( ” 


x 








properties and on whether the core is taken to be in a state of 








o.,t plane stress or in a state of plane strain [see Eqs. (1.3.3) and 
6, (4,2 y ‘Tay (X h) (1.3.4), reference 4, page 7]. It should be noted that Ey in Eqs, 
z 3 (3) and (4) is either E or E/(1 — v?), depending on whether the 
Try (x,-B) state is one of plane stress or plane strain. . 
Equilibrium of stresses in the core is satisfied by introducing 





a,0$)> f » t, Airy’s stress function ¢, such that 
<7 


o, = 0*y/dy" / 


a. a 
Y-— a, = 0*¢/Ox (11) 
a + —— v( > Zz aa ao oaialilial 


cor 











Compatibility of strains is satisfied if 


Fic. 3. Forces acting on top and bottom faces. dy dy Wy 
— 2K = () (12) 
ox‘ Ox?dn? on 
juncture of the core and faces. The sign difference in these two where 
equations is due to strains being taken at the bottom of one face 
and top of the other. In most previous work, the first two terms 
have been neglected, and the conditions have been written as 


n= ey 
4/ 

€ = WV Q22/ay 

K = (1/2V anaes) ((1/psy) — 2ar2! 


u(x, h/2) or (Ou/Odx) (x, h/2) = 0 (7) 
A solution of Eq. (12) which satisfies Eqs. (1), (2), and (9) is 
u(x, —h/2) or (0u/dx) (x, —h/2) = O (8) given by 


Additional conditions that are specified are that the ends of the gy = (A sinh ny + B sinh my + C cosh ny + 
faces are simply supported—that is, D sinh vy) sin (max/L) (18) 


h h h h where 
aro. ) ov(o, ~ ‘) ov, ) ox 1, - ) : 
2 2 2 2 me (Eth 4 z=!) 
— = — = — = = ( y= | 
\ 


oe 9 9 
~ 


ox? Ox? Ox* ov? L 


(9) mre ( kK +1 Ix - ') 
— oe wight 
L\NV 2 V2 


The solution of the problem as next given is similar to that of 


reference 4. from which 
ma\? ; . 7 _ mrx 
melee 3 (A sinh ny + B sinh wy + C cosh ny + D cosh ry) sin 
mr ro ; max 
ty = 7 (1A cosh ny + v2 B cosh wy + »,C sinh ny + » D sinh ry) cos 7 
rl } (14) 
ou j mra\? ; 2 mra\? . ] _ maxx 
= 4] 7a + ay | (A sinh ny + Ccosh ny) + | m%an + [- dy |(B sinh vxy + D cosh vey)? sin — 
ax (t L L f L 
oo /L Jace 27012 ‘L)*a22 : . MmaAx 
ad [? ay + (mx/L)ae (Aen oy + Coah ns) + v"aiz + (ma = (8 cosh wy + D sinh “| aia - - | 
Vv) Vo } 


The substitution of Eqs. (14) into Eqs. (5) and (6) and the addition and subtraction of these equations yield 


j{ a\? mah\? Lh ht (mah\? [n(h/2) Pars + (mah/2L)2a22 | a 
AW o) mF Vor) Flop 27 aor n(h/2) water Nd 
J h\? mah \* LAhioh t [mah\? [vo(h/2) ai. + (marh/2L)*a22 
ss at oie © Mel SO es Oe 4 vs(h/2) 
J h\? mrah\?2 1h kh t {mah\? [v(h/2)}ai + (mrh/2L)?ax hl 
di sk aie) Raed es ee ni(h/2) a 


J h\? mh \? = aa t mah\? [ve(h/2)]2a12 + (mah 2L)%a22 h) , 
ol. dead Siete |. is Oe 9 yo(h/2) a ee ee 


The substitution of Eqs. (14) into Eqs. (3) and (4) and the addition and subtraction of these equations yield 


coth vy, st =Q (15a) 
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: mah \? 
| t { 2] 1 mr w(h/2)]}*ay2 4+ 1 
{<] y coth » + - Ert - P 2] coth » + 
| nh “ o h 12 I * “ 
vi(h/2) 
mah 
tt h 2 i mr \? [ve(h/2)]*%ai2 4+ h Ree 
B«1 + v coth v + Ey — 21 coth » = 0 (l6a 
h*2 2" Al 12 I 2 
v(h/2 
h : mah\? } 
| t oh h oF 4 mr \? nu\> aye + = a h 
Cé1 4 v tanh yp + Eyt —PpP = 2L tanh »; 
| h 2 2 h{_12 # 2 
v(h/2 


h ; mah \? ] 
| t h h 2 ] ma \? v a t : ) a h 
D<1 4 v , tanh v + Eyt ) - Pf 2 2L tanh » a l6b 


| .- 2 2 h{12 


It is seen that the equations separate into two (sets—Egqs. 15a) and (16a) containing unknowns A and B and corresponding to an 
antisymmetrical wrinkling mode, and Eqs. (15b) and (16b) containing only the unknowns C and D and corresponding to a symmet 


rical wrinkling mode. These equations would not have separated with the same sign in front of o,(x, h/2) and o,(x, —h/2) in Eqs 
3) and (4). The wrinkling criterions are obtained by setting the determinant of the coefficients of A and B and of C and D in thes« 
equations equal to zero (to permit values of A and B or of Cand D other than zero Then 
l h h l h nh 
l+-~yp coth v 1 + y coth pv 
h 2 2 h 2 2 
h h 
a, + y; coth v : as + y2 coth v ; 
l {mak\? _  /t\? lh 2 2 " 
Org = Ey T (Iva) 
3 2 h re h h 
B, coth », B» coth v 
9 
h h 
a, + y; coth v as + y. coth » 
9 
t h h t h h 
1 + v tanh p 1 + i tanh vy 
h 2 2 h 2 2 
h 
a, + yy, tanh yy a» + ye tanh v 
1 {mah i Lh 2 < - 
Ors =: Ey + (17b 
3 rt i h oe h j 
8, tanh v Bs tanh v 
h h 
a + y, tanh v : ae + ye tanh pv 
where 
a = |v, 2(h/2)]?a,, + (mah/2L)?a, 
3 = }[m, o0(h/2)]?ay. + (mah/2L)2a.2} 
m1, (h/2 

1 = (1/2E,) (h/t). o~(h/2) — (t/h) (mxh/2L)?,, 

Eqs. (17a) and (17b) are readily applicable to those sandwich structures for which K > 1 For A -1, both » and vy. are pure 
imaginary quantities, and the equations must be revised by converting the hyperbolic functions of a pure imaginary number into their 
trigonometric equivalents. For |A| < 1, » and v2 are complex, and the complex functions in the right sides of Eqs. (17a) and (17b) 
must be transformed into their real equivalents. For |K| = 1, the criterions can be obtained by finding the limit of Eqs. (17) as K 
—~ 1. In all cases the resulting values of o.- are real. In calculating values of o.,, the right sides of Eqs. (17a) and (17b) must be 
minimized with respect to the quantity mah/2L, which is present in »(h/2), vo(h/2), a1, a2, Bi, 82, y: and ye, for given values of the 
material constants and geometric properties of the sandwich structure For a finite value of L, it is necessary to take m equal to 
the integers 1,2,3,.... However, since the wrinkling wave length is small compared with the length Z of the column, it is customary 


to consider L to be infinite and to consider the quantity L/m to be a wave length A, which may vary continuously. The minimization 
in this case is probably best done numerically than analytically. However, it might be possible to solve the equation for the analytical 
minimum for the quantity //t, thus obtaining corresponding values of o., and h/t for given values of the material constants and the 
buckle wave-length parameter mah/2L. 

Comparisons may be made with the results of other analyses of the problem. Eqs. (17) are presumably similar to those of refer 
ences 5 and 6, which investigate the problem with identical assumptions, according to reference 4 rhe solution of reference 4 is ex 


pressible as 


1 {mxh\? E $\3 lh (1/a,) — (1/a 
ag * | " uf : 
3 ya h 2t (B:/e) coth »,(h/2) — (B2/a2) coth v(h/2 
1S 
1 f/ mah\? E t \? lh (1/a;) — (1/ae 
Cecrs = Lf . a . 
3 at h 2¢t (B:/a;) tanh »(h/2) — (82/ae) tanh w(h/2) 
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In this solution the effect of transmitted shear stresses on bending of the faces is neglected, and relative displacements between th 


core and faces due to stretching and bending of the faces are permitted 


The corrected solution of reference 1 is obtained by placing y 


and 2 equal to zero in Eqs. (17), in which case the actual location of the juncture of the core and faces is partially taken into account 


by considering the effect of shear stresses on bending of the faces but allowing relative displacements between the core and face 


Thus, the problem stated in reference 1 has both the symmetric 
and antisymmetric wrinkling modes found in other investigations 
and is readily solvable in terms of simple functions. It should 
be noted that wrinkling takes place in sinusoidal waves as im- 
plied by Eringen’s analysis but of which a false impression may 
have been given by Fig. 1 of his paper. For other boundary con 
ditions for the curve and faces, the elaborate analysis made by 
Eringen of the boundary value problem with mixed boundary 
conditions would be necessary, but it is doubtful that an asym- 
metric buckle pattern would result so long as both faces were of 
equal thickness, of the same material, and subject to the same 


boundary conditions. 
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Spin Effect on Base Pressure of Cone-Cylinders 


R. Lehnert and S. M. Hastings 
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T° THE WIND-TUNNEL TESTING of spin-stabilized missiles, the 
effect of model rotation on its aerodynamic characteristics 
A few systematic experiments with ro- 
Simi 


usually is disregarded. 
tating models have been carried out at subsonic speeds. ? 
lar data at supersonic speeds were preliminary in character.* 

Experiments were recently carried out at the Naval Ordnance 
Laboratory to determine the spin effect at zero angle of attack 
on the stability of the laminar boundary layer and on base pres- 
sure at supersonic speeds (Magnus effects are outside the scope 
of this investigation). A family of five cone-cylinders all having 
25° apex angles and 2-in. diameters, varying in total length from 
L/D = 6 to 10, were tested at a Mach Number of 2.86 with 
spin rates varying from 0 to 43,000 r.p.m. The maximum r.p.m. 
corresponded to a value of approximately '/, for the ratio of cir- 
cumferential velocity to translational velocity. Spark-schlieren 
pictures were taken in order to determine the character of the 
boundary layer and wake flow. It has been shown previously 
that base pressure depends on the status of the boundary layer 
at the base,* ® and therefore some correlation between the effect 
of spin on boundary-layer transition and on base pressure was 
expected. 

Experimental values of the ratio of base pressure to free-stream 
static pressure are plotted versus Reynolds Number for various 
spin rates on Fig. 1 (Reynolds Numbers are based on free-stream 
conditions and the length, along the surface, from the cone tip to 
the cylinder trailing edge). From spark-schlicren photographs, 
it was found that the minima of the curves were obtained when 
When the 


boundary layer at the base is laminar (two short models, L/D = 


boundary-layer transition occurred at the model base 


6 and 7), increasing the model spin rate produces a considerable 
decrease in the base-pressure ratio. Spark-schlieren photographs 
show that the decrease in base pressure was accompanied by an 
upstream movement of transition in the wake mixing zone with 
increasing spin rate. This decrease becomes less pronounced 
when the boundary layer at the base becomes less stable (in- 
creasing body length), and the trend reverses when the boundary 
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layer at the trailing edge has become fully turbulent, but the de- 
pendency is less pronounced. 
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The minima of the base-pressure ratio versus Reynolds Number 
urves occur at lower Reynolds Numbers with increasing spin 
rate. In addition, increasing the model spin rate from 0 to 43,000 
rpm. shifted the average transition location on the model up- 
stream by approximately one model diameter, corresponding to a 
change in Reynolds Number of about 0.5 * 10°. 
shift of the transition location due to spin can be seen in the photo 


The upstream 


graphs of Fig. 2, which are given as an example 


Consequently, 
we conclude that spin destabilizes the laminar boundary layer, 
thereby promoting transition. 

It was feared that vibrations induced by the spinning of the 
model might have caused a destabilization of the laminar bound 
ary layer. The model vibrations incurred during the spin test 
were therefore measured and simulated in frequency and ampli- 
tude Base 


transition location were then measured on the nonspinning model 


with an electromagnetic vibrator. pressure and 


und found to be independent of model vibration 
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Buckling Loads for Columns of Variable Cross 
Section 


B. E. Gatewood 
Mechanics Department, USAF Institute of Technology, 
Patterson Air Force Base, Ohio 


December 24, 1953 


Wright- 


if iemeonguners HAS DISCUSSED buckling of bars of variable cross 
section and has given tables for buckling coefficients for cer- 
tain taper ratios and for certain types of cross section. This note 
presents curves for the buckling coefficient for all taper ratios 
and for any moment of inertia variation between constant and 


the sixth power. From Fig. 1, the taper ratio is defined as 


a = a/(a +l) (1) 


and the parameter m for moment of inertia (J) variation is de- 
fined by 


I = I,(x/a)™ (2) 
where for solid section m = r + 3s with x’ the width variation 


the thickness variation. These variations are restricted 
to give minimum J about same axis for entire length. 


With the coordinate system of Fig. 1, the differential equation 


and x 


for the deflection y is 


Pa™ 
x"y" + gy=0, g? = (3 
El, 
This can be written as Bessel’s equation 
tv” + tr’ + (2 — n*?)v = 0 (4) 
wher« 
I 
y= Vx, t = 2nqx™ nu n= (5) 
. m— 2 
rhus 
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y = Vx [(AJ,(t) + BY, At (6 
in which J-,(¢) may replace J),(t) if m is not a whole number 
For m = 2, 

y = Vx(A exp jilg? — (1/4)]'? log (x/a)} 4 
B exp | —ilg? — (1/4 log (x/a)} 
With boundary conditions y = 0 at x = aand y’ = Oatx = 


a + /, the determinant for calculation of the buckling coefficient 


can be written as 


J 2pn) Vnti (2Qpn a/?") — Int (2pn a/?*™) Y,(2pn 0 (7) 
where J—,, may replace J’, (m not integer) and 
p? = Pa?/EI, = q? a 1’) (2 


For given m and a, p can be obtained from Eq. (7 The results 


are shown in Fig. 1 for all a from 0 to 1 and m from 0 to 6, where 


wr? REI, x? sas 
P.,. = J k = p*( l1—a)a™ “_ solid curves 
} i? } 
(9 
RET 
Pig = (22*) ae Qr2*k = p?, dashed curves 
a* 


Fig. 2 shows a cross plot of k against m with a as parameter 
Note that, for any m which makes m a multiple of one-half 


1), Eqs. (6 


For m = 


(such as 4 and (7) reduce to trigonometric 


expressions 4, Eqs. (6) and (7) are 


y = x[A cos (q/x) + B sin (q¢/x)]) 


ai 


Boyd? obtained this result of Eq. (10) by direct change of vari 


(10) 
tan p(l — a) = 


able and made tests on steel pin-ended columns (both ends ta 
pered) form = 4,3,and1. The test results agree closely with the 


theoretical buckling load 
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Fic. 1. Variation of buckling coefficient with taper ratio a. 
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Fic. 2. Variation of buckling coefficient with parameter m 


For design with given load, given length, and given type sec- 
tion, an interesting question is what value of the taper ratio a 
gives the column of least weight. If A = Aji(x a)", then the 
average area is 

Ay 1 — alt t 


0 = (11) 
(h + 1l)a’ l—a 


whence, from Eq. (2), the moment of inertia corresponding to 


: m/h 
| (12) 
— gy 


(13) 


this average area is 


Ao m/h l1— al” 
TI => I; = In 
A; (h+1)(1 


Take klz = Kh, Eq. (9), whence 


P , \m/h 
- (h+1)(1 — a)”/" 
K=f 

1 ae a! +1 

Now m and hare related through the width and thickness equa- 
tions, x” and cx’ of the column. For solid columns hk = r 


+ s,m =r-+ 3s, and since r > 0, s > 0, Fig. 1, 


1<m/h<3 (14) 


Fig. 3 shows values of K, Eq. (13), in regions of peak values for 
For given r and s, or m and h, read maxi- 
then Eqs. (9), (11), 


A, may be restricted 


m/h = 
mum XK and corresponding a from Fig. 3; 
(12), and (13) give the column dimensions. 
by P/A, < yield stress of material; r and s are restricted in that 
For 


1, 2, and 3. 


minimum J should be about same axis over entire length. 
constant load and constant length, the weight of the tapered 
column to the weight of the constant cross-section column is 


W/W. = K7*™ (15) 


Fic. 3. Maximum values of K for minimum weight. 


For round thin-walled columns, m/h = 3, and Eq. (15) checks 


Feigen’s*® results on tapered thin-walled columns for m = 


1m 
3, and 4, 


Recently, several authors‘~? have considered stepped col- 


umns. In some cases the stepped column may be approximated 
by a continuous tapered column such as discussed above. A 
curve can be put through three points on the stepped column to 
determine a and m, whence & can be read from Fig. 1 or 2. For 
instance, in the example in Ram and Rao’s paper,’ use the three 
points J,;/J; = 1 at the wall, J2./J; = 0.81 at 3L/8 from wall, and 
I,/I, = 0.36 at the tip to get a = 0.17, m = 0.56, whence k = 
0.82 or (9r?/4)k = their result is 2.019. Also, columns 
with different taper ratios in several regions may be approxi- 


2.02; 
mated by one taper ratio. 
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